Prom Hadrons to Nuclei: Crossing the Border 



Silas R. Beane", Paulo F. Bedaque\ Wick C. Haxton'' 
Daniel R. Phillips'"^, and Martin J. Savage" '' 

° Department of Physics, 
University of Washington, Seattle, WA 98195-1560 

^ Institute for Nuclear Theory, 
University of Washington, Seattle, WA 98195-1560 

Department of Physics, 
Ohio University, Athens, OH 45701 

Jefferson Laboratory, 
12000 Jefferson Avenue, Newport News, VA 23606 

sbeane, bedaque, phillips, savage, haxton@phys.washington.edu 

The study of nuclei predates by many years the theory of quantum chromody- 
namics. More recently, effective field theories have been used in nuclear physics to 
"cross the border" from QCD to a nuclear theory. We are now entering the second 
decade of efforts to develop a perturbative theory of nuclear interactions using 
effective field theory. This chapter describes the current status of these efforts. 



1 Introduction 

Nuclei are fascinating objects. In a nucleus the protons and neutrons, col- 
lectively known as nucleons, are bound together by the strong nuclear force. 
At a fundamental level these interactions are described by Quantum Chromo- 
dynamics (QCD), a theory of quarks and gluons carrying color charges that 
are asymptotically-free at short-distances. However, the quarks and gluons 
in a nucleus are very far from being asymptotically-free. Instead they com- 
prise individual, colorless, nucleons, which largely retain their identity in the 
many-body system. The color-singlet nucleons are then bound to each other 
by what can be thought of as "residual" QCD strong interactions. This sketch 
of nuclear dynamics from the QCD point of view — brief as it is — makes it 
clear that from this standpoint the nucleus is an incredibly complicated, non- 
perturbative, quantum-field-theoretic, infinite-body problem. 

This makes it truly remarkable that the properties and interactions of 
nuclei can be largely understood within the following framework, which we 
shall refer to here as the "traditional" picture of a nucleus. The A nucleons 
in a nucleus are understood to be non-relativistic particles interacting via a 
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quantum-mechanical Hamiltonian consisting of two-body potentials and three- 
body potentials which have considerably smaller, but still important, effects: 

N , |2 

^ = + YlVNN{r^,r,) + J2 V^N{r^,r,,rk). (1) 

i—l i<j i<j<k 

This Hamiltonian is then used in the Schrodinger equation 

7J|Vj(')) = (2) 

to determine the nuclear wave functions IV'^*-') and energy- levels E^'K With 
modern computing power the energy-levels of this Hamiltonian can be calcu- 
lated with ^ 1% precision for nuclei as large as A = 8 and soon calculations 
for twelve-body nuclei such as ^^C should become available. The results of this 
program are in impressive agreement with experiment. With wave functions 
in hand, properties of the nuclides such as magnetic and quadrupole moments, 
electromagnetic form factors, and the response of light nuclei to weak probes 
such as neutrinps can all be calculated. For a thorough review of these suc- 
cesses see Ref. El. 

However, the success of a picture described by the Hamiltonian in Eq. 
raises a number of other questions, among which are: 

• Why are nucleons the appropriate degrees of freedom inside a nucleus? 

• Why are two-body interactions Vnn dominant? And how should they 
be described? 

• How do we construct a three-nucleon interaction V^n? 

• Are all contributions to observables equally large, or is there some set of 
small parameters which can be used to organize the calculation? 

All of these issues are part of one overarching question, namely: 

Can we both understand how the Hamiltonian in Eq. (|l|) emerges 
from QCD to dominate the dynamics of multi-nucleon systems, and 
systematically refine it to accurately describe all nuclear observ- 
ables. 

One kind of answer to this question would be provided by computations 
of the properties and reactions of nuclei directly from lattice QCD. For in- 
stance, consider what happens when a proton captures a neutron and forms 
a deuteron, in the process emitting a photon: the radiative capture process 
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np — > ^7- A lattice calculation of np d-y would provide an unambiguous 
cross-section in terms of the quark masses and Aqcd- However, such a cal- 
culation would be very difficult to perform. For instance, a cartoon of one 
contribution to np — > c?7 in terms of perturbative quarks and gluons is shown 
in Fig. |l|. Sources of quarks and gluons that have non-zero overlap with the 
proton, neutron, and deuteron as well as a source for the photon would be 
used to generate the reaction amplitude. One key problem with using lattice 
QCD to calculate this reaction is that a significant amount of effort goes into 
forming the asymptotic states themselves — before one even begins computing 
their interactions. Indeed, since the deuteron has such a small binding energy 
and so is quite extended compared to the nucleon, considerable effort will be 
required to generate the deuteron on the lattice. The deuteron is yet to be 
formed in any lattice calculation, and so, to put it mildly, a comprehensive 
lattice study of this simple hadronic process is not imminent. Therefore, given 
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Figure 1: A "cartoon" of a contribution to np —> d-y in QCD. 

the current state of affairs in lattice QCD, serious efforts to determine proper- 
ties of multi-nucleon systems must start with neutrons and protons as the basic 
degrees of freedom. This has the considerable advantage that one is working 
in terms of the asymptotic states which are those detected in the laboratory. 

Clearly an important first step in describing the dynamics of a multi- 
nucleon system is to describe the interaction between just two-nucleons, i.e. the 
two-nucleon scattering amplitude. This has been accomplished very success- 
fully by fitting the nucleon- nucleon (NN) potential Vn n ■ There are many NN 
potentials on the market, but the most accurate models are now tuned to fit the 
over 3000 data which constitute the database for NN scattering below labora- 
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tory energies of 350 MeV ^ In all of these models the long-distance part of the 
interaction, corresponding to inter- nucleon distances of order 2 f m l/rriT^, is 
described by the model-independent exchange of one pion — a sixty-five year old 
idea, due to YukawaB. At shorter distances ~ l/(2m7r) some potentials include 
the exchange of two pions. This too can be derived in a model-independent way 
from chiral perturbation theory (xPT). In contrast, once one looks at the force 
between nucleons separated by ^ 1 fm, there are a variety of short-range com- 
ponents present in the different NN potentials currently in use. Some assert 
that the short-range interaction arises from the exchange of vector and scalar 
mesons, while others look to the naive quark model for inspiration. More pop- 
ular, and more effective^ than either of these models, ace potentials which fit the 
data essentially to within experimental accuracy □'BeI, using a short-distance 
part of Vmn that is essentially a best fit, with no dynamical assumptions other 
than continuity, differentiability, and the desire to include all possible opera- 
tor structures which can arise in the NN system at low energies. This part 
of Vmn occurs in the region where the nucleons themselves have substantial 
overlap, and r ^ A^^ with A^ is the chiral-symmetry breaking scale. It is 
clear that without input from lattice QCD, assertions about the underlying 
dynamics in this region are nothing more than guess-work. Thus the use of 
these "best-fit" short-distance NN interactions is very much in the spirit of 
effective field theory. However, once we move beyond the NN interaction and 
start considering three-nucleon forces and the construction of current opera- 
tors for use in multi-nucleon systems the picture becomes significantly more 
complex. As nuclear systems are inherently non-perturbative it is difficult to 
see which parts of a computation are essential to the description of a particular 
observable, and which parts are not. Effective field theory (EFT) allows one to 
both systematize calculations in light nuclei, and simplify the picture of these 
nuclei down to its essential ingredient^. 

In this chapter of the Handbook of QCD we will describe the progress that 
has been made in constructing and applying EFT's to multi-nucleon systems, 
and the interactions of these systems with "external" probes such as 7's, W^'s, 
Z^'s, and TT. Since Weinberg's pioneering effortslZl'DEl in the early 1990s to set up 
a systematic description of multi-nucleon interactions using EFT much effort 
has been gone into realizing this goal. At present these nuclear EFTs are far 
from being able to determine, for instance, the B{E2) values for transitions 

"In Ref.S this was chosen as the uppermost energy for which data would be included in 
the NN database because the pion-production threshold, at Tij^j, = 280 MeV, begins to play 
a role there, and so the NN collisions are no longer elastic. 

''There is an old saying: Those who do not know history are doomed to repeat it. However, 
there is a slightly more modern version that may be more relevant: Those who do not 
understand history will not recognize progress. 
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in ^^Mg, but important progress has been made in understanding systems of 
two and three nucleons as well as in discussion of probes of the two-nucleons 
system. We will also discuss the progress and problems facing EFT's attempts 
to discuss larger nuclei. 

An EFT is defined by the most general Lagrange density consistent with 
the global and local symmetries of the underlying theory, a declaration of the 
regularization and rcnormalization scheme, and a well-defined power counting 
in terms of small expansion parameters. This issue is somewhat controversial 
in the presence of bound- or quasi-bound states because there are differing 
criteria for what constitutes a rigorous EFT. On the one hand there is cutoff 
EFT. In cutoff EFT the power counting is essentially dimensional analysis at 
the level of the potential: operators are ordered according to inverse powers 
of the momentum space cutoff. Since the quantum mechanical potential is 
not an observable quantity, the justification for the power counting is neces- 
sarily a posteriori. One computes to a given order in the expansion and then 
one performs an error analysis for observables such as phase shifts in order to 
verify that cutoff dependence is removed order-by-order in the expansion of 
the potential. This method has much in common with lattice gauge theory 
with improved lattice actionslij. Lattice errors are removed order-by-order in 
the lattice spacing by including operators in the action with increasing powers 
of the lattice spacing. On the other hand, there is EFT as it is tradition- 
ally practiced in a purely perturbative setting, like chiral perturbation theory 
(xPT). The power counting is in the scattering amplitude itself and there- 
fore applies directly to observable quantities. In this case one can explicitly 
prove rcnormalization scale independence at each order in the EFT expansion 
of observables. Moreover, one can make an a priori estimate of the theoretical 
errors which arise from working to a given order in the perturbative expansion. 
The distinction between these two criteria becomes an issue only when there is 
a bound- or quasi-bound state near threshold. In these cases at least one local 
(as in a contact operator) or non-local operator (as in a particle exchange) must 
be summed to all orders and so the EFT is intrinsically nonperturbative. The 
nonperturbative nature of the NN interaction leads to several unresolved for- 
mal issues in a cutoff EFT of NN interactions. These issues will be discussed 
below. 

The underlying motivation for using EFT to describe systems with well- 
separated lengths scales can be found in the Chapters by Leutwyler, by Meifiner 
and by Manohar, and we will not repeat it here. Instead, we begin our tour 
of EFTs in nuclear physics by looking at the simplest nuclear system, A = 2, 
at momenta well below the pion mass. At these very low energies the only 
relevant degrees of freedom are the nucleons themselves. We will show how to 
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systematically describe both NN scattering and the coupling of electroweak 
probes to the NN system. These calculations of very-low-energy properties 
provide analytic results for all observables. They must be compared, process- 
by-process, against the results of the "traditional" approach. Subsequently, we 
move away from the A = 2 system and discuss EFT for the ^ = 3 system — also 
at very low energies. The EFT treatment of the three-nucleon problem provides 
some truly novel features, including the emergence of a type of renormalization- 
group (RG) behavior not seen before in a single-parameter system. 

We then turn to higher energies, where the pion must be included as 
an explicit degree of freedom in the Lagrangian. Attempts to include pious 
in perturbation theory are presented and their current status discussed. We 
then discuss Weinberg's original suggestion for dealing with the pions non- 
perturbatively. This is to apply chiral perturbation theory (^PT), as formu- 
lated in the Chapters by Leutwyler and Meifiner, to the NN potential and 
other two-nucleon irreducible operators. We display some successes of this ap- 
proach, and discuss some unresolved formal issues. Finally, we discuss recent 
progress in incorporating the ideas of EFT into more traditional methods of 
calculating nuclear properties. 

2 Two-Nucleon Systems at Very Low Energies 

Ideas akin to EFT appeared in nuclear physics as early as the 1940s, when 
Bethe intmduced Effective Range Theory (ERT) as a description of NN scat- 
tering EllEj. Bethe showed that at low energy the experimental data did not 
depend on dynamical details, but only on a few numbers which summarized 
the impact of strong interactions on observables. ERT is an attempt to de- 
scribe both elastic and inelastic processes involving two nucleons, provided 
that all momenta involved are less than the pion mass. By construction, 
strong-interaction data is reproduced in ERT to arbitrary precision in this 
low-momentum regime. However, ERT fails to converge to the measured am- 
plitudes for processes involving photons or other "external" probes. As the 
energy of experiments increased, ERT become less relevant, and more compli- 
cated descriptions in the form of potential models, such as the Hamiltonian 
(|l|), became popular. In this section we will show how to derive an EFT which 
will reproduce ERT, and hence the NN data, to arbitrary precision, provided 
that p < m7r/2. EFT goes beyond ERT, as it systematically includes operators 
which enter when one calculates the response of the NN system to external 
probes. We will give explicit examples of analytic calculations for low-energy 
reactions which can be carried out to ~ 1% accuracy. 
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2.1 Naive Dimensional Analysis Really is Naive 



What scales might be reasonably expected in the two-nucleon sector if one 
was simply presented with the QCD Lagrangian? Neglecting quark masses, 
QCD is a theory with only one length scale, Aqcd- In the hadronic sector 
this scale appears as the scale of chiral symmetry breaking, ~ 1 GeV, and 
the nucleon mass M. At low energies, the effective theory of QCD will be 
one whose matter content includes only non-relativistic nucleons and massless 
pions. Neglecting factors of An the only mass scale in the theory is A^ ^ 

The inclusion of the quark masses in the QCD Lagrangian explicitly breaks 
chiral symmetry. The pion acquires a mass « 140 MeV, and so introduces 
another "small" length scale. What then, is a reasonable scale to expect for 
the binding energy of nuclei? Naively one might expect that the nuclear size 
will be set by the pion mass. However, one must remember that the pion is 
derivatively coupled, and so it is not clear that the pionic dynamics are strong 
enough to be responsible for bound-state formation. Even accepting that the 
pion is the particle whose exchange generates the non-perturbative interaction 
which leads to the formation of the deuteron, one would not then expect it 
to be bound by only B = 2.224644 ± 0.000034 MeV. This corresponds to 
a pole in the NN scattering amplitude at momentum |k| = ij — with 7 = 
y/MB = 45.7025 MeV — in the ^Si channel '|: significantly less than our prior 
expectation. In the ^Sq channel the situation is even more striking. There 
is no bound-state in the ^5'o channel, but a pole on the second energy-sheet 
represents a quasi-bound state. In terms of momentum, this pole sits at |k| ^ 
(— i)(8 MeV). The unnaturally low energies at which this and the deuteron 
pole reside are an experimental fact. Such threshold states can only result from 
a fine-tuning between long and short distance physics, and so their existence 
can never be understood using naive dimensional analysis (NDA). 

In an EFT such fine-tuned bound states result in operators that are naively 
irrelevant (in the RG sense), playing a crucial role. Consequently, power count- 
ing based on the NDA we are familiar with in the zero- and one-nucleon sec- 
torscJ does not apply. In fact, as is discussed below, the two-nucleon system is 
very close to a non-trivial fixed point in the RG-flow. Another way to see this 
is to remind ourselves of something obvious: the theory is nonperturbative. 
Nuclei are bound states of nucleons and so will never arise at any finite order 
in perturbation theory. Unlike the baryon-number zero and one sectors, the 



'^Strictly speaking this is not true, since the mass difference between the nucleon and the 
A is Ma - M = 232 MeV. 

"^Here we use the nuclear spectroscopic notation ^^+^Lj, where S is the total spin of 
the A'^A'' system, J denotes its total angular-momentum, and L denotes its orbital angular 
momentum. 
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EFT appropriate to nuclei must admit non-perturbative structure somehow 
and somewhere. Thus, an infinite number of loop-diagrams must enter at the 
same order in the EFT expansion. 



2.2 Low-Energy NN Scattering Data and Unusually Large Length Scales 

Over the years there have been literally hundreds of experiments performed 
to measure NN cross-sections and polarization observables below the pion- 
production threshold. In 1993 the Nijmegen group performed a phase-shift 
analysis (PSA) of these experiments □. This analysis summarizes all of this 
data in terms of the phase shifts and mixing parameters that describe scat- 
tering in the different spin-angular momentum channels. Clearly any theory 
which reproduces these phase shifts will also reproduce the data. Conversely, 
although the mapping from data to phase shifts is not mathematically unique, 
the Nijmegen group's imposition of "reasonable" physical constraints appar- 
ently makes it quitij difficult to obtain phase shifts which differ radically from 
those found in Ref.El. We will always compare the EFT phase shifts and mixing 
parameter(s) with the Nijmegen PSA and not directly with data. 

At low energies the NN cross-section is dominated by S-wave scattering. 
Let us recall the behavior of the phase shift in a strongly-interacting system. 
Firstly, the S-matrix of two non-relativistic nucleons scattering in the ^5*0 chan- 
nel with center-of-mass momentum |k| = Afriab/2 is 

,2, 2^\k\ .\k\M 

where 6 is the phase shift, and Tjab is the kinetic energy of the incoming nucleon 
in the laboratory frame. We use a normalization such that explicit computa- 
tion of Feynman diagrams, in an arbitrary frame, gives the non-relativistic 
amplitude A. From elementary scattering theory, |k| cot 6 is an analytic func- 
tion of external kinetic energy with a radius of convergence bounded by the 
nearest t-channel singularity. In the case of NN scattering this is the pion cut, 
which occurs at |k| = Within this radius of convergence we may use 

the expansion: 

|k|cot<5 = + irolkp + ri|k|4 + .... (4) 

a 2 

While Eq. corresponds to expanding about |k| = 0, |k| cot J can, of course, 
be expanded about any point within the region of convergence. In general, the 
size of the coefficients ro (the effective range) and ri (the shape parameter) 
is determined by the range of the interaction between the nucleons, which is 
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~ l/mTrlij'cJ. One might also expect the scattermg length a to be "natural" 
in this sense. However, a is not constrained by the range of the underlying 
interaction, and, as we have already discussed, the scattering lengths for NN 
scattering in both the ^5*0 and ^Si channel are significantly larger than I/tott- 
In the channel the empirical values of the parameters a, rg, and ri are^ 

a ^ -23.714 fm, tq = 2.73 fm, ri = -0.48 fm^. (5) 

We could also expand about |k| = in the spin-triplet S-wave channel, 
the '^Si. However, since the deuteron is present in this partial wave it is 
convenient to-jjostead make the expansion of |k| cot (5 about the bound-state 
pole positionBEJ. In fact, if one truncates the expansion in Eq. at order 
|kp, it is straightforward to relate the location of the pole, |k| — i-f, to the 
scattering length and effective range. 



1 

^ ro 

Neglecting relativistic effects, and expanding around the deuteron pole, |k| = 
ij one obtains: 

|k| cot 6 = ~j + \pa (|kp + 7^) + W2 (|kp + 7^)2 + . . . , (7) 
withlli: 

7 = 45.7025 MeV = (4.318946 fm)~^ , 

Pd = 1.764 fm, W2 = 0.389 fm^. (8) 

Finally, we should mention that since the NN interaction is not spher- 
ically symmetric it mixes states of different orbital angular momentum. In 
particular, the deuteron is not solely a bound state in the '^S'l channel, but 
instead represents a bound state in the coupled ^Si — ^Di channels. Below we 
shall show that this mixing effect is sub-leading in the low-energy NN EFT, 
but still numerically significant. A more detailed discussion of this mixing is 
presented below. 

"^Throughout, we quote only the strong-interaction part of the phase shift and give the 
parameters in the np channel. There are corrections for charge-independence breaking and 
charge-symmetry breaking which shift the value of both a and ro if one considers nn or pp 
scattering. 



l-\ l- 



2ro 



(6) 
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2.3 You Can't Choose Your Neighbors: Living near a Fixed Point 

We now begin our attempts to build an EFT which reproduces this low-energy 
NN data. The location of the deuteron pole and other low-energy NN 6b- 
servables are input for the EFT, and so it might seem that such an EFT can 
accomplish very little. However, as explained above, a low-energy pole in the 
scattering amplitude necessarily involves short-distance physics, and so the 
EFT must have this informatimi as input. It simply cannot be expected to 
yield it as output. (See Ref. EZI for an early and explicit implementation of 
this idea.) Putting the bound state in "by hand" in this way does not vitiate 
the ability to make predictions with the EFT. Once the physics of the low- 
energy pole is included in the theory we can use symmetries to predict other 
observables, such as the cross-section for np d'y. 

The extreme situation where 7 = corresponds to a bound state at thresh- 
old. If, in addition, the mass scale of the potential is much greater than the 
incident center-of-mass momentum the theory is approximately scale invari- 
ant. This leads to a somewhat peculiar EFT, which exhibits a non-trivial fixed 
point. 

Our goal is to recover such an NN scattering amplitude from a Lagrange 
density with local four-nucleon momentum-independent operators alone: 

iV - Co {N^PNy (n'^PN) , (9) 

where P is a spin and isospin projector. Note that in the normal xPT counting 
the contribution of the term is suppressed. However, in a non-relativistic 
two-body bound state its contribution is of the same order as the contribution 
from the iDo term in Eq. (Q). In other words, the kinetic energy of the NN 
bound state is the same order as both its potential energy and its total energy. 
Further, the size of the interaction term can be estimated by NDA, 
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2M 



where A is the "high" scale in this theory, which is ^ m^r. 

The Lagrange density in Eq. (^) is sufficiently simple that the sum of 
all Feynman diagrams contributing to NN scattering, shown in Fig. |^, can 
be found exactly. Note that in constructing the NN amplitude in this way, 
contributions from an irrelevant operator are being resummed. If the size of Cq 
were in fact given by NDA, as in Eq. (|lO|), each additional loop in the bubble 
sum would be suppressed by one power of fc/A. Thus, without any further 
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work we sec that the NDA scahng of Eq. (^0|) cannot produce a low-energy 
NN bound state. 



x-xx-:xro<-- 



Figure 2: Feynman diagrams eontributing to the scattering amplitude from the Lagrange 

density in Eq. (0). 

If, m spite of this, we proceed to sum the bubble chain depicted in Fig. || 
we obtain an amplitude: 



A 



— [ Co + CqIqCq + CqIqCoIqCo + 



-1 



(11) 



where Iq is the momentum-independent single-loop integral, which in the NN 
center-of-mass has the form: 
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(12) 



with |k| = \JME the three-momentum of each nucleon in the center-of-mass 
frame. Note that this /o is linearly divergent, and so we have evaluated it in 
n space-time dimensions. In order to go further, a subtraction scheme must 
be specified; different subtraction schemes amount to a reshuffling between 
contributions from the vertices and contributions from the the UV part of the 
loop integration. The usual MS scheme corresponds to taking the n — > 4 limit 
of the expression in Eq. (|l^), obtaining: 



(13) 



In this case each successive term in the bubble-sum in Eq. (g_l|) involves the 
product CoM|k|, and so, in order to produce a bound state of momentum 7 it 
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is necessary to have the coefBcient Cq scale as Co ^ l/{M'y). This unnatural 
scaling in MS, and^the problems it leads to, were first discussed by Kaplan, 
Savage, and WiseEa. Co is simply unnaturally large: in fact it is ill-defined for 
the case of a threshold bound state. 

There were many discussions as to how to deal with this problem in 
EFT, with some advocating the use of cutoff regularization, instead of di- 
mensi on al rfigularization , so that the, iinnflituxallj large value of Co would not 
arise EjMeIeI. A number of papers oBa^y^EEE] discussed the peculiar fea- 
tures which arise in this EFT when one naively implements different regular- 
ization schemes. Ultimately though, the issue here was not associated with the 
use of a specific regularization scheme. Instead, ilic-jsplutipa to this difficulty 
lay in the development of a new power counting c3'Ej'L!j't!iJ'E3. Subsequently it 
was shown how to implement this power counting on a diagram-by-diagram 
basis, .using a general subtraction procedure called Power Divergence Subtrac- 
tionBEl (PDS). 

In PDS both the n — 4 and n = 3 poles are subtracted from each bubble. 
The logarithmic divergence in n = 3 — which corresponds to the power-law 
divergence in n = 4 — is then included in the expression for the bubble. PDS 
thereby keeps track of linear divergences (see Ref.Ea for extensions of this idea). 
This allows for a fine-tuning between the coefficient Co and the linear diver- 
gence which produces a shallow two-body bound state. In PDS, the integral 
/q in Eq. (O) is defined to be 



lo ^ -(^^j(/. + z|k|). (14) 

A similar expressicui|is,abtained with a momentum cut-off0 or with a momentum- 
space subtraction EjEZlEa. If PDS is employed the bubble sum (O) yields 



M 



47r 1 

M ^+/.+z|kr ^''^ 

Matching to the leading-order ERT amplitude in the '^Si channel, given by 
Eq. (^ with |k| cot 5 ~ —7, we find 

47r 

In the limit of vanishing binding energy, it is then clear that 

fi^lfi Coifi)] = 0. (17) 
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At thfi fixed point the theory is confomiaUy invariant (in the two-body sec- 
tor)E3. Therefore, in order for the four-nucleon interaction in Eq. (|^) to gen- 
erate a system with a shallow bound state (i.e. 7 '--^ 0), Co must be near a 
non-trivial ixed point, as given by Eq. (|l7|). The existence of this non-trivial 
fixed pointcJ is indeDcndent of the particular regularization and renormaliza- 
tion scheme chosen EH. 

The explicit introduction of the subtraction scale /i into the expression for 
the amplitude results in Co having its scale set by fi and not by 7. Ultimately 
all of the couplings associated with higher-dimensional two-nucleon operators 
also scale with /i, as well as with the underlying scale, A. A useful power- 
counting emerges if fi is taken to be of order the extcijiitd mnjpfj]it.iuTT,-Lk| . This 
leads to a regulator-independent power-countingESl^J'llatlJBilliilEjtlJ'M, known 
as Q-counting. In Q-counting 7 and the external momentum k Q are both 
denoted by Q, and Feynman diagrams are analyzed in terms of the power 
of this parameter that they carry. The result is an expansion in powers of 
Q/A, where A is the underlying scale which sets the size of the higher-order 
coefficients in the effective-range expansion, i.e. A 1/re ~ m,r- Q-counting 
is regulator-independent, but if PDS is used then each diagram has a definite 
order in the Q-expansion and gauge invariance is explicitly preserved. These 
advantages greathi. simplify higher-order calculations of processes other than 
NN scattering BEI. 

This completes the discussion of the Lagrangian in Eq. (|) in the NN 
system. Now we generalize this analysis to the most general S-wave NN in- 
teraction that is consistent with Galilean invariance. In this case the sum of 
all graphs contributing to NN NN can be determined exactly if PDS is 
used to define the divergent loops. Suppose the tree- level amplitude for NN 
scattering is 

Arcc = -(f) ""E^2.(/.) \kr. (18) 



r=0 



Here the coefficients C2r(/^) are linear combinations of couplings in the La- 
grangian for operators with 2r gradients which contribute to S-wave scatter- 
ing. If this vertex — instead of the vertex Co — is used for each interaction in 
the bubble chain of Fig. || the resultant amplitude is 

The yU-independence of this amplitude gives RG flow equations for each of the 
coefficients C2r(M). Matching to the ERT amplitude, Eqs. (||) and (0), for NN 

^Consequently fi is also denoted by Q. 
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scattering, requires that for fi ^ j ^ Q the couphngs C2r(M) scale as 



Thus the rth term in the sum for Ati-ca (^2^^^' , has a Q-counting order of Q^~^. 
The counting rules for this low-energy effective theory of NN scattering are 
then as follows: 

• Each nucleon propagator scales as l/Q^; 

• Each loop scales as Q^; 

• Each vertex C2rk^^ scales as Q^~^. 

In this way each NN loop contributes one factor of Q, but the vertex Cq 
scales as and so, by construction, the leading-order contribution in this 

EFT, is the complete bubble-chain shown in Fig. |[ with Co operators at all 
the four-nucleon vertices. The leading-order amplitude, A'--~^\ is then given by 
Eq. (p^). Contributions to the NN amplitude which scale as higher powers of 
Q come from perturbative insertions of derivative interactions, C2rk^^ , f > 1, 
but these are always "dressed" by the leading-order (LO) amplitude A^~^\ 
since the product of .A^"^-* with a loop is of order 

2.4 Low-Energy Nucleon- Nucleon Interactions: ERT EFT(7f) 

The previous section contains all that is needed in order to set up the EFT 
for low-energy NN interactions. We will call this theory EFT(7f) since it 
does not contain pious as explicit degrees of freedom. EFT(7f) reproduces the 
ERT NN scattering amplitude A order-by-order in a Q-expansion. If this was 
all that it did, then EFT(7f) would represent no more than ERT. However, 
EFT(7f)'s power becomes apparent when external probes of the NN system 
are considered. One of the very interesting results obtained in EFT(7f) is that 
ERT generally gives the LO contribution to observables where the NN sys- 
tem is coupletLtXL low-momentum external probes, such as photons and weak 
gauge bosonsE3o. This justifies the use of ERT to estimate, for instance, the 
magnetic moment of the deuteron. However, ERT computations of deuteron 
properties are not always successful, since ERT only gives the LO contribution 
to these observables, and subleading corrections can be numerically significant. 
Operators in EFT(7f) exist which are absent in ERT, and these operators con- 
tribute to electroweak observables. Including such operators allows EFT(7f) 
to produce closed form, analytic expressions for electroweak processes such as 
np — > d'y and vd I'd which are accurate at the 1% level. 
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2. 5 NN Scattering in the ^Sq Channel 



Nucleon-nuclcon scattering in the ^5*0 channel is described by the Lorentz 
invariant one-body Lagrange density 



_ -^0 

2M 2M 

and the two-body Lagrange density — up to N-^LO: 



(21) 



Co 



'n'^P^nV (N'^OIN) + h.c. 



" ^4 2 

where is the spin-isospin projector for this channeh [J ~ 0, 1 — 1] 



V8 



0'2 T2T 



Tr 



cob 



(22) 



(23) 



and where the derivative operators are 



^2 



1 r- 

4: 

1 r- 
16 



p"e' - ifep^'S' + 6B'p"6' - 45'p''0 + & V 



(24) 



The subscript on the coefficient denotes the number of derivatives in the 

operator, and D = D D . We have employed the covariant derivative, 
Dp, = {do + ieQAo, —V -I- ieQA), in defining the two and four derivative op- 
erators in Eq. ( p4| ) where is the electromagnetic field and Q is the electro- 
magnetic charge matrix. Contributions to the Lagrange density in Eq. (22) 
including more than four derivatives can be found in Refs.llj'cll. Hereafter rel- 
ativistic corrections to observables will not be included in expressions, unless 
explicitly stated otherwise. It is always straightforward to include the thkd 
term in Eq. (pl|), which generates these corrections, in perturbation theory Ej. 

At LO the scattering amplitude found from Eqs. (21) and (|2^) is that 
shown in Eq. (|l|). Matching this to the amplitude from ERT, Eq. (|), 
gives 



th C^'^o^Ai) = 



(25) 
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Note that in Q-counting tq ^ 1/A, while 1/a, k ^ Q. Using this counting and 
expanding the ERT amphtude defined by (j^) and in powers of Q/ A we 
can match to the NLO amphtude in EFT(7f). This yields: 

= (26) 
Carrying out the same procedure at N^LO, the matching gives: 

It is interesting to note that onhr.jthe sum can be constrained 

by NN scattering. RG scalingtillS^I of the operators in Eq. ( p2[ ) indicates that 

while the contribution from is N^LO, the contribution from IS 

N'^LO. The time-ordered product of two ''Cj '^"^ operators does not induce the 

momentum structure of the operator. However, inelastic processes can 

be sensitive to different combinations of ^Cl^°^ and 'fc['^°\ 

The convergence of EFT(7f) can be tested by an examination of the be- 
havior of the phase shift (5*^ order by order in the Q-expansion. The results 
are shown in Fig. |^. It is clear from Fig. ^ that in the ^5*0 channel the ex- 
pansion must be carried out to N*LO or higher to describe the strong inter- 
actions for |k|< 100 MeV. Also, it is clear that EFT(7f) will not be valid for 
|k|> 150 MeV. At these momenta the full |k|cot(5 description of the phase 
shift starts to deviate from the observed phase shift. Given the location of the 
t-channel pion-exchange cut in the NN amplitude this deviation is expected. 



, (27) 



2.6 NN Scattering in the Coupled ^Si — '^Di Channel 

Scattering in the ,7=1 channel is somewhat more complicated because both 
L = {^Si) and L = 2 {'^Di) states contribute. The S-matrix describing 
scattering in the coupled channel J = 1 system is written as 

a _ f e'2*°cos2£i ie*('5«+'52)sin2ei\ , , 

- [ie'(So+S2) s[n2ei e'^^^ cos2si ) ' 

wheije-|We use the "barred" parameterization of Ref . 0, which was also used in 
Ref. e3. It will follow naturally from the EFT(7f) that ei is supressed by 
compared with (5q"\ and therefore, up to N'*LO, the ■^Si — element of this 
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Figure 3: The phase shift 5^ as a function of the center of mass momentum |k|. The 
dashed curve corresponds to the LO result, the dotted curve is the result at N'^LO, the thick 
solid curve corresponds to N'^LO. The thin solid curve corresponds to the phase shift from 
ERT, and the dot-dashed curve are the phases of the Nijmegen PSA. 



matrix may be written as: 

Soo = e'^'" = 1 + J' . . (29) 
cot Oo — I 

The phase shift Sq has an expansion in powers of Q, Sq = S^'' + S^^^ + S^'^ + 
where the superscript denotes the order in the Q expansion. 

Now, taking the expressions in Eqs. (^) and (j^), and expanding in powers 
of fc/A and 7/A, we obtain the LO and NLO contributions to A, which are 

^ ^ 1 ,(0) ^ ^ Z,-l 

M 7 + i|k|' M 7 + i|k|' ^ ' 

where Zd is the residue of the NN scattering amphtude at the bound-state 
pole. In this way, the position and residue of the deuteron pole are correct in 
the NLO amplitude. In terms of the parameters 7 and pd, Zd is 

^ (31) 



1 - IPd 



Consequently Zd — 1 is of order Q. Zd is also related to the asymptotic S- 
state normalization of the deuteron wave function by Ag = 2jZd. Using the 
Nijmegen PSA value for As we get Zd = 1.69. 



Forming the logarithm of both sides of Eq. (29), using Eqs. (^ and 



and expanding in powers of Qj it is straightforward to obtain the LO and NLO 
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results 



TT — tan 



4"(|k|) = -^V^|k|, (32) 



7/ ■■ ■■ 2 

which are shown in Fig. ^. The N^LO piece of the phase shift, d^f^\ also 

180 




100 



200 



300 



Figure 4: The phase shift <5o as a function of the center of mass momentum |k|. The 
dashed curve corresponds to , the dotted curve corresponds to 5^'^ + 5^^^' , the sohd_£;urve 



corresponds to 5^ 



(0) 



s(2) 



lidcu 

and the dot-dashed curve is the Nijmegen PSAq. 



depends only on Zd and 7. Up to N^LO, the shape parameter term, W2 does 
not contribute to the S-wave phase shift. In fact, W2 is also numerically small 
(about a factor of 5 smaller than pd), but even if it were not, it enters the Q- 
expansion only at high orders. The convergence of the Q-expansion is clearly 
demonstrated in Fig. ^ 

The Lagrange density describing scattering in the ^5*1 channel is the same 
as that in the ^5*0 channel, but the C"s and the spin-isospin projectors are 
different. It is thus the £1 of Eq. (|2l|) plus 



{N^P'Ny {N'^O^N) + h.c. 



h.c. 



(33) 
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where is now the spin-isospm projector for the [J = 1, / = 0] channel, 

1 



V8 



0'2Cr T2, 



(34) 



and where the derivative operators are 



^ 16 



Matching the amphtudes in Eq. (^| 
by the Lagrangian in Eqs. ( pT]) and 



.(35) 



to the Q-jExpanded amphtude generated 
3|), gives y, up to NLO 



'-'0,-1 
•-"2,-2 



47r 



1 



M(a* 
27r 



7)' 
1 



'-'0,0 



27r 7(Zrf- 1) 
M (/i-7)2 ^ 



M 7(^ - 7)2 



(36) 



Here, the coefficients themselves have been expanded as a power series in Q, 
e.g. t'C^'^'^ = ''c£!\^ ''C0'^ + ''cjf'^ so that the location of the 



deuteron pole is unchanged order-by-order in the Q expansion. 

Turning our attention now to the off-diagonal elements of the S-matrix in 
Eq. (p8|), the El parameter defines the amount of mixing between the '^Si and 
'^Di channels. It has a momentum expansion in terms of the constants E^^"^ 
and e[ 



(4) 



£l 



El 



(2) 



E- 



(4) 



+ 



(37) 



are fit to the Nijmegen PSA □, and are found to be E^ 

(4) _ 



(2) 



) and e['^ 

0.386 fm^ and E'i^' — —2.800 fm'*. As far as the power counting is concerned 

it is important to note that the e['^^ 's arfi, set by physics at the scale A and so 
contributions of order must vanishcj. Thus the first term here is 0{Q), 
and the second 0{Q^). Corrections to Eq. ( ^7| ) are suppressed by two further 
powers of Q. The LO and N^LO contributions to ei are shown in Fig. ^. 

In order to calculate ei up to N^LO we only require one insertion of -4*^°"^, 
dressed with the leading-order amplitude for — scattering (^0|). The 
Lagrange density describing the first two orders of S-D mixing interactions iscSl 



£(sd) _ _j- 



(sd) 
ij,xy 
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Figure 5: The mixing parameter ei as a function of the center of mass momentum |k|. The 



dashed curve corresponds to el^ and the dotted curve corresponds to e^"' + ej 



.(0) 



1) ^ 



determined in Eq. (|37[). The dot-dashed curve is the Nijmcgen PSA 



h.c, 



(38) 



where 



ij,xy 



71-1 



4 V 



(39) 



The tree-level amplitude for an ^Si ^Di transition resulting from this La- 
grange density is 



A tree 



(40) 



Like the S-wave coefficients, the coefficients that appear in Eq. (|3^) have their 
own expansion in powers of Q, e.g. ^C^^'^'' = '^c'^^f'i + ^Cq'^q^ + ■■■■ 

It is convenient to introduce the quantity rjsd- This is the asymptotic D/S 
ratio of the deuteron, and may be extracted from the ^Si — ^Di scattering 
information by analytic continuation of the amplitude to the deuteron pole. 
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In fact, 



Vsd 



tan (2£i 



sin(5o - h) 



(41) 



|k|=i7 



The value is then found from a continuatkin of the Nijmegen PSA into the 
unphysical region: 77,^ = 0.02543 ± O.OOOOtH. 

Matching the ^Si — ^Di scattering amphtude to the Q-expansion of the 
off-diagonal elements of the S-matrix in Eq. ( p8| ) evaluated at the deuteron 
pole, gives the coefficients in Eq. (Esl) in terms of {r/sd, e[^^}, 



6V27r 



Vsd ^^..2(„_^A ' "^^0,0^ - Vsd {Zd 1) ^,J^(^„^)2 



^0,-1 - ^72(^-7) 
^C^% + ^C^% = (Z,-l) 



3V27r 



M-i^{^i - 7)2 



(42) 



In this way one ensures that that all of the asymptotic physics of the deuteron 
state: its binding energy, and its asymptotic S and D-state normalizations, are 
reproduced in EFT(7f) at low orders in the Q-expansion. 

Finally, the four-nucleon operator that first contributes to scattering in the 
'^Di channel involves four derivatives. It is also renormalized by the operators 
contributing to ei . Analysis of the powers of Q which appear in the ^Di phase 
shift leads us to conclude that it contributes to the observables we examine 
in this chapter at a much higher order than we will be concerned with, and 
therefore we do not consider it further. 



2. 7 Processes Involving the Deuteron: Formalism 

When considering clastic or inelastic processes involving the deuteron it is 
convenient to define an interpolating field for the deuteron I?i(x) = PiN {x) , 
where Pi is the projector defined in Eq. (|3^). The Green's function associated 
with this interpolating field is then just that for the channel: 

dx e-^(^*-''"') (0|T \v\{x)Vj{()) \ |0) = G{E)5,,, (43) 



where, by Lorentz invariance, the propagator only depends on the energy in 
the center of mass frame, namely E = E — + . . . with the ellipses referring 
to relativistic corrections. 

The Green's function G is straightforwardly related to the amplitude A in 
the channel. However, some relationships between observables are clearer 
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when expressed in terms of this deuteron-deuteron Green's function. Since 
there is a pole in this channel we may always write G{E) as: 



G(B) = .|<^|. (44, 

where B is the deuteron binding energy and the numerator Z is a. function 
which is smooth near the deuteron pole. When evaluated for E k, B, Z gives 
the wavefunction renormalization Z^. This is related to the Zd defined in 
Eq. (pl|) by simple factors: 



Z{-B) = Z^, 



dE 



87r7 

E=-B 



It is now useful to define "Co-irreducible" Green's functions, since observ- 
ables can be computed solely in terms of these quantities. A Co -irreducible 
Green's function is the sum of all possible diagrams which do not involve the 
leading-order interaction ^C^ and contribute to a given n-point function. 
The irreducible deuteron-deuteron function is denoted by E. It has the expan- 
sion shown in Fig. ^ 

^^^^^ = ^^^^l^y^ 



Figure 6: The perturbative expansion of S. The first row has the LO result, where (gi 
represents an insertion of the deuteron interpolating field. The second row is the NLO 

contribution, arising from a single insertion of the '^Cqq^' and ^C2^2 operators. The 
third row is the N^LO contribution, arising from two insertions of the 

operators, and a single insertion of the ^C^l'f'^ '^C^^,-! ^ '^^i^^-i ^'^^ '^C'i^'-s'- 

A simple consequence of this definition is that the relationship between 
the full C and E is: 

^ ^ /(3g,) ■ (46) 
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Figure 7: The expansion of of the full 2-point function G in terms of the irreducible 2-point 

function S. 



This can also be seen graphically in Fig. |^. 

Recall that we wish to choose the coeiBcients 
so as to reproduce the Green's function (^4|), with the function Z having the 
precise value (|45| ) at the deuteron pole. This imposes the following conditions 
on E: 



T.{E = -B) 



•-"0,-1 



1 dS 



(E = -B) 



(47) 



These conditions hold at each order of the Q-expansion, and so restrict the 
number of independent higher-order couplings which appear in the deuteron 
channel. 

We are now in position to compute the matrix element of an electroweak 
current between two deuteron states. We first define the three-point function 



G^(i?,£^',q) 



■IT 



I?J(x)J,^^(0)I?,(y)J |0) , 
(48) 

where = {E' — i?, k' — k) is the momentum transferred to the deuteron 
system. is related to the desired form factor via the LSZ formula 



(k^J|Je^^„|k,^) = G-\E)G-\E)G%{E,E 



E. E -^-B 



(49) 



where G{E) is defined in Eq. (43). As promised, this formula may now be 
reexpressed in terms of S and the "Co-irreducible" 3-point function, F^. Again, 
it is then a simple matter of definition that: 



GtAE.E.ci) 



r^AE,E ,q)G{E)G{E) 



) 

This relationship is depicted in Fig. ||. Making use of Eqs. (|50|), (^), (| 



(50) 



and (|4S|), S-matrix elements can then be found in terms of the Co-irreducible 
vertices of the theory 



{k',j\J!^jk,i) 



ri;(^g,q) ' 

di:iE)/dE 



(51) 



E, E 



23 



Figure 8: The expansion of full three-point function in terms of the irreducible two- 

and three— point functions S, r**. 



Since E and T both have a simple perturbative expansion in terms of Feynman 
graphs the right-hand side of this equation can be straightforwardly calculated 
to an arbitrary order in the Q-expansion, thereby providing a result for the 
current matrix element appearing on the left-hand side. 



3 Phenomenological Applications of EFT(/r) 
3.1 Electron- Deuteron Scattering 

A deuteron with four- momentum and polarization vector is described by 
the state |p, e), where the polarization vector satisfies p^^e'^ = 0. It is convenient 
to choose the basis polarization vectors so that in the deuteron rest frame 
ef = (Sf. We then write |p, i) = |p, ef), and these states have the normalization 
condition (p', j|p, «) — (27r)'^(5'^(p — p')'5jj- In terms of these states |p, i) the 
nonrelativistic expansion of the matrix element of the electromagnetic current 
is 



(p',jVem|P:*> 

(p',j|JLIp,*> 



1 



2Mj 



^qI?') q.q. 



q^% 



2Md 



2Mrl 



[Fc{q^)5^,{v + p')' + FM{q^) (<5|q. - 5fq 



+^^s('z^)(q^q 



-q^% ) (P + p') 



(52) 



Relativistic corrections to these formulae enter at N^LO. In Eq. (^) = 
Qq — |qp, is the square of the four-momentum transfer. Momentum conservation 
implies q = p' — p. The dimensionless form factors defined in Eq. ( |5^ ) are 
normalized such that 



i^c(O) = 1, 



2Md 



Fm{Q) = -^FQiO) = iiQ, (53) 



where — 0.85741 is the deutepjipjpagnetic moment, and /ig = 0.2859 fm^ 



is the deuteron quadrupole moment 
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As an example, the calculation of the electric-quadrupole form-factor of 
the deuteron will be outlined. In EFT(7f) the physics of this form factor at LO 
is exactly the same as in ERT. The photon couples to the D-state component of 
the deuteron wave function, which in EFT(7f-) appears through a perturbative 

insertion of ''Cq^!^\. Some of the NLO corrections, those that serve to give the 
deuteron its correct asymptotic S-state normalization, could be computed in 
ERT. However, there is also a contribution beyond ERT, depicted in diagram 
(d) of Fig. 1^, where a quadrupole photon mediates an S-to-S transition in the 
deuteron. 




(d) 



Figure 9: Some LO and NLO contributions to the deuteron quadrupole form factor. Dia- 
grams of the form of (a) contribute at LO and higher. Diagrams of the form of (b) and (c) 
contribute at NLO and higher. At NLO there is a contribution from a local counterterm, 

diagram (d). 



At LO, the contribution from diagram (a) in Fig. ^ is, after performing 
the relative-energy-integrations , 

d"-^kd"-H (k + §)-(k + §)''-^|k+f|2 
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f d" 
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J (27r)3 
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(2^)3 (27r)3 [|1 + §|2 -t- 72] [|k + f P 7'] [|kP + 7'] 

gab 



n — 1 

1 _ 1 

2 n-l 



M«cn {^) (^) 



_2 

n-l 
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|k|2+72 [|k+a|2+^2] [|k|2+72] 
gab 



n — 1 



m3 
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where q is the three-momentum transfer to the deuteron, and n is the num- 
ber of space-time dimensions. The integrals appearing in Eq. (|5j) have been 
regulated with PDS. There are no n — 3 poles in the k integration in Eq. (54) 
because the integral is convergent. Including the hermitean conjugate of di- 
agram (a) and wavefunction renormalization the LO quadrupole form factor 
is 



2V2|k|3 



-47|k| + (I672 + 3|k|2) tan-i ( ^ 



(55) 



At NLO, there are contributions from all of the diagrams in Fig. |. They 
are straightfoward to compute, and we will not present the calculatiOjU-here. 
The contribution from a four-nucleon-one-quadrupole-photon operator dia- 
gram (d) in Fig. ^, has a coefficient ^Cq defined by the Lagrange density 



C 



2 ^N^P^Nf (N^PjN) 



1 

n- 1 



n — 1 



-5'^ V • E 



,(56) 



where E is the electric field operator. 

At N^LO the only correction comes from a modification to the one-body 
charge operator. This modification introduces the finite charge-radius of the 
nucleon into the calculation. Thus, the deuteron quadrupole form factor in 
EFT(7f) is, up to N^LO: 



2V27|qP 
4\/27|q| 



4|q| 7 



47|q| + (3|q|' 



(3|qp + 167^) tan-i 



167^) tan"^ 



M 

47 



47 



(57) 



where r^fi is the isoscalar nucleon charge radius. Note that we have elimi- 
nated the bare coupling '^Cq in favor of the RG- invariant quantity SjjLQ. The 
value of 5^Q is chosen so as to reproduce the observed deuteron quadrupole 
moment, — Fq{Q)/M'^ — S^q + J^^^ . The value that does this is 6fiQ — 

-0.0492 fm^ 
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Figure 10: The deuteron quadrupole form factor. The dashed, dotted and soUd curves 
correspond to the LO, NLO and N^LO predictions. The dot-dashed curvereprresponds to a 
calculation with the Bonn-B potential in the formulation ofcJ. 



Once the quadrupole moment is fixed we have a prediction for the 
dependence of Fq. The form factor computed at LO, NLO and N^LO is 
shown in Fig. afong with a potential-modeLcalculation performed with the 
Bonn-B potential in the formulation of Ref. cZl. The second and third terms 
in Eq. (^7|) would be reproduced in ERT, as long as '^Si — ^Di mixing and 
the contribution from the finite size of the nucleon were included in ERT. The 
first term in Eq. (|57| ) is a contribution beyond ERT, and corresponds to an 
insertion of the lowest deuteron quadrupole counterterm not constrained by 
NN scattering alone. The presence of this counterterm at NLO indicates that 
fiQ is sensitive to short-distance physics — a fact confirmed by results gleaned 
from potential-model calculations, and calculations in EFTs with pions (see 
below). The absence of the appropriate short-distance physics in potential 
model calculations is presumably|-.the reason that the best of these calculations 
underpredict /ig by about 5 %tM. 

In contrast, we can define the deuteron's quadrupole "radius", via: 

At N^LO this quantity is given by: 

('i> - + - 1-398 f-- (59) 

independent of the quadrupole counterterm. (r^) is, not surprisingly, less sen- 
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sitive to short-distance physics than /ig. Therefore, we expect that potential- 
model calculations of this quantity will agree substantially better with observa- 
tion and with each other, than potential-model calculations of the quadrupole 
moment will. 

Similar calculations can be performed for the charge and the magnetic 
form factors. At LO, NLO and N^LO the contributions to the charge form 
factor are 



F^^\\C\) = Stan-f^) , 4^^(|q|) = (Z.-l)(F<^?)(|q|) - 1 



4'^(|q|) = -g(4,o> |qP4°^(|q|) ■ (60) 

The sum of the LO, NLO and N^LO contributions is exactly what one obtains 
in ERT. At the next order, N'^LO, there is a contribution from a four-nucleon- 
one-photon interaction that represents physics beyond ERT. The deuteron 
charge radius that arises from the sum of the form factors in Eq. (^) is, 
including the leading relativistic effectEj, 



(^?)^"^ = (-^.o> + ^ + ^2 - 4.565 fm^ (61) 

which is to be compared with the experimental value (r^) — 4.538 fm C3n 
The EFT(7f) calculation agrees with data (and with potential-model calcula- 
tions too) to within ~ i%. 

The magnetic form factor of the deuteron is dominated by contributions 
from the single nucleon magnetic moments. The Lagrange density describing 
the magnetic interactions of the nucleons is 

Ci.B = -^iV^Ko + «iT3) fT • B iV , (62) 

ZM 

where ko = ■^{np + Kn) and ki ~ ^{i^p — i^n) are the isoscalar and isovector 
nucleon magnetic moments in nuclear magnetons, with Kp = 2.79285 NM and 
Kn = —1.91304 NM. The magnetic field is conventionally defined B = V x A. 
The magnetic form factor of the deuteron at LO and NLO is 

4?(|qn ^ 2.o^tan-(M 



F«(|qn = (Z,-l)Fl?(|qn - 2^o(^. - 1) + L, ^ll^ 



(63) 
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Figure 11: The form factor A(q'^) (dominated by the square of the charge form factor) 
measured in electron-dcuteron scattering. The dashed, dotted and solid curves correspond 
to the LO, NLO and N^LO form factors. The solid circles correspond to data. 



Note that at LO the result is just the charge form factor muhiphed by the 
isoscalar anomalous magnetic moment. The first piece of the NLO result has 
a similar interpretation, but another piece also appears at this order. This 
second contribution represents the effect of a two-body current on Fm- The 
coefficient, L2, which appears there is the coefficient of a four-nucleon-one- 
magnetic-photon term in the Lagrange density EjE3 

C - -e^L2 ^e,,kiN^P^Ny (N^P.N) Bk + h.c. (64) 

The value of this coefficient can be found by demanding that the the deuteron 
magnetic moment is reproduced at NLO. Since the NLO result for hm in the 
Q-expansion isc3 

Mm = ^ I Kp + K„ + ^La 1 , (65) 

we reaiiire that if Li is evaluated at the scale [i = to^ then it must have the 
value BE3^L2(m^) = -0.149fm'*. This is significantly smaller than the naively 
estimated size of ^ 1 fm*. 

3.2 Deuteron Polarizabilities 

The electromagnetic polarizabilities of a hadron measure its, deformation in 
external electric and magnetic fields (for an overview see E^l). The electric 
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polarizability of the '^'^l/^Y/Or',. ctBOi has been investigated thoroughly with 
NN potential models L3£3C3'EII. To a very high precision one obtains a^o = 
0.6328 ± 0.0017 fm"^ with those techniques. This is not surprising since the 
electric polarizability is dominated by the long-range behavior of the deuteron 
wavefunction. Thus, if a model is tuned to reproduce this asymptotic part of 
the deuteron wavefunction then the predicted electric polarizability should be 
very close to nature. 

Similarly in EFT(7f) aso can be calculated to high precision. The con- 
tributions from the first four orders in perturbation theory give an electric 
polarizability of the deuteron which is (including a relativistic correction) 
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(66) 



where Dp = tC\''°' +2 ^cf^'^ -I- f ^C^ -1.51 fmMs a combination 

of P-wave coefficients. The same combination of P-wave scattering parame- 
ters contributes to the El amplitude in np d-y. Higher contributions are 
estimated to be < 1%. 

In fact, in general, insertions of higher-dimension operators will also pick 
up factors of Zd when computed to higher orders and therefore we provocatively 
write 
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= 0.6325 fm-" 



(67) 



The ellipses denote higher-order terms, which include the N^LO contribu- 
tion from '^Si — ^Di mixing proportional to 'tL ~ 6.5 x lO"*. We estimate that 
the uncertainty in the result shown in Eq. ( p^ ) is ±0.002 fm^, which includes 
the effect of this S-D mixing. The first contribution from interactions that are 
not constrained by NN phase shifts or gauge invariance occurs at N^LO, two 
orders higher than we have calculated, and hence contributes an uncertainty 
of ±0.0007 fm^ 



3.3 np d7 

The capture of low-energy neutrons by protons to form a deuteron is a reaction 
which is not only of fundamental interest in nuclear physics, but also represents 
a critical input to astrophysical simulations. At low energies the total cross- 
section is dominated by Ml-isovector capture from np continuum states to 
the deuteron. This is the energy region relevant for Big-Bang-Nucleosynthesis 
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(BBN), and ttese isovector amplitudes have been computed to very high or- 
ders in EFTE3'Eil. Meanwhile, the much smaller isoscalar amplitudes which 
contribute to this transition have been computed to NLO in order to make 
predictions that are to be compared with imminent experimental measure- 
mentsEJ. 

At low energies it is useful to make a multipole decomposition for the 
amplitude np dj 
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Only the lowest partial waves are shown, corresponding to electric dipole cap- 
ture of nucleons in a P-wave with amplitude Xei^ , isovector magnetic capture 
of nucleons in the ^Sq channel with amplitude Xmiv : isoscalar magnetic cap- 
ture of nucleons in the ^Si channel with amplitude Xmis^ and isoscalar elec- 
tric quadrupole capture of nucleons in the '^Si channel with amplitude Xe2s ■ 
Higher multipoles are suppressed by additional powers of k, the photon mo- 
mentum. Note that Un is the neutron Pauli spinor and Up is the proton Pauli 
spinor, while |p| is the magnitude of the momentum of each nucleon in the 
NN center-of-mass frame. The photon polarization vector is £(-y), and €(^4^ 
is the deuteron polarization vector. For convenience, we define dimensionless 
variables X, as follows: 
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where Xt^Li are the amplitudes other than Xeiv In terms of the amplitudes 
given in Eq. (^8|) and Eq. (^) the unpolarized cross-section is 
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(70) 



where a is the fine-structure constant, and |v| is the velocity of the incoming 
neutrons. Close to threshold Xmiv is the dominant contribution to this cross 
section, being larger than the other multipoles by several orders of magnitude. 
Therefore a measurement of a in this kinematic regime does not constrain 
Xeivj Xmisj or Xe2s- 
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The isovector Mly amplitude up to NLO is 



\X 



Mil 



Zd - f'O 



(72 + IPPf 



^ + |P|2 

ai I I 



0(i 



■ 7 



Ki 27r i 

ai 



7 



(71) 



where the first term can be obtained in ERT0. However, ERT cannot repro- 
duce the cross-section for this process, giving 10% less-jthan the experimental 
result for the thermal neutron capture cross-section Eil for incident neutrons 
with \v\ = 2200 m/s of Uth = (334.2 ± 0.5)mb. In thft_txaditional picture this 
10% discrepancy was explained by Riska and Brown E2l as being due to two- 
body currents. These mechanisms contribute significantly to the cross-section 
for thermal neutron capture at threshold. Below we will discuss work by Park, 
Min, and Rho, ps^u^h shows that this explanation emerges naturally in the 
EFT with pionsNy. 

In contrast, in EFT(7f) the pion is regarded as "short-range" dynamics 
and this two-body current appears as the isovector contact term L„p. The 
constant L„p is an RG-invariant combination of parameters. L„p is related to 
the couplings that appear in the Lagrangian via 
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where the Lagrange density describing four-nucleon-one-magnetic-photon in- 
teraction is 



C 



= e ''Li {N'^ P, N)\n'^ P" N) 5""^ B, + h.c. 



(73) 



The coefficient of this operator, which involves four nucleon fields and a mag- 
netic photon, is not related to NN scattering data by gauge invariance. Hence 
it must be fitted to data on the interactions of the magnetic photon with the 
two-nucleon system. The value for Lnp is determined by the requirement that 
we reproduce the cross-section a*"^. 

Once the value for the Mly amplitude is fixed at threshold its energy- 
dependence is determined by information from the NN system. But it ceases 
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to be the dominant contribution once incoming np energies of a few MeV are 
reached. The El amphtude becomes increasingly important as the energy of 
the emitted photon increases. This amphtude has been computed up to N^LO, 
and the result is: 



(7' + |Pp) 



(74) 



With this result in hand, and the value for Lnp fixed at threshold, the photo- 
dissociation cross-section for the inverse process, 7c? — > np can be determined. 
The results are shown in Fi g, p^ . 

The amplitudes in Eq. ([7lD and Eq. ( fz^ ) have been computed to one higher 
order by Rupalcll. This necessitates the inclusion of relativistic effects and a 
counterterm for the El multipole. The result is then a calculation of np — > dj 
which has an error of ~ 1%. Since the physics in these EFT calculations is very 
similar to the physics of potential-model calculations of this process it would 
seem that these "traditional" approaches can also calculate this reaction with 
comparable precision if they reproduce a*"^ accurately E3. Recent calculations 
confirm that this is indeed the case. 




Fiffiirp 1 9.: The photodissociation cross-sectifln for -yd —f—mfr-iThe solid line results from 
Eq. {UW evaluated with the amplitudes Eq. (plh and Eq. (MH. Lnp is determined by the 
cross-section for np 1 1 1^ 7, ct*''. The dashed lines denote the ~ 3% theoretical unf«iitainty 
estimated in RefE3. Rupak has further reduced this uncertainty to below 



We now turn out attention to polarized neutron capture. For a neutron of 
polarization rj^ incident upon an unpolarized proton target, the cross-sections 
for the production of right-handed and left-handed circularly-polarized photons 
are different. Defining the asymmetry (6') to be the ratio of the difference 
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and sum of these two cross-sections, we can express this observable in terms 
of the amphtudes X^^j-, 

Al^S^) = Vn[{P-y{Ml) + P^(E2)) co&e + P^{El) fiin^e], (75) 

where 



^/me\XMi^X*j 1 Re[XMiv^Pi 1 

\^Mlv\ \^Mlv\ 



Where we have neglected the small Mis, E'^'S and Ely amplitudes in the 
denominators. _ 
The isoscalar E2 amplitude has been computed to N^LO in Ref.cJ, 



= ^^^^-^{l + -{Za-l)-—{Za-lY\+—E\^>.m) 
4V2 10 V 8^ ' 32^ ' ) ^ ^ ' 

The deuteron quadrupole moment counterterm appears here, and once it is 
fixed there are no free parameters in this expression. Since Xmiv '^as com- 
puted to N^LO, this leads to a series for the ratio of isoscalar E2 to isovector 
Ml matrix elements which has the numerical values 



Xe2s 



= - [ 1.565 + 0.693 + 0.209 ] x 10"^ = -2.47 x 10"". (78) 



Note that this series appears to be converging well, with an expansion param- 
eter of about one third. 

The convergence of the isoscalar Ml amplitude is harder to assess. At 
NLO this amplitude is: 

4^l-V2^^L,(^-7)^ (79) 

Once again a counterterm that was fitted to a deuteron static property — in 
this case the magnetic moment — makes an appearance. While this is formally 
the leading contribution, the smallness of suggests that the contribution 
given in Eq. ( [79| ) might not dominate over higher order terms, and the A/I5 
amplitude might not be predicted well by EFT(7f-) at this order. Numerically, 

= -5.0 X 10-^ (80) 

with an uncertainty that we naively estimate to be of ~ 60%. 
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Using these amphtudes to compute the photon polarizations P-y, we find 

P^(Ml) = -7.1 X 10"'', P^{E2) = -3.5 x 10"'^, (81) 

giving a total of Py = —1.06 x 10"'^ in the forward dicection. This is approxi- 
mately 2/3 of the experimentally determined value olE3 p^'^v'^ = — (1.5±0.3) x 
10"'^. Given the large uncertainty in the calculation of the MI5 amplitude, 
and the uncertainty of the measurement, the two are consistent at this order. 
The valued of P^(Ml) = -2A x lO"'' agrees with the result obtained us- 
ing a Reid soft-core potentialEZl. However, it is significantly smaller than the 
ERT value of P-y(Ml) = —9.2 x 10"''. However, given the large uncertainty 
in the Mis amplitude of EFT(7f) these three values_are consistent. They are 
also consistent with the calculation of Park et al. Lj. This calculation will 
be discussed below, since it is a calculation in an EFT with explicit pions 
which uses the counting of WeinbergQ'El. The EElE(7f) value of P^{E2) quoted 
in Eq. ( ^l| ) agrees well with that found in Ref. E3. Note that the zero-range 
value of P^{E2): P^{E2) = -2.4 x 10"^, is reproduced at LO in EFT(7f), but 
it is modified by the presence of the deuteron quadrupole-moment countert- 
erm. Hence, the EFT(7f) value of P^{E2) ends up lying somewhere between 
this ERT result, and the result obtained in a Reid soft-core calculation EZI: 
P-y{E2) = -3.7 X lO-'^. 

3.4 Neutrino- Deuteron Interactions 

Weak interaction processes involving the deuteron are central to current re- 
search efforts in nuclear physics. In addition to the accelerator based programs 
that will elucidate the flavor structure of the nucleon, such as the SAMPLE ex- 
periments at Bates, the interactions between neutrinos and the deuteron form 
the core of our efforts to learn about the neutrino and look beyond the stan- 
dard model of electroweak interactions. Both in the production process, e.g. 
pp —> de'^Vf., and in detection at the Sudbury Neutrino Observatory (SNO), 
e.g. vd vnp, charged and neutral current weak interaction matrix elements 
between the deuteron and continuum states are required. 

Over the past few decades much effort has gone into calculating the pro- 
duction r nenha nism pp de'^v^, both from standard non-relativisticjauantum 
mechanicO'CHi, and from sophisticated potential model technique^. Over 
time, there have been a series of progressively more sophisticated calcula- 
tions of neutrino-deuteron scattering, lately with the goal of providing SNO 
withpiugh-nre/pision theoretical predictions of the cross-sections in various chan- 
nels OO'EjEaEi-ipIhe most sophisticated calculations, byi-YiHg, Haxton and 
Henley (YHH) BEa and by Kubodera and Nosawa (KN) blo, including the 
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effects_opii£:s£in-£xchange currents (MEC's), agree with one another to within 

Recently,, the.. process pp — > de^iy^ has been examined withJilFT by Kong 
and Ravndal E3Eill§3S and by Park, Kubodera, Min and RhoEJ. Elegant ex- 
pressions and numerical values for the weak capture cross-section that are con- 
sistent with previous estimates have been obtained. The cross-section depends 
somewhat on the value of one four-nucleon-one-weak-gauge-boson interaction, 
with an energy- independent coefficient e3 Li^a, defined by the axial currents 

Af^ = Li^A {N^PkN)^ (N^P^N) + h.c. , (82) 

and the vector currents 

V^^'-"' =2tL'^{N^P,N)H^, + ^,)iN^PkN) + h.c. . (83) 

Only Li^A contributes at NLO as it is the only two-body operator connecting 
the ^5*0 and "^Si channels. 
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Figure 13: Inelastic v(V)d cross-sections psia function of incident uijj) energy. The solid 
curves in the left graph are from referencellj, while the dashed curves are the EFT results 
at N ^L t p, ifit with Li ,4 = 5.6 fm^. The solid curves in the right graph are the results of 
YHHLZIO, and the dashed curves are the N^LO EFT resultsEl fit with Li_a = 0.94 fm^. 
In both graphs, the dashed curves all lie right on top of the solid curves. 

The neutrino-deuteron break-up cross-sections have been determined @0 
and it is found that Li a needs to be known at the 10%-level in order to perform 
a ~ 1% calculation of the cross-section. The YHH and KN numerical results 
can be recovered with different choices of Li.^, as shown in Fig. confirming 
that the difference between the potential-model calculations is short-distance 
in origin. Therefore, to predict the break-up cross-section with a precision 
of better than ~ 5%, one has two options. Firstly, compute the /3-decay of 
tritium, and use this to determine the EFT counterterm Li a, or equivalently 
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the MEC's in potential models |. Secondly, one can perform an experiment 
to measure one of the break-up cross sections to high accuracy, and thereby 
extract Li^a, or the MEC's. An important input into determining if neutrinos 
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Figure 14: The ratio of charged current to neutral current cross-sections in inelastic v — d 
and V — d scattering, as functions of incident energy. Shown are ratios of the EFT N-^LO 
results with Li.yi = —20 fm^ (dashed) and L^ a = 40 fm'^, ^gmonstrating the insensitivity 
of the ratio to the value of Li ^Ej- 



are changing flavor as they move through the interior of the Sun and to the 



Earth is the ratio of charged current to neutral current cross-sections. Fig. 14 
shows that this ratio, -Reft, unlike the cross-sections themselves, is relatively 
insensitive to the counterterm Li^a (or equivalently, the MEC's). 



a{E) = a{E) + b{E)Li^ 



(84) 



The cross-sections for two neutral current and one charged current neutrino- 
deuteron break-up channels are listed in Table |l|, parameterized by two energy 
dependent functions, a{E) and b(E), and the energy independent counterterm 
Li^A, with The magnitude of Li^a is estimated to be ~ 5 fm^ by NDA and 
is consistent with the values favored by potential models. The uncertainty 
introduced by omitting higher order effects is conservatively estimated to be 
~ 3% at threshold and ~ 2% at 20 MeV (using Li^a = 5.6 fm^ the cross- 
sections of ref. Lj are reproduced within ~ 2%). It can be seen from Table | 
that if Li A is constrained by a ~ 1% measurement of the v^d — > e~pp cross- 
section, then the uncertainty in the other channels would also be ~ 1%. 



s With certain assumptions alfoiit three-body forces, this method of fixing MEC's has 
been implemented for pp —> de^uS^- 
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E^^T, (MeV) 
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4 


0.0287 


0.00034 


0.0282 


0.00034 


0.146 


0.0016 


6 


0.188 


0.0024 


0.182 


0.0024 


0.574 


0.0067 


8 


0.514 


0.0070 


0.492 


0.0069 


1.34 


0.017 


10 


1.02 


0.015 


0.968 


0.014 


2.48 


0.032 


12 


1.72 


0.025 


1.61 


0.025 


4.01 


0.054 


14 


2.60 


0.039 


2.41 


0.038 


5.95 


0.082 


16 


3.69 


0.056 


3.37 


0.054 


8.31 


0.12 


18 


4.97 


0.077 


4.49 


0.074 


11.12 


0.16 


20 


6.47 


0.102 


5.76 


0.097 


14.39 


0.21 



Table 1: Neutrino-deuteron inelastic cross-sections from ref. l3, parameterized by ct(E) = 
a{E) + b{E)Li^A, where a{E) and b{E) have units of 10~*^cm^. ux denotes a standard 

model neutrino, u,,, Vf^ or v-r- 



3.5 ^-Deuteron Compton Scattering 

An area of nuclear physics that is stiU in its infancy is ^d Compton scattering. 
It is only recently that experimental measurements have been made of the dif- 
ferential cross-section for 70? scattering, but this has been limited to just three 
energies, = 49, 69 and 95 MeV. This process is important because it may 
provide the only rigorous way to determine the polarizabilities of the neutron. 
The amplitude for 7^ Compton scattering from an unpolarized deuteron target 
is 

A - -i^ [ Ti £ • e' k' • £k • e' ] £d ■ Ed , (85) 

where e and e' are the polarization vectors of the initial and final state pho- 
tons, and Sd and e'^ are the polarization vectors of the initial and final state 
deuterons. It is straighforward to compute the contributions to both form 
factors, Ti and T2, order-by-order in the Q-expansion. 

The differential cross-sections at E^ — 49 MeV and 69 MeV are plotted in 
Fig. |l^. At E^ = 49 MeV the theoretical expression is in very good agreement 
with the data at N'^LO. However, at E^ = 69 MeV there are significant differ- 
ences between data and the N'^LO cross-section, particularly at back angles. 
This is no surprise, as the theory should have failed at much lower energies. 
It is clear that the electric isoscalar polarizability of the nucleon, as computed 
at one-loop in chiral perturbation theory, makes a significant contribution to 
the 7c? scattering cross-section. There is good reason to believe that precise 
measurements of the cross-sections at low energies will allow for a good deter- 
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Figure 15: The differential cross-section for 7c( Compton scattering at = 49 MeV and 
69 MeV in EFT(7f). The dashed curve is the LO result. The dotted curve is the N^LO result 
with the nucleon polarizabilities, a and /3 set equal to zero. 'This solid curve represents the 
full N^LO result. The data is froniEZl. 



mination of the isoscalar nucleon polarizability and hence the neutron electric 
polarizability. 

4 The Three-Nucleon System 

The scalmg of the three-body force is the main issue to be understood in 
extending the results of the previous section to the three-nucleon system. De- 
spite recent progress this is not a fully understood topic and a complete power 
counting at all orders does not presently exist. A number of partial results are 
known however. The existence of the non-trivial fixed point, and the conse- 
quent non-perturbative resummation of the two-body interactions changes the 
scaling of the different three-body force terms from their NDA values. In most 
three-body channels this change suppresses the effect of the three-body forces 
and makes them appear in higher orders than naively supposed. However in 
the S"— wave, spin 1/2 channel (where the triton and '^He are) the opposite 
happens. The three-body force in this channel is enhanced and becomes a LO 
effect. The special physical property shared by no other channel is that, in this 
channel, Fermi statistics and/or centrifugal barrier repulsion does not block all 
three particles from being at the same point in space. One can then expect, 
on a intuitive basis, that this channel will have a larger sensitivity to short 
distance physics. We will see that this expectation is indeed fulfilled. 

When considering three-body systems, in addition to the one- and two- 
body Lagrange densities in Eqs. and ( p^ ) we have also three-body force 
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terms 

£3 = -DoiN^N)^ + ■■■, (86) 

where the dots represent terms involving more derivatives. 

Using NDA the coefficient D2r of a three-body force with 2r derivatives 
should be of order 

47r 

which would imply that their effect is very suppressed compared to the two- 
body interactions. Close to the non-trivial fixed point however, NDA does not 
have to hold and I?2r can be enhanced by powers of Aa or A/7. The natural 
size of D2r(A) is assumed to be set by the size required to appear in the low 
energy expansion at the same order as the cutoff dependence of the sum of two- 
body diagrams. Of course, for a particular value of the cutoff D2r(A) could be 
much smaller, but any small change in the cutoff would restore the "natural" 
value for D2r(A). This same criterion applied to the low energy expansion 
around the trivial fixed point is the usual NDA argument leading to Eq. (p7|). 
The proximity of the non-trivial fixed point leads to the necessity of a non- 
perturbative resummation of some of the interactions. The ultraviolet behavior 
of an infinite sum of diagrams can be very different from the behavior of each 
individual diagram and this affects the scaling of the D2r(A) coefficients. In 
all but one channel, i52r(A) will be even smaller than the NDA suggests. In 
the 5— wave, J — 1/2 channel (where the triton and •^He live) though, the 
opposite occurs and Dq is a leading order effect. 



4-1 The Spin ^ Channel 

In the spin 3/2 channel SEl the spins of all three nucleons are parallel and 
Fermi statistics forbids the three nucleons to occupy the same point in space. 
The three-body force term Dq cannot thus contribute in this channel and only 
three-body terms with more derivatives might contribute. Also, only the 
two-body interactions appear. Making the provisory assumption that the effect 
of the three body forces will be of higher order than our calculations we are 
left with a Lagrangean involving only one- and two-body terms: 



C1+C2 = 



lOn H 

' 2M 



(iV^p'7V)t(7V^pa(y _ y)2^) + h.C.j+--- 
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It is very convenient in three-body calculations to rewrite the theory described 
by Eq. ( p8| ) by introducing a dummy field with the quantum numbers of the 
deuteron (referred to hereafter as the "deuteron" ) 
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yd [d*^(7V^P''iV) +h.c.] 



(89) 



By performing the gaussian integration over the deuteron field and treating its 
kinetic term as a perturbation we recover Eq. ( p8[ ) with time derivatives instead 

of space derivatives in the term proportional to ''C^ (and some higher order 
terms). This time derivative can be traded by space derivatives by use of a 
field redefinition ^ 

(90) 



N + 



The parameters in Eq. (|8|) and Eq. (|8^) are related by 



A, 



(91) 



With this choice of parameters Eqs. ( p8| , p9| ) describe equivalent theories, i.e., 
they both give rise to the same S-matrix. Notice that after the introduction of 
the deuteron field there is a large reparametrization invariance in the theory. 
We could for instance, change the coefficient yd and add an explicit two-nucleon 
interaction is such a way that after eliminating the deuteron field the same 
Lagrangean was recovered. Thus the Lagrangean in Eq. (|8^) is not the most 
general one involving nucleons and deuteron fields. The important point is 
that after the elimination of the deuteron field we recover the most general 
Lagrangean involving nucleons only. From Eq. ( |9l| ) and Eq. (20) we arrive 
at the scaling ~ {AttK)/M and ~ (47rA/i)/Af. The scaling above shows 
that the kinetic term of the field is indeed sub-leading compared to the term 
proportional to A^. 

At tree level the d* propagator is given by the constant i/Ad- The loop 



corrections shown in Fig. 16 contribute to the same order as the tree level 



term. Summing up these graphs gives 



47ri 



47rAd 



SI 



(92) 
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Figure 16: The leading order deuteron propagator containing the sum of bubble graphs. 



The leading order (T field propagator is, up to a constant, the two-nucleon 
scattering amplitude in the ^5*1 channel. 

We consider now the scattering of a neutron on a deuteron with total spin 
S' = 3/2 (the more experimentally accessible proton-deuteron scattering would 
include Coulomb effects that would mask the strong interactions we are inter- 
ested in here). We use the field as the interpolating field for the deuteron. 
The graphs that contribute to this process at leading order are depicted on 
the top of the Fig. A simple power counting shows that they are all of the 
same order. The tree level diagram, for instance, contains two dNN vertices 
(~ A/M) and one nucleon propagator (~ M/Q'^), for a total contribution of 
order ~ A/Q^. The one loop propagator contains two additional dNN ver- 
tices, two extra nucleon propagators, one field propagator (~ 1/My'^Q) and 
one loop (~ Q^/M), thus the one loop graph is also of order ~ A/Q^. All 
the graphs shown in the first line of Fig. |l^ are of the same order and must 
be resummed at LO. This is very similar to the two-nucleon scattering case 

where all graphs involving only the leading interaction proportional to ^Cq '^^^ 
are leading. The difference resides on the fact that, while all the graphs for the 
LO two-nucleon scattering are easily summed in a geometric series, the same 
cannot be done in the three-body sector. The best we can do is to write an 
integral equation (shown in the second line of Fig. |l^ whose solution, obtained 
numerically, equals the desired infinite sum^. 

With the kinematic shown in Fig. the S*— wave amplitude satisfies 

toik,p) = -'^Ko{k,p)~ - [ dqq''to{k,q) / Koiq,p) 

(93) 



''Clearly the validity of this integral equation goes beyond all orders in perturbation 
theory even though the derivation shown here is based on diagrams. 
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Figure 17: The first line shows the infinite sum of graphs giving the LO neutron-deuteron 
scattering in the spin 3/2 channels. The second line is the graphical representation of the 
Fadeev equation that sums them up. 



where E = 3fc^/(4Af) — is theJx)tal energy. The projection on other 

partial waves lead to similar equationsEj. Notice that when p = k, all external 
legs are on-shell. This equation is just .the Faddeev equation for the case of 
contact forces derived first in reference by different methods. Eq. ( |93| ) is a 
one dimensional integral equation (the variable k appears only as a parame- 
ter) and its numerical solution is computationally trivial. A little care must be 
taken with the singularities of its kernel however. Although the pole in the real 
axis due to the deuteron propagator is regulated by the fact that the energy 
should be thought as having a small positive imaginary part and the logarith- 
mic singularity occurring above threshold is integrable, both cause numerical 
instabilities. An efficient|-s5:ayj|af dealing numerically with those instabilities is 
the Hetherington-SchickM'L3€3 method. 

To obtain the neutron-deuteron scattering amplitude we have to multi- 
ply the on-shell amplitude to{k,p = k) by the wave function renormalization 
constant, 



with 




(94) 




(95) 



where 



1 





(96) 



apo iA(p) po=-^J,i3^ 



i3=0 



87r7 
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Unlike to{k,p), the scattering amplitude To(fc) depends on the parameters 
and Ad only through the observable 7 — + ^J.■ 

At NLO, we have contributions coming only from insertions of the kinetic 
energy in the deuteron propagator. For book keeping purposes it is convenient 
to split the constant into A^,-! + A(j,o with 

Adfi is subleading compared to the leading order A^,-! operator and is treated 
in perturbation theory. With this choice the position of the deuteron pole is 
not changed by the NLO corrections to the deuteron propagator and Arf__i 
does not have to be refitted. 




Figure 18: The NLO contributions to nd scattering in the quartet channel. 

The NLO corrections to neutron-deuteron scattering are shown in ^'E Fj^ 
They are simply the range corrections in first order perturbation theoryOL3o. 
The graphs in Fig. ^ do not mix different partial waves and give the contri- 
bution 



00 

itl(k) = J ^q^MK^)? 



(98) 



. 3k2 _ 7d _ q' ^^ 



where the correction to the deuteron propagator isXd(pojP) = Po'^'iM' 
The amplitude in Eq. ( |9^ ) has to be solved for numerically as to{k,p) is known 
only numerically. The diagrams in Fig. ^ are ultraviolet finite. The degree of 
convergence however is not what one might naively guess, as can be seen by 
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the following argument. Each one of the diagrams in Fig. behaves for large 
values of the external momenta p ^ 7 as Their sum however has a faster 

decay. To see this consider Eq. (|9^) for large values of p ^ fc ~ 7. In this case 
the integral is dominated by tli£.^mgion of integration p ^ q where the kernel 



simplifies and we are left with @E 

00 

toik,p)^~^ fdqto{k,q)-ln(^±^^^) . (99) 

VSttJ p \p-pq + q^J 



Using the ansatz tti{k,p) ^ p^~^ we find an equation that determines the power 
describing the decay of to{k,p) for large p 

^ ' ^ - 0. (100) 



V3s cos(7rs/2) 

The negative solution (the positive ones lead to divergent integrals, cannot 
describe solutions of the integral equation and are discarded) is s ~ —2.166. 
The loops in Fig. |l| converge in the ultraviolet as A^^^" A ^-^^ as opposed 
to the naive behavior A~^, valid at each order in a perturbation expansion 
in the strength of the kernel. This fact has consequences for the estimate of 
the three-body forces since it changes their running. As we will see below, in 
the triton-^He channel the opposite happens and the three-body force in that 
channel is enhanced compared to NDA estimates. 

To complete the NLO calculation the wave function renormalization con- 
stant Z is found from 



1 1 1 Zi 1 . d 



— 1 



Z Zo + Zi Zo Z^ Zo dp, 







iJd(po,p) -,2 ^ (101) 



PO = T#,P=0 



and thus Zi ~ Zq. The NLO amplitude in the quartet channel is therefore 
given by 

T{k) = Z t{k,k) ~{Za + Zi)(to(k,k)+ti{k,k)) 

~ To{k) + Zohik,k) + ZMk,k)=To{k)+Ti{k) . (102) 

The phase shift for each partial wave is extracted by expanding both sides of 
the relation 

Tik)^Mk)+Uk)-'ij^J^, (103) 



Q with 5 = (5^"^ -|- S^^^f + . . . and keeping only linear terms, 

^(°) = ilnfl + ^To(fc)) , ^(^) = lA(foW.+ ^iW) 
2i V Stt V i Stt (1 + ^To(fc)) 
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At N^LO there are two kinds of contributions. One comes from the same 
insertions of deuteron kinetic energy and A^*^^ appearing at NLO, but now 
inserted twice. The other comes from inserting the LO operator describing the 



tensor force appearing in Eq. ( |105| ). The leading tensor force is generated by 
the term 

C'^^ = y'U'^N{{V - V)*(V - - \5''{V - vf)P'N (105) 



in the Lagrangean Eq. ( p9D after the gaussian integration over the deuteron 
fields is performed. This operators mixes spin and angular momentum and 
produces a splitting between the different three-body amplitudes amplitudes 
with the same spin S and angular momentum L but different values of the to- 
tal angular momentum J. This splitting is very important in determining spin 
observables like the vector analyzing power Ay. However the N^LO calculation 
of the amplitudes is only a leading calculation of the splittings between differ- 
ent values of J and cannot be very precise. Further, the tensor force appears 
only linearly and after averaging over the different values of J, its contribution 
vanishes. Here we will restrict ourselves to the averaged phase shifts and, con- 
sequently, drop the tensor force contribution at N^LO. The N^LO calculation 
(including some higher order contributions that can be ignored) is simply per- 
formed by replacing the LO deuteron propagator in Eq. (^2|) by the propagator 
which includes all effective range effects 

iA*^(p) = - , . (106) 

^^Vd _^ + lp,{Mpo + + Mpo 

in the integral equation in Eq. (]9^). 

The calculation for the spin 5 = 1/2 channels proceeds in an analogous 
way to the one in the spin 3/2 channels discussed above. The main difference 
is that two-body interactions in the spin singlet channel also contribute. We 
will introduce another dummy field T"", now with the quantum numbers of two 
nucleons in the ^Sq channel 



92 

(-ido + ) - At 



r + yt [t'''f{N^P''N) + h.c] (107) 



with the constants yt and Ai satisfying relations analogous to Eq. (|9l|) . Again, 
by integrating over the "dibaryon" t"^ we recover the ^Sq two-body interactions. 
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Figure 19: The coupled equations in the spin 1 /2 channels. 



The spin zero dibaryon field contributes in intermediate states of neutron- 
deuteron amplitudes. The analogue of Eq. (93) is then a system of two cou- 
pled integral equations for the d + N d + N {itii{k,p, e)) and d + N 
t + N{itt{k,p,e)) amplitudes. 



tdip, k) 



ttip, k) 



Koip,k) 



dq 

dq 

q^ — fc^ 



Ko{p,q)[td{q,k) + itt{q,k)] (108) 



Ko{p, q)[itd{q,k)+tt{q,k)] 



The results of LO, NLO and N^LO calculations 



of the phase shifts 



for the lowest partial waves (averaged over the different values of J) are shown 
in Fig..2p. For comparison we include the results of some old phase shift anal- 
ysesE§E3(that unfortunately do not contain er.rni:. bars) and a recent extensive 
numerical calculation using modern potentialsll2£lll22l. In most cases the N^LO 
EFT computation seems to be convergent and accurate at the precision one 
would expect {{'l/rnT^Y ^ 3%). This theoretical error is probably an overes- 
timate because the N^LO corrections come from one insertion of the shape 
parameter, that is unnaturally small, three insertions of the range correction, 
that were already included in the N^LO calculation as discussed above aijd ijine 
P— wave insertion. There is a very precise experimental determination 



the quartet S*— wave scattering length 



of 

6.35±0.02fm. The N^LO result, 



6.33 ± 0.05 fm, compares well with this value and with calculations 



103 



EFT 
"3/2 

using second generation potential models 

Since the scattering lengths ad, at, the deuteron binding momentum I/7 
and effective range rot of the two body system in the S-waves are all input 
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Figure 20: Real part of the first partial waves in the quartet (left) and doublet (right) 
channels of nd scattering versus the center-of-mass energy. The dashed line is the LO, the 
dot-dashed the NLO and the solid line the N^LO result. 
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to the EFT calculation, we expect that any potential model which reproduces 
these numbers will agree-Jidth our results within the estimated error of an 
N^LO calculation, ~ 4%E23. Higher order calculations in EFT can be easily 
performed and involve only well known two-body interaction parameters. They 
too should agree with the latest generation of potential models. Only at the 
order where three-body forces first contribute will differences between EFT 
calculations and realistic potential models arise. In some observables, like spin 
observables in neutron-deuteron scattering, the contribution from three body- 
forces appears at lower orders and it is here where differences between EFT 
and potential models will first become apparent. 



4-2 The Triton Channel 

We know turn to the spin 1/2, 5*— wave channel. The essential physics point 
that distinguishes this channel from all the others is that, in the triton-'^He 
channel, all three nucleons can occupy the same point in space while in all 
the others the centrifugal barrier and/or Fermi statistics creates an effective 
repulsion between the three particles at short distances. This suggests that the 
ultraviolet behavior of this channel should be similar to the case of three spin- 
zero bosons in a S— wave that also does not have any short distance repulsion. 
It turns out that the similarities are more than qualitative and the spin-zero 
system has all the novel features we want to illustrate. We will discuss first 
the spin-zero boson case to avoid unessential complications due to the coupled 
equations in the triton-'^He channel. After the conceptual issues are discussed 
in the bosonic case we will turn to the triton channel. 

The integral equation analogous to Eq. (^3|) for the three bosons case is 

y^M ^ fp"^ +pk + P - ME 
2pk Vp^ -pk + k^ - me 



0(J 

2 f 1 

— / dqq^ to{k,q) Ko{q,p) 



The difference from the quartet case (Eq. (|93|)) is that the kernel in Eq. (109) 



has the opposite sign and twice the strength. This is enough to make-ilie 



Eq. (109) ill-posed and it can be shown that its solution is not unique tZIEa. 



Some physical insight into the situation can be obta|iii^d by going to coordinate 
space. For that, let us perform the transformation 



104 



p sm{y/3/ App) _ /i 



Xir,p)= dptoik,p)^ . ^r. ^-.v^-^"^. (110) 

4 
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The wave function x(r, p) satisfies 



Qj.2 Qp2 



+ ME]x{r,p)^0, (111) 



with the boundary conditions 



X(r,0) = , |-;^(0,p)-7x(0,p) + ^^x(^P,^p) = .(112) 



The last boundary condition is equivalent to Eq. ( [1091 ). Eqs. (111,112) could 



also be obtained by considering the three particle problem in the wave-mechanics 
formalism and using Jacobi coordinates. The variable r corresponds to the 
distance between two of the particles and p to the distance between the third 
particle and the center of mass of the other two. Using polar coordinates 
R = i/r^+p^, tana — r/p, we can see that this co mplic ated boundary con- 



dition is simpler for i? <C a = I/7. In this region Eq. ( |112| ) becomes 



^x(i?,a = 0) + ^x(^,7r/3) = 0. (113) 

In the region i? <C a we have then a separable problem with the solution 

X(i?,a) = ^i^,(i?) sin(s,(7r/2-a)), (114) 

i 

where the Si are the solutions of 

^_ 8 sin(W6) (115) 
y/3s cos(si7r/2) 



Eq. (118) has two imaginary solutions ±zso, with sq ~ 1.006. The equation 
for the radial part Fq{R) is that of a two dimensional Schroedinger equation 
with an attractive potential 

As is well known the quantum mechanics problem with an attractive \/R^ 
potential (or more singular potentials) is not well defined and its hamiltonian 
is not bounded from below. There are arbitrarily deep bound states so the 
ground state describes the "fall to the center" . Of course this conclusion is 
not valid as soon as the bound states start probing distances of the order of 
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the cutoff. It indicates though that the spectrum wiil be dependent on short 
distance physics. 

Another way of analyzing this system is to consider the equation describing 
the asymptotic behavior of the half off-shell amplitude, as we did in the quartet 
case in Eq. (p9|). In the case of three bosons we have 



pto{k,p) ^ [ dqqto{k,q) In ( ^9 ^^'^T^2 ) ' (l^'^) 



where we introduce a parameter A taking values A = 1 in the bosonic case and 
A = —1/2 in the quartet case. Besides the symmetry under scale transforma- 
tions p tQ{k,p) ^ a p t(j{k,ap) Eq. ( |117| ) has a symmetry under inversions 
p to{k,p) ^ to{k,^) that will be crucial later. Using the same ansatz 
to{k,p) ^ as before we arrive at an equation determining the asymptotic 
behavior of to{k,p) 

l-4^^^^=0. (118) 
V3s cos(s7r/2) ^ ' 

For A > Ac = — 0.4135 Eq. ( |118D has two imaginary solutions s — ±iso, 
with sq ~ 1.0062 in the bosonic case A = 1. This means that for p ^ fc,7, 
to{k,p) is a linear combination of p^''"'^ and p^^^o^i or 

toik,p) = - cos(so log(p) + 6). (119) 
P 



The amplitude A is determined by matching the asymptotic solution in Eq. (|119| ) 
to the full solution in the infrared region. The phase however remains undeter- 
mined. This is another way of seeing that there is a one parameter family of 



solutions to the equation Eq. (109). Of course this undctcrminacy is an artifact 
of extrapolating the low energy effective theory to the ultraviolet region and 
is eliminated by the introduction of an ultraviolet regulator. 

A perturbation theory on the strength of the kernel (A) corresponds to the 



diagrams on the top line of Fig. \^ The behavior shown in Eq. (11£) is missed 
in all orders of this expansion as can be seen by expanding Eq. ( |119D in powers 
of the kernel strength A ( the dependence on A is implicit in sq) 

j\ 

- cos(so log(p) + (5)=0 + 0A + 0A2 + -- - (120) 
P 

This explains why a perturbative analysis of the ultraviolet behavior of to{k,p) 



misses the oscillatory behavior in Eq. (119) 
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Figure 21: Half ofT-shell (K-matrix) amplitudes in the bosonic case for different values of the 
momentum cutoff A in units of the two-body scattering length a2 ■ 



In Fig. ^ we can see numerical solutions of Eq. ( |l09[ ) (or more preci sely, the 
K-matrix corresponding to the amplitude to{k,p) satisfying Eq. ( |l09|) with a 
principal value instead of an ie prescription) for k = obtained with different 
values of the momentum cutoff, A. The presence of A breaks the inversion 



symmetry and selects a phase in Eq. (119), and further, different values of A 



result in very different values of to(k = 0, fc = 0). This A dependence of a 
physical quantity indicates something is missing in our calculation up to now. 
As it appears in the three-body amplitudes it should be absorbed by the cutoff 
dependence of a three-body force. We see now that we were unjustified in 
dropping the three-body force term from the integral equation in the case of 
bosons. To reinstate the three-body force shown in Eq. ( |8^ ) we include the 
term 

C, 3 = dUN^N (121) 

in the Lagrangean and choose H = Do ^./^i ■ The integration over the field d 
generates the three-body force in Eq. ( |8q). Graphs containing the three body 
force are now included as shown in Fig. |22| 

The new graphs with a three-body interaction change the kernel to 



qp 



In 



p'^ +pq + q'^ - ME\ H 



p^ -pq + q^ - ME ) A^ 



(122) 



It is not currently known how to dimensionally regulate the three-body integral 
equation, and from now on we will use a sharp momentum cutoff instead (some 
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Figure 22: Equation for the three-body ampUtude including the three-body force. 



progress in this direction was made in the case of a similar equation inE££l). To 
be more precise, we wiU be using two distinct cutoffs: one regulating the two- 
body loops in the deuteron propagator that is taken to infinity and another 
regulating the three-body loops. That is why the kernel of the three-body 
equation contains the renormalized form of the deuteron propagator. This 
procedure is justified since keeping an unrenormalized form for the deuteron 
propagator amounts to keeping some sub-leading terms in the kernel that are 
suppressed by powers of g/A. These terms are relevant only at the cutoff scale 
g ~ A and can change only non-universal features, like the value of the bare 
coupling constants, but not the final physical results. 



For i?(A) ~ 1 the effect of the three-body force in Eq. (122) is important 
only for g ^ A thus our analysis of the asymptotic behavior of the amplitude 
remains valid in the intermediate region 7, p A. A simple argument 
now shows that H{A) can be chosen in such a way as to maintain the low- 
momentum part of the amplitude to{k,p <C A) and, in particular, the on-shell 
point p = k, cutoff independent (up to terms o f ord er k/A). For that we define 
a scale 7 ^ ^ ^ A and rewrite equation Eq. ( |l09| ) as 



4(fc,p) = K{k,p) + - r dqqh^{k,q)Ko{q,p) 

Jo 

+ dqqh^{k,q){Ko{q,p) + ^) + 0{p/A,^/ A). {123) 

The superscript A reminds us that (fc,p) is the solution of the integral equa- 
tion with a particular cutoff A. Taking another cutoff A' we can write the 



equation for <g {k,p) as 



fj- 



t^\k,p) = K{k,p) + - dqqX'{k,q)Koiq,p) 
^ I' dqqX(k,q){Ko{q,p) + ^) 



TT 
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A' 



[ H{A') 



1 



dqq to {k,q) , . 



If H{A) can be chosen in such a way as to make the terms in the last two 
Hnes of Eq. (124) vanish, the equations satisfied by tQ{k,q) and ig {k,q) are 



the same and, since the equation with a finite cutoff has an unique solution, 
the amplitudes are the same. Using the asymptotic form of the amplitudes 

4(/c, 7 < p < A) = 4 (/c, 7 < p < A') = - cos(so log(p)/A), (125) 

where A is an arbitrary scale fixing the asymptotic phase we find an approxi- 
mate form for H{A) 

sin (so log [A/AI — arctan fl/sol) 

H{A) = ; J i (126) 

sin (so log [A/AJ + arctan [I/sq]) 



which assures A independence. The precise form of H{A) in Eq. (126) might 
dependent on the particular shape of the cutoff, but its qualitative behavior, 
in particular the existence of a limit cycle in RG flow, are general. 

We can check numerically the argument leading to Eq. ( |126| ). For this we 
can fix (arbitrarily) the value of the three-body scattering length and determine 
numerically the value of H{A) needed to reproduce this same scattering length 
for any given value of A. This can always be done since we have one free 
parameter {H{A)) to fit another (the three-body scattering length). Using 
these determined values of H{A) we can check that the amplitudes at finite k 
are indeed cutoff independent (up to terms of order fc/A,Aa). The result of 



this procedure is shown as the dots in Fig. g^. The agreement with Eq. (126) 
lends support to the arguments leading to it. 

The appearance of a new parameter (that can be taken to be the value 
H{A) at some particular value of A or the scale A) that is not determined by 
two-body physics is the main result of this analysis. It implies that different 
models having the same two-body scattering length can have quite different 
three-body physics. This is shown in a striking manner by plotting the pre- 
diction of a number of different models with very similar low energy {p < ttIt^) 
two-body physics for two three-body observables like, for instance, the three- 
body bound state energy (B^) and the sC|attering length of a particle colliding 



with the bound state of two others (as) It is known from the numerical 
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Figure 23: R iiTin 
Eq. (126D 



ng of the three-body force around a limit cycle. The solid line shows 
and the dots the numerical procedure described in the text. 



studies of potential models that those predictions can differ wildly among dif- 
ferent models, but not to scatter all over the x B3 plane cither. They tend 
to lie along a reasonably well defined line (Phillips line, see Fig. For the 
point of view discussed here each point on the Phillips line correspond to a 
different value of the three-body parameter A. 

The existence of the Phillips line is usually understood as the result of the 

' off-shell extension 



, LO' 



This is equivalent 



dependence of the three-body observables on the particuJ 
of the two-body data described by the different potentialsE 
to the explanation discussed above since field redefinitions change both off- 
shell behavior of amplitudes and strengths of three-body forces. Looking at 
those "off-shell ambiguities" as values of short distance counterterms has the 
advantage of being a systematic way of mapping them out and estimating their 
sizes. In particular it becomes clear that the "off-shell ambiguity" related to 
the three-body bound state cannot be fixed by two-body processes involving 
external currents since those are related to different counterterms. 

We now have the undej:atanding necessary to tackle the triton-'^He chan- 



nel, described by Eq. (108)1123. There is an interesting limit where Eqs. (108) 
decouple. The two spin and isospin states of the nucleons can be put in a 4 
representation of spin-flavor SU{A). This symmetry is broken in leading or- 
der by the difference between the scattering lengths in the '^Si (7) and ^Sq 
(a) channels but is an approximate symmetry of the system when the typical 
momenta is Q ^ 7,1/a since in this case we can take 7 ~ 1/a ~ 0. This 
approximate symmetry is observed in light nuclei and has interesting conse- 
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Figure 24: Results for the doublet S— wave noutron-deuteron scattering length for different 
models having very similar two-body scattering lengths and effective ranges, denoted by the 
dots. The solid line is the EFT prediction over a range of values for A. The cross indicates 

the experimental value. 



quences 
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In the 5/7(4) limit the system of equations Eq. (|l08| ) break up 
into two equations: one for the sum t+(k,p) = td{k,p) + tt{k,p) and the other 
for the difference t-{k,p) — td{k,p) — tt{k,p). The three-body force in Eq. (86) 
contributes only to the bosonic-like equation for t- {k, p). It is actually the only 
three-body force with no derivatives one can write down, as a simple argument 
shows. Fermi statistics mandates that the three nucleons should be in a com- 
pletely antisymmetric representation of 5/7(4). Consequently, the three-body 
force transforms as 4(8)4 = 1©15, but only the singlet is separately invariant 
under the spin and isospin subgroups. The three-body force in Eq. (|8^) is 
generated by the inclusion of the term 



2MH{K) 



A2 



ylN\d'a')\d^(j^)N + y^^N\eT'')\t''T^)N 



+ -ytys[N\d'a')\eT'')N + h.c: 



(127) 



Upon integrating the field and out we recover Eq. (86). Fitting A to 



the value of the experimental doublet scattering length ol^i^ we arrive at the 



alid curve in Fig. g5|. Also shown in Fig. are the phase shift analysis of 
3€3 and the experimental value of the doublet scattering length. Using those 
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Figure 25: Doublet S— wave neutron-deuteron phase shifts. 

parameters the prediction for the triton binding energy is of i?3 = 8.03 MeV, 
to be compared with the experimental value B^^^ = 8.48 MeV. The dashed 
line was obtained by using the experimental value of the two-body scattering 
length, instead of the deuteron binding energy, to fix ^Cq^^\ while keeping 
the three-body force constant. Both procedures differ only by higher orders 
terms so their difference provides an estimate of the higher order corrections. 
Comparison with experiment suggests that the low energy expansion in the 
triton is under control since the errors are within the expected value for a 
LO calculation. A more solid statement can be made only when higher order 
calculations become available. A natural question to ask then is whether the 
three-body force needed to give the triton its experimental binding energy is 
"natural". For cutoffs of the order A ~ m^r , 171^^) ~ 1 suggesting that 

the relative shallowness of the triton does not require any further fine tuning 
besides that in the two-body sector. 

Compared to the two-body sector the study of the three-nucleon prob- 
lem with EFT methods is still in its infancy. A number of conceptual and 
technical problems have to be understood before high order calculations can 
be performed. Among them are the complete power counting to all orders 
and the scalinfriif the sub-leading three-body forces (some recent progress was 
report in Ref.EiJ). 



57 



5 NN Interactions at Higher Energies 

In order to make progress toward describing multi-nucleon systems, such as 
nuclei, an EFT must be constructed that wiU describe the NN system at 
momenta |p|^ TOtt- It is obvious that one must include the tt- field explicitly 
in this kinematic region. However, a power-counting must be developed that 
is both systematic and convergent. In what follows we will describe the two 
different power-counting schemes that currently exist in this kinematic regime. 
Weinberg's power-countin^'Ei is based on a momentum and chiral expansion 
of the NN potential, in which the momentum independent, local four-nucleon 
interaction and one-pion-exchange (OPE) both enter at leading order. Higher 
derivative operators and the light-quark mass matrix occur at higher orders 
and are inserted in perturbation theory. Calculations that have been performed 
to date, indicate that the perturbative expansion appears to be converging 
in the two-nucleon sector, as we shall discuss subsequently, however, there 
are formal issues that remain to be resolved. The jury is still out in the 
three-nucleon sector, but encouraging results are being obtained. In contrast, 
KSW power-counting provides an explicit framework in which to systematically 
expand S-matrix elements as opposed to an expansion of the (unobservable) 
potential. Calculations are renormalization scale independent order-by-order 
in the expansion parameter. However, a-calculation of S-wave NN scattering 
at N^LO by Fleming, Mehen and Stewarti£l indicates that, in fact, treating the 
pions perturbatively is not converging for momenta |p|^ They were able 

to identify the one large contribution that occurs at N^LO, and found that it 
persists in the iriq — > chiral limit. Hopefully this problem can be understood 
and solved in the near future. 

This is essentially the current status of the field. In what follows both 
power-countings will be discussed, both their strengths and currently under- 
stood weaknesses. Emphasis will be placed on the Weinberg counting due to 
its present numerical successes and lack of numerical failures, despite some 
unresolved formal issues. 



6 Perturbative 7r-Exchange and its Problems 

Let us recall how pions couple to nucleons. The Lagrange density describing 
the interactions between pions and nucleons is 

£ = ^TrLi^ED^S^ + ^-^N^a ■ (^D^^ - ^^D^)iV + ■■■ , (128) 
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where the pion fields are contained in a special unitary matrix, 



E = = exp 



V2iU 



n 



TT' 



■VV2 



) 



(129) 



with / = 93 MeV and gA ~ 1-26 is the axial coupling constant. The ellipses 
represent operators involving more insertions of the light quark mass matrix, 
meson fields, and spatial derivatives. The leading-order contribution to the 
NN potential due to the exchange of a single "potential pion" fl is: 



where a'' denotes the spin operator acting on the fcth nucleon. It is clear 
that this interaction does not have any dependence upon the renormalization 
scale ^. Since the Co operator scales as l/{Mfi) it would appear that as 
the system is evolved into the infrared, the short-range contact interaction will 
dominate over one-pion exchange. From the standpoint of conventional nuclear 
physics this is a shocking claim. The tensor force in the ■^Si — ^Di channel is 
known to be strong, and this casts doubt on the applicability of a perturbative 
treatment of the pion exchanges. The question that must be addressed is 
what is the kinematic regime where this hierarchy is actually present. The 
observation that there is a kinematic regime where pions can be .treated in 
perturbation theory was first made by Kaplan, Savage and WiseE2l'E2l (KSW), 
and the associated power-counting worked out in these papers is commonly 
known as KSW-counting. KSW-counting has the very attractive feature that 
analytic expressions for any scattering amplitude can be found. 

6.1 NN Scattering with Perturbative Pions 

The NN scattering amplitude has now been computed up to N in KSW 

counting. The LO contribution to the amplitude is exactly the same as in the 
pionless theory, since, as mentioned above, the pion scales as while Cq, 
and the whole bubble chain contributing to y^*^^^^ scale as Q^^. Thus, once 
more: 



'The variation of the pion propagator due to the energy transfer between the two nuclcons 
is negligible compared to the variation due to the three-momentum transfer. 




(130) 




(131) 
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Figure 26: Feynman diagrams contributing to A'^A'^ scattering in the ^5o channel. 



Meanwhile, the Feynman diagrams contributing to the NLO amphtude 



are shown in Fig. M and their contribution is the sum Aq ~ X'q^ + -4q 



•^0 ^ •^o 



{IV) 



J\q . The local operators at this order involve either two 



spatial derivatives, with coefficient " , or one insertion of the light quark 
mass matrix, with coefficient D2 ° ■ Thus, the amplitude is 



A, 
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(132) 



The resulting phase shifts at LO, NLO0 and N^LoEil are shown in Fig. ||. 
It should be noted that these results are for a particular fitting procedure 
and appear to be converging very well. Detailed discussions of various fitting 
procedures, and whether the pion can truly be treated in perturbation theory 
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in the ^Sq channel, can be found in Refs. E§EJOHOE3IH; 
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Figure 27: Fit to the ^Sq phase shift S. The sohd hne is the Nijmegen fitQ to the data. The 
long dashed, short dashed, and dotted Hnes are the LO, NLO, and N^LO results respectively. 



A similar calculation has been performed in the ^Si — ^Di couplerLcl 
but we will not present the expressions for the amplitudes herelEI'ti^ 
resultant function yields the '^5i phase shifts shown in Fig. Ea. From Fig. 28 
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Figure 28: The ^Si - ^Si phase shift, <5o and ^Si - ^Lh mixing parameter, si for A'^A'^ 
scattering. The solid line is the Nijmegen multi-energy fitQ, the long dashed line is the LO 
effective field theory result, the short dashed line is the NLO result, and the dotted line is the 
N-^LO result. The dash-aiotted line shows the result of including the parameter (^5 (defined 
in Ref.lli3) which is higher order in the power counting 



it is clear that the N^LO phase shift is deviating significantly from the data 
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at momenta |p|^ m7r/2. Given that one expects the perturbative series to 
converge to the data at successively higher orders, it would appear that the 
KSW-expansion is in trouble. 

This conclusion is reinforced by looking at the ^Si—^Di mixing pararneter, 
ei, defined above. This observable has also been computed up to N^LCo. The 
results at LO, NLO and N^LO are shown in Fig. GH. 

The curves in Fig. ^ suggest that KSW power-counting is failing in the 
•^^i — "^Di coupled channels. The origin of this failure in the '^S'l-channel at 
N^LO has been tracked down to a term that survives in the chiral limit 



2 M fMgy ^ 



77 ipr 



(133) 



where '^^^ 



is the LO ^Si - ^Si amplitude. The N^LO amplitude in 
Eq. (133) is large because of the coefficient of 6 which is to be compared with 
the expansion parameter of ~ 1/3. The fact that this correction survives in the 
chiral limit indicates that it comes from the short distance part of potential 
pion exchange. The problem one encounters in this coupled channel is really 
due to the large strength of the tensor force. Of course, in the ^Si — "^Si ampli- 
tude this problem only arises at N^LO, since that is where two-pion-exchange 
diagrams first appear. A method to deal with the large non-analytic terms 



of Eq. (133) has yet to be uncovered, although there are hopes that the solu- 
tion will come from a better understanding of the renormalization of singular 
potentials. 

The phenomenology of KSW power-counting in the two-nucleon sector at 
relatively low orders has been extensively explored in a number of processes 
which are of current experimental interest. However, none of the calculations 
of these processes were carried out to an order sufficiently high that the strong 
N^LO correction to "^Si scattering played a role. Therefore, it may be that 
these results represent only low-order terms in a series which ultimately does 
not converge. None of the results obtained using KSW power-counting can be 
taken too seriously until the convergence issue is resolved. The fact that the 
KSW power-counting allows for all calculations to be performed analytically 
with dimensional regularization means that it is very easy to perform gauge 
invariant calculations for processes involving photons and^ther gauge fields. 
In particular, ed scattering has been calculated up to NLOlHI, jd —> jd Comp- 
ton scattering has beep co mputed (both for unpolarized and tensor-polarized 
targets) up to NL0E£3r£il, the polarizabilities of the deuteron have been com- 
puted to high-orders E£j, parity violation in np — s- dj radiative capture has 
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been explored 



and also parity violatioi; 
anapole moment has been investigatedliHi 



fd scattering due to the deuteron 
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Figure 29: Neutron-deuteron phase shifts in the S-wave quartet in LO (dashed), NLO (full), 
NLO in the pionless theory (dot-dashed) and N^LO in the pionless theory (dotted). Also 
shown are the results of PSA of neutron-deuteron (squares) and proton-deuteron scattering 

(stars). 



Calculations with perturbative pions have also been performed in the 
three-body sector E£3. In the J = 3/2, S-wave channel there are no unknown 
coefBcients appearing at NLO, and thus there is no dependence of the result on 
any fitting procedure. The LO calculation involves only two-body contact in- 
teractions. The NLO corrections come from one potential pion exchange, taken 
in perturbation theory. They are of two kinds: corrections to the deuteron 
propagator involving only two nucleons and true three-body diagrams, where 
the three-nucleons participate. The diagrams have to be computed numeri- 
cally since the half off-shell amplitude to{k,p) (solution of Eq. (p3|)) is known 
only numerically. As discussed before, the contribution from the tensor force 
vanishes in linear order in the S-wave so there is no contribution from the ten- 
sor part of the one-pion exchange. The results for the phase shift are shown 
in Fig. together with the extrapolation from the result obtained with the 
theory without pions and phase shift analysis, where it is available. The results 
indicate convergence up to very high energies but are barely distinguishable 
from the pionless theory. This is probably a consequence of the high anoma- 
lous dimension of the amplitude to{k,p), that decays to zero very rapidly at 
large p. 



7 Don't Mess with Texas: the Weinberg Program for Physics of 
the Two-Nucleon System 
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1.1 The Proposal 

One way of dealing with the non-perturbative nature of the pion exchanges 
uncovered by higher-order calculations kuKSW counting is to follow Weinberg's 
original proposal for the NN system 013. In this proposal one uses baryon 
xPT to compute the NN potential Vnn perturbatively. The S-matrix is then 
generated by the Schrodinger equation. In this approach, power-counting is 
not manifest in the amplitudes for individual processes, but there is a naive a 
priori estimate of the error in a particular calculation, which must be checked 
a posteriori. 

Weinberg's separation of all graphs contributing to some process in the NN 
system into graphs which are "two-particle reducible" and those which are not 
makes power-counting much simpler. Since, by definition, the "two-particle 
irreducible" (2PI) graphs do not contain the low-energy scales B (deuteron 
binding energy) and /M . The 2PI graphs are precisely those where the 
energy of all internal legs is of order . It follows immediately that the 2PI 
graphs have the power-counting of ordinary xPT. 

The "Weinberg program" for calculating properties of the NN system in 
xPT then centers around the construction of these 2PI graphs. Its implemen- 
tation has three stages. Firstly, we must consider the case where the NN 
system is not subject to any external probes, in which case the NN irreducible 
diagrams just yield the potential V (e.g. the left-hand graph of Fig. |3^.) Con- 
structing an NN potential using the rules of xPT will then yield the NN 
interaction, V'^"'\ up to a given order n in the chiral power counting. 

The second step is to obtain the full NN t-matrix by iterating these irre- 
ducible graphs together with the free two-nucleon Green's function, Gq: 

Go(i?;p) = — ^ -. (134) 

(77 is a positive infinitesimal.) This free Green's function can, of course, have a 
small energy denominator that does not obey the usual xPT power-counting. 
Iteration of V with Go gives contributions to the t-matrix such as the graph 



shown on the right-hand side of Fig. 30. This simply means that we follow the 



usual practice in the NN system, and construct the full NN t-matrix, T, from 
the chiral NN interaction and the free two-nucleon Green's function Gq, 

via the Lippmann-Schwinger (or, equivalently, the Schrodinger) equation 

r = "t^(")-f 1/(")GoT. (135) 



This point of view is similar to that seen when low-ene: 
tions are derived using Wilsonian RG arguments (see Refs. 



fective ac- 
IHZI'tHa for more 
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Figure 30: One- and two-pion exchange graphs. One-pion exchange is in the potential and 
(uncorrelated) two-pion exchange is generated by solving the Schrodinger equation. 

explicit investigations of this connection). These techniques are widely ap- 
plied, and have, for instance, been used to obtain improved actions for lattice 
gauge theory. In all cases the effective action is obtained using a perturbation 
theory in the ratio of low-energy scales to the Wilsonian RG cutoff scale. That 
Lagrangian is then solved to all orders. The error of such a calculation is ex- 
pected to be dominated by the lowest-dimensional operator omitted from the 
calculation. This allows one to roake an a priori estimate of the error, which 
can then be verified a posteriori^. In the same way the Weinberg program 
employs xPT to generate a nonrelativistic, particle-number-conserving Hamil- 
tonian for the nuclear system. This Hamiltonian is valid up to a given order 
in the chiral expansion and can be employed in any many-nucleon system at 
low energy. Once written down it should be solved to all orders. 

With the NN t-matrix in hand implementation of the third step of the 
Weinberg program can proceed. In this stage we consider the interaction of 
some external probe: a pion, a photon, or some weakly-interacting particle, 
with the NN system. Provided this probe carries momentum of order TOtt it 
can have its interaction with the NN system represented as a sum of irreducible 
diagrams which forms a kernel ifpiobc for the process of interest. Again, the 
power-counting of xPT applies to this kernel, since all nucleon energies ap- 
pearing in it are, by definition, of order to^. Naturally, this is only true if the 
external probe energy and momentum are both of O(to^), but, subject to that 
limitation, xPT allows us to calculate the kernel -ftTprobc up to a given order 
in the chiral expansion. The full amplitude for this reaction is then found by 
multiplying this kernel Xprobc by factors which take account of the initial and 
final-state interactions of the NN pair. In other words the amplitude is 



^= (l + TGo)ifprobc(l + GoT). 
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(136) 



Here T is calculated from Eq. (135), using a F^"^ which is also computed 
in xPT. For consistency the expansion for the interaction with the external 
probe should be taken to the same chiral order as the order n used in the 
chiral expansion of the NN potential 1/^"'. Once this is done the result for the 
amplitude A should be the unambiguous prediction of xPT for the process of 
interest up to that order in the chiral expansion. 

Up until now our discussion has been completely general, and could apply 
to the scattering of particles from the NN bound state, the deuteron, the 
break-up of the deuteron by a probe (e.g. — > np) or to the emission of 
particles from the NN state (e.g. proton-proton bremsstrahlung) . However, if 
we are especially interested in probes of the deuteron we may take the residue 
of Eq. (136) at the deuteron pole in both the initial and final-state. We then 
see that the amplitude for the interaction of the probe with the deuteron is 



^= (^lifprobcIV'), (137) 

where the wave function ji/;) is to be calculated by solving the Schrodinger 
equation for the potential T^*^"-' : 

{e-^^\^)=V^-^\^). (138) 

This program is a close relative of that pursued by nuclear physicists in 
the NN system for many decades. It centers around the same objects: poten- 
tials, wave functions and transition operators, and, like calculations using NN 
potential models, it treats the interaction of the nucleons non-perturbatively. 
The main difference to more traditional potential-model calculations of NN- 
system properties is the use of xPT to calculate the potentials and kernels that 
are the building blocks of the nuclear-physics calculation. Below we investigate 
this connection a little further, comparing and contrasting the xPT approach 
to the physics of the NN system with physics based on more phenomenolog- 
ical potentials. Crucially though, the use of xPT to derive the basic building 
blocks that go into potential models allows for simplification and systemati- 
zation of the calculations performed using the successful phenomenological 
potential models. 



7.2 The Weinberg Program: Power Counting 

The power counting for contributions to either the NN potential, V , or the 
kernel for some process, i^probc, is just that of xPT. In other words, we perform 
NDA on the kernels of processes, keeping track of factors of p and m^r in the case 
of V , and p, m,r, and k — the probe energy /momentum — in the case of Kp^ohe- 
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We generically denote all these "small" quantities as P. The expansion is then 
one in powers of P/ (A^, M). A diagram which contributes to V has an order 
n in this expansion which is calculated by just multiplying together the powers 
of P appearing in the diagram as follows: 

• Each vertex from counts as P"; 

• Each pion line counts as P^^; 

• Each loop counts as P^; 

• Each nucleon propagator counts as P^^. 

Note that these last two rules must be modified if we consider two-nucleon- 
reducible diagrams as well, since in such diagrams the nucleon propagator can 
scale either as P~^ or as P~^, depending on whether the energy present is 
"small" or of order mT^. 

For kernels A'probc one additional rule must be added: 

• Two-body contributions to A'p^obc carry an extra factor of P^. 

This rule, which means that two-body contributions to the reaction mech- 
anisms are systematically suppressed in this low-momentum regime, can be 
thought of as arising from the additional "loop" which occurs when a two-body 
reaction mechanism is sandwiched between NN wave functions. Equivalently, 
it occurs because there is one less momentum-conserving delta function in a 
graph with a two-body contribution to ATprobc- Either way the rule is again a 
straightforward consequence of NDA. 



7.3 Nucleon- Nucleon Potential 



The earliest discussion of constructing a NN potentialin the context of Wein- 
berg power counting was that by Weinberg in Refs. iM, where was com- 
puted. Higher-order calculations were carried out by Ordonez et al., where 
energy-dependent^iotentiaJs which represented the results for V^°^ , F , and 
were derived E£3'E£j'lHll. The energy-dependence of this potential was not 
intrinsic to ^ind could be transformed away. In Ref. E23 Epelbaum et 

al. did this, using unitary transformations to turn these energy-dependent 
potentials into energy-independent ones ^. By construction, this transforma- 
tion leaves all S-matrix elements unaffected, and so there is .no.. physical dif- 
ference between the potential of Ref. IH3 and that of Refs. P^H^^qjisil ^pj^g 

JTn Fjpf. Mj a, similar rnnncr.t'inn hptwppn energy-dependent and energy-independent forms 
of the two-pion-exchange NN force was established, thereby elitetalating the difference rfcm- 
tween old work on two-pion exchange by Brueckner and Watson [£3 and Taketani et ai. E£j. 
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olating NN potential can be derived from xPT in an analogous fasli- 
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but we shall discuss only the charge-symmetric, charge-independent 
piece of the NN interaction here. 

The Lagrangian of xPT is used to generate an NN potential, V as follows. 
The interaction V^"' contains all two-particle-irreducible NN NN graphs 
that are of chiral order n or .kpiver. At leading order (n = 0), the chiral NN 
potential has only two pieces ElEl: 

^(q) = -tI 2"^,"' (^1 ■ ^2) + Cs + Ct<Ji ■ a2. (139) 

where Cs and Ct are coefficients of contact interactions that do not involve 
spatial gradients. 

At any chiral order an analogous separation into a long- and short-distapx::£ 
NN EFT potential can be made. Thus we will follow the authors of Rcf. E2f 
and write V = V^^ + Kihort At low orders in the expansion Vshmt has a straight- 
forward counting in powers of P, with the higher-derivative operators naively 
suppressed by powers of P/A. 

Regardless of the scaling of the short-distance part, the long-distance piece 
is given by one-pion exchange plus irreducible two-pion exchange plus three- 
pion exchange etc. The dimension of the contributions of these multi-pion- 
exchange graphs produce in V^r increases by two for each additional pion loop, 
as can be seen by simple f^r counting. Consequently, irreducible multiple pion 
exchanges are suppressed by powers of 111^^/ A and p/A. Having made a chiral 
expansion for both V^hort and , we must truncate both at the same order in 
the chiral expansion parameter (p, m7r)/A. 

The contributions to V^hort at zcroth and second-order in xPT are just 
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(p'IKhortIP) = Cs+Ct<71-<72 

(P'iKhortIP) = + + (6q' + U^'){<Jl ■ CJ2) (140) 

+0-2) • (q X k) +C6(c^i • q)(o-2 • q) +^^7(^1 • k)(cr2 • k) 

where k = p' — p and k = (p' + p)/2. The translation into the short-range 
potential that is operative in a particular partial- wave is then a straightforward 
matter of taking matrix elements of these operators. We will not present the 
explict formulae here, but two points are immediately apparent: 

1. This short-distance potential operates only in S-waves, P-waves, and to 
mix the ^Si and ^Di states. 

2. In the and ^Si channels there is a LO and NLO contribution, with 
the NLO piece being proportional to {p^ + p'"^) by Hermiticity. 
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3. In the next five partial waves {^Pi, ^Pq, ^Pi, ^P2, and ei) there is a 
straightforward set of linear equations which relates the combination of 
short-distance parameters that appears in any one partial wave to the 
parameters appearing in the potential. 

This combination of facts means that by performing a fit with one param- 
eter in the partial waves ^Pi, ^Pq, '^Pi, '^Pj, and ei, and two parameters in 
the partial waves ^5*0 and ^Si we may completely determine the short-distance 
potential. 

At NLO the long-distance, pionic part of the potential is generated by the 
one-loop diagrams shown in Fig. These diagrams are divergent in xPT 
and require renormalization, but the required renormalization may be made 
by absorbing the divergences in the constants — > ^7 given above. Strictly 



speaking then, the relationship of the ^1 ^7 which appear in Eq. (141) 
to the parameters which appear in the Lagrangian involves both the linear 
transformation mentioned above and this subtraction of certain forms arising 
in the calculation of the graphs shown in Fig. pl^. After renormalization, the 



(a) (b) 
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Figure 31: The two-pion-cxchange graphs which are included in the A'^A'^ interaction at 

NLO. 



long-distance part of V^^ at NLO is 
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(p'|yi^)|p) = -L(|q|)-^[(ai.q)(a2-q)-q^(ai.a2)] (141) 
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L(|q|) [Aml{bg\ - Agl - 1) + ci^{23g\ - 10g\ - 1) 



'384^2^4 



4to2 + q2 
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with 



L(|q|) = i^yi;^ 



log 



V4to2 + q2 + |q| 



2to, 



(142) 



At N^LO, 0{P^ / A^), there are no new short-distance operators. There 
are, however, new contributions to Vtt at this order. Where the Weinberg- 
Tomazawa term entered in the "triangle" and football graphs (diagrams (c) 
and (d) of Fig. and contributed to V^"^^ we can now make one insertion of 

— (2) 

an operator from the pion-nucleon Lagrangian C). ^ and generate contributions 



to 



(3) 



+ 



with 



(3) 

After renormalization Vv reads 
(p'l^i^'lP 



35i 



167r/4 [ 16Af(4m2 +q2) 



9l 



2to^(2ci -C3) -q 



C3 



16M / 



(2m^ + q2)A(|q|) 



5127rM/; 



4mJ + q2 

+ (4m2 + 2q2 - (4^2 _^ 3^2^) (3^2 + q2)A(|q|: 
((ai • q)(a2 • q) - q2(ai • as)) (2m2 + q2)A(|q|) 



3^ 



32^/4 



(n • fs) (((Ti • q)(cr2 • q) - q^(CTi • 0-2)) 



C4 



— ^ (4m2+q2)_^(10™2+3q2) 



4A-f y 

■ i (ai + <J2) • (p ' X p) {2ml + q')^(|q|) 



^(|q|) 
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4l^^(ri • r2) ^ (ai + ^2) • (p ' x p) (4m^ + q^)A(|q|), (143) 



^(|q|) = 



2™^ 



(144) 



The free parameters ci, C3, and C4 are related directly to data in the ttN 
system. _Xlae_p's have been determined in recent fits to -kN scattering data 
in xPT t23>lll3. In principle then, there are no undetermined parameters in 
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("3") 

V^TT • The fit to the data extrapolated to the region inside the Mandel- 
stam triangle is used by Epelbaum et al. to fix ci, C3, and C4. Interestingly, 
Sugawara and Okubo wrote down a TPE potential which-cpntained a large 
part of the correct TPE presented above as early as 1960 tld. Stoks andrEi,- 



jken also generated all of the two-pion exchange interactions written abovelll3, 
although their potential also contained many other, more phenomenological, 
short-range contributions. Although the fact that xPT connects nN data to 
NN data in a mqd pj|-i njd pnen d ent way makes it similar in philosoplw-to these, 



and many other tl2 
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!d3, works, it was really only in Refs. 



L3C 
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that a 



complete expression for the full chiral TPE was first written down. 

This chiral TPE potential was also calculated by Kaiser et al. as part of 
a study we will discuss in more detail below til. Kaiser et al. observed that 
the combination V'^'^^ + V^'^^ provides strong scalar-isoscalar attraction. This 
attraction is governed by the LECs ci and C3 . This suggests that if two-pion 
exchange is calculated in a fashion consistent with chiral symmetry then there 
may be no need for the presence of an explicit cr-meson in the NN potential. 

This potential 
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must now be iterated in the Lippmann-Schwinger equation, Eq. ( 135 ). How- 
ever, it is a poorly-behaved potential from a quantum mechanical point of view, 
since it is very singular at short distances, and the corresponding Hamiltonian 
is thus unbounded from below. This is perhaps not surprising, since we know 
that the short-distance physics in xPT is not the full short-distance physics of 
QCD. At distance scales ~ 1/A^ mechanisms not encoded in xPT will begin to 
play a dynamical role, and will regulate the behavior of. the singular potential 



y ['^1 . In order to calculate with the potential Refs. introduced cutoffs 



on the potential before iteration, defining a regulated potential 

(p'I^«Ip) = /a(Ip'I)(p'I^Ip)/a(|p'|), 



(146) 



where /a(|p|) is some regulator function that goes to zero sufficiently fast as 
IpI goes to infinity, and obeys /a(0) = 1, e.g.. 



/A(|p|)=exp(-|p|VA4). 



(147) 



This regulates the potential 
then quite simple. Wj-th. 
distance coefficients t3ll£3l2i 
observables. 
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Practical implementation of renormalization is 
cegulator mass. A, fixed, one tunes the short- 
so as to reproduce some specified low-energy 



71 



With renormalizable theories, the next step would be to take the regulator 
to infinity. Here A should be of the order of the high energy scale A^. There 
are two reasons for this. Firstly, NDA applies when A ~ A^. Secondly, in 
the Weinberg scheme subleading terms are resummed to all orders and may 
become large for large A as the counterterms required to renormalize these 
contributions are not present. 

The algorithm for generating predictions from the regulated third-order 
chiral potential adopted in Ref. t23 is then: 

1. Choose a regulator mass A. 

2. Fit the constants Cs, Ct, andpfi-^y to the S and P-wave phase shifts 
and ei from the Nijmegen PSAS up to Tiab = 100 MeV. 



3. Predict these same phases from Tiab = 100 MeV to Tjab 
predict other phase-shifts over the full energy range. 



300 MeV, and 





Figure 32: Predictions for the ^Sq and phase shifts from Ref. E3. The dotted, dashed, 
and sohd curves represent the results from the LO, NLO, and N^LO potential. At NLO and 
N^LO the two parameters in this channel are fitted to the phase shifts below Ti^b = 100 MeV. 
The squares represent the Nijmegen PSA phases. The exponential regulator was employed 
with A = 500 MeV at NLO, and A = 875 MeV at N^LO. 



In general the results for the S and P waves (where at least one fit pa- 
rameter was present in each case) are quite good, although the behavior of the 

phase shift causes some concern. Some of these results are displayed in 
Figs. |3^ and 

We stress once again that there are no free parameters to fit in partial- 
waves with L > 2. Thus, the phase shifts calculated there are parameter-free 
predictions of this approach. In the D-waves, the results are in quite good 
agreement with the Nijmegen PSA although some deviation begins to be seen 
around Tiab = 200 MeV (see Fig. 35). Epelbaum et al. also note that these 
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Figure 33: The ^Pq a-nd phase shifts as calculated in Ref. In this case there is one 
parameter at NLO and N^LO, which is fitted to thenhase shifts below Ti^b = 100 MeV. 

Legend as in Fig. p2l 




0.05 0.1 0.15 0.2 0.25 0.3 0.05 0.1 0.15 0.2 0.25 0.3 



Figure 34: The mixing parameters ei and 62 from Ref. There is one fit parameter in 
ei. The results for Ti^b > 100 MeV in ei, as well as aM of the results for €2 are predictions. 

Legend as in Fig. p2l 
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phases are quite sensitive to the value of the cutoff chosen. This may be because 
at next order, in V'^^\ contact interactions enter which will affect the D-wave 
NN channels. Meanwhile, the mixing parameter £2 is not well reproduced 
above Tjab = 100 MeV, while the results for the triplet F-waves are really 
quite poor. It seems that for these smaller NN phase shifts a higher-order 
calculation is needed for an accurate description. 





Figure 35: The ^Di and ^D2 phases from Rcf. 
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In this case there are no free parameters. 
Legend as in Fig. p2l 



As the orbital angular momentum is increased the Born approximation 
works increasingly well, and there is, in fact, little need to iterataihe potential. 
In this way the study of Epelbaum et al. reduces to that of Ref.EII in the higher 
partial-waves. In that work Kaiser et al. calculated the following perturbative 
NN amplitude 



(148) 



In other words, these authors employed the potentials V^^ and at tree 
level, and OPEP in second Born approximation (i.e. iterated once) as the 
iViV amplitude. The graphs included are those of Figs. ^ and ^ together 
with the triangle and football graphs with insertions from C^^j^. This is an 
approximation to the full NN t-matrix obtained by iterating the chiral NN 
potential V^^'^-' . It works increasingly well as the angular momentum of the NN 
state increases, since the centrifugal barrier ensures that only the long-range, 
weak, pact of theJVA^ force is probed in these waves. Comparison of the results 
of Refs.tZl andt^S indicates that iteration begins to be important below L — 4, 
but that for i > 4 the amplitude ( |148[ ) does a good job of reproducing the 
phase shifts. 

Indeed, in these peripheral partial waves one-pion exchange alone, i.e. 
the interaction T^i*^^ uniterated, reproducer the experimental data reasonably 



well — as has been known for many years tla. For L = 6 the effects of the parts 
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Table 2: Results for the PWA with different long-range interactions and R = 1.4 fm. 





^Parameters 


• 

A mm 


OPE 






31 


2026.2 


OPE H 






28 


1984.7 


OPE H 


h - 




23 


1934.5 



of Tport other than V^i"' are really quite difficult to see below Tiab = 200 MeV. 
However, for the G and H-waves the inclusion of V^^'' and V'^^'' improves the 
agreement with the data. 

A calculation of the next-order piece of the chiral NN potential, y^'*-' will 
bring two new ingredients. Firstly, a number of short-distance operators will 
arise. The role that these might play in the ^S'o channel has been considered 
in Ref.tl2l. In general though, a number of channels will have to be considered: 
these operators will have to be fitted to the S, P, and D wave phase shifts, as 
well as to mixing parameters. Secondly, many two-loop .contributions to the 
potential — some of which have been calculated in Refs. P°lr^H — will need to 
be added to the interaction. The calculation of the full F'-^-' is a Herculean 
task. Nevertheless, we can get some feeling for the size of higher-order effects 
by including more freedom in the short-distance part of the interaction. 

This is essentially the direction taken in a recent calculation by the Ni- 
jmegen grou pE£3 . Rentmeester et al. calculated the piece of the chiral potential 
of Eq. (145) which acts at radii r > — i.e. the potential's long-range 
tail — and incorporated it into the Nijmegen PSA machinery. Schematically, 
they worked with a potential 



V{r) =Vsi,ortir) e{R-r) + 



Li=0 



-R), (149) 



where Vera accounts for electromagnetic corrections to the long-distance part 
of the potential. As for the short-distance part, Kshort is only operative at 
radii r < R, and it depends on many free parameters. These parameters are 
tuned so as to reproduce the experimental phase shifts. Fixing the boundary 
R at 1.4 fm, and the constant ci = —0.76 GcV^^, Rentmeester et al. tuned 
23 parameters in T^hort and the ttN LECs C3 and C4 to fit 1951 pp scattering 
data below Tiab = 350 MeV. The resultant minimum is very close to the 
theoretical ideal of one per degree of freedom (see Table ||). 

This is a strong indication that the chiral TPE is the intermediate-range 
part of the NN force. This conclusion is reinforced by what happens if the 
chiral TPE is removed from the interaction. In that case the minimum 
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shoots up appreciably, even if a number of new parameters are introduced in the 
short-distance interaction. The final indication that the chiral NN potential is 
reproducing the physics of the NN system in this phase-shift analysis is that 
the values of C3 and C4 obtained from this fit are C3 = —4.99 ± 0.21 GeV^^ 
and C4 = 5.62 ± 0.59 GeV~^. These, together with the value of ci, agree with 
the values extracted for these LECs from the irN data, within the significant 
uncertainties associated with fitting in the ttN system. 



8 External Probes of the Deuteron 



In scattering processes involving the deuteron EFT has proved more successful 
phenomenologically than traditional potential models. Part of this success is 
due to the fact that the relation between exchange currents and potentials is 
dictated by the rules of field theory and is therefore unambiguous in EFT. Per- 
haps more importantly, EFT relates scattering processes involving the nucleon 
alone to nuclear scattering processes in a controlled manner; all of the phe- 
nomenology of xPT in the single nucleon sector makes its appearance at some 
level in the nuclear sector. Hence, counterterms which are difficult to measure 
in experiments with nucleons can sometimes be extracted in a systematic way 
from nuclear experiments. This is particularly important given the absence 
of suitable neutron targets. Weinberg power-counting has led to fruitful caaa. 
putation of many pionic and photonic probes of the two-nucleon system 
including electron-deuteron scattering, electromagnetic deutesan, form factors 
neutral pion photoproduction on the deuteron ailhreshold 



electroproj 
ing lengthJ 
capture B 



£±ion on the deuteron at threshold 



L57 



155 



15C 



158 



neutral pion 
pion-deuteron scatter- 
is-well as pn radiative 



Hj. We will consider 



Compton scattering on the deuteronE 
and the solar burning process pp — > de^v 
each of these processes in turn. 

On the technical side, it should be noted that typical nucleon momenta 
inside the deuteron are on the order of m^, and consequently, we expect no 
convergence problems in the xPT expansion of any low-momentum electro- 
magnetic or pionic probe of the deuteron. We will use wavefunctions gener- 
ated using modern nucleon-nucleon potentials as well as wavefunctions recently 
computed in xPT . Generally we find that any wavefunction with the correct 
binding energy gives equivalent results to within the theoretical error expected 
from neglected higher orders in the chiral expansion. 
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8.1 Electron- Deuteron Scattering 

For simplicity we will work in the Breit frame in calculating matrix elements 
of the form, 

(j q/2|J,'',Jz -q/2) (150) 

where a virtual photon of three-momentum (\ = qz couples to a deuteron 
electromagnetic current when the deuteron is in states with specific magnetic 
quantum numbers i and j, as discussed earlier. The electromagnetic structure 
of the deuteron may be parametrized in terms of the three form factors, Fc{q'^), 
FMiq^), and Fqiq^) defined in Eq. (||). 



8.2 The deuteron current 

In general, the current to which the electron couples in electron-deuteron scat- 
tering may be expressed as: 



(P' pVJP P> = (27r)3jW(p,q)<5(p'-p 



q/2)+jW(p,p';q), 



(151) 



where p and p' are the initial and final relative momenta in the deuteron state, 
and q = P' — P is the three-momentum of the photon exchanged between the 
electron and the deuteron. The one- and two-body currents discussed here 
are depicted in Fig. ^ This is a useful decomposition because Weinberg 
power-counting predicts that the two-body currents are generically suppressed 
relative to one-body currents. 




Figure 36: One- and two-body current contributions to electron-deuteron scattering. The 
square and rectangle indicate, respectively a -yNN vertex, and the two-body NN — > NN-y 
kernel. The blobs are deuteron vertex functions. 

For on-shcU nucleons the most general form that the single-nuclcon current, 
can take, consistent with Lorcntz invariance, time-reversal invariance, and 
current conservation, is: 



J^Hp,q) =eAru(p + 3q/4) 



Fi{q'h,. + F2{q')^<J. 



u(p-q/4), (152) 
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where cat is the nucleonic charge. For the deuteron case, Fi and F2 are isoscalar 
combinations of the proton and neutron Dirac and Pauh form factors (i = 1, 2) 



E,iq') ^ F^^\q') + F^-\q') . (153) 
It is straiehtforward to calculate these isoscalar Dirac and Pauli form fac- 



tors in xPT l^b"H^b^H^b^ At LO the only contribution comes from C'^J'jj, and 



corresponds to an ^0 photon coupling to the charge of the nucleon. At NLO 

(2) 

the current receives contributions from the El and Ml vertices of 
Thus the first two orders of the calculation are trivial, and lead to the cur- 
rents: 

J^'^(p,q) = e , j«(p,q) ^ e (|± + ^(«., + + 1 w) .(154) 

Corrections at next order will include relativistic effects, and the effects 
of nucleon structure. The former are suppressed by ratios of momenta to the 
nucleon mass, which we denote generically by the mnemonic: /M"^ , p ■ 

Q/M'^, q^/M^. Corrections at order 6^ include thejisual spin-orbit and Dirac- 



Foldy relativistic corrections to the charge operator[L22l. Meanwhile, the second 



expansion will be governed by the parameter g^/A^. In the isoscalar form fac- 
tors these structure corrections are purely due to short-distance counterterms, 
since a photon cannot couple to two pions in an isoscalar state. 

We now turn to the task of quantifying the effect of two-body currents 
on the observablea under consideration. The power counting here was first 
explained in Ref.E3 and has beep-further worked out for the exchange currents 
which contribute to Fm in Ref. 
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The contribution of two two-body pionic mechanisms to the deuteronic 
current is depicted in Fig. It is straightforward to count the powers of 
momentum or pion mass in such graphs. These current contributions scale 
at least as P"^ |. Therefore we expect that such meson-exchange currents will 
be suppressed by P^/A^ relative to the lowest-order ^NN interaction, which 
scales as P. 

Since j+ begins only at 0{P^) it might seem that meson-exchange currents 
are only suppressed by one power of P in the deuteron form factor Fm , but this 
is not the case. In fact, although the leading effect of meson-exchange currents 
is nominally O(P^), the meson-exchange current which occurs at that order 
is isovector in character and thus vanishes in the deuteron. It is well known 
that this 0{P^) isovector current is the leading MEC correction to processes 
such as electrodisintegration of the deuteron, and 7c? — > np (see, for instance, 



'^Here we are counting the charge e as of order P in order to simplify matters. 
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<ri-<±i- 

Figure 37: The contribution to ed scattering from the LO pionic current. Here all vertices 
are from the LO chiral perturbation theory Lagrangian. The graphs with nucleons one and 
two interchanged also contribute. 



Ref.E£3'E3). In practice this results in the leading effect from meson-exchange 
currents coming at 0{P^), or N'^LO in the deuteron form factors. In such 
an 0{P'^) exchange current the -fir contact interaction with the nucleon has 
a coefficient that is related to known constants by the requirements of gauge 
invariance and reparametrization invariance. We note that this correction 
to the charge operator is phenomenologically important in potential-model 



1 
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calculations of Fc 

At some order above 0{P^) operators enter whose coefficient is not con- 
strained by symmetries, but instead is sensitive to short-distance physics. 
These have only a small effect in Fc and Fm for 1 GeV, but they may 
affect Fq markedly, since the natural size of Fq at q = is much smaller than 
that of Fc and Fm- A natural expectation for ^qI^^q = Qd is: 

Qd-^. (155) 

T 

Here rjsd — the asymptotic D-to-S normalization — enters because an S-wave 
deuteron cannot emit a quadrupole photon via any one-body current mech- 
anism. Thus the power-counting for the contribution of one-body currents to 
Fq is exactly the same as for similar contributions to Fc- However, two-body 
currents can connect the S-wave pieces of the deuteron initial- and final-state 
wave functions one to another and still result in the emission of a quadrupole 
photonEj. The lowest-order Lagrangian which can generate this effect is that 
written down in Eq. (|56|). This two-body operator does not enter in the ex- 
pansion of the deuteron current until O(P^) in Weinberg counting. On the 
other hand, crucially it contains no powers of rjsd- Hence, numerically it can 
be important at much lower order than naively expected. Counting only over- 
all powers of p, M, and A, we see that the contribution of this operator to Qd 
will be of order 
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A/A2 A2' 



(156) 



Although this is formally of order P^, numerically it is as important as N^LO 
corrections to Q^. 
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From Ref. 



131 



138 



we have both NLO and N^LO NN poten tials , derived 
from chiral perturbation theory. The Schrodinger equation, Eq. (138), can be 
solved with these potentials to generate deuteron wave functions. We can then 
use the systematic chiral expansion of the deuteron electromagnetic current 
to calculate the electromagnetic form factors of the deuteron. These results 
are compared with results from a sophisticated potential model, the Nijm93 
potential. In both cases we employ the nonrelativistic impulse approximation 
with point nucleons, since this represents a calculation to NLO in the ed scat- 
tering kernel. Before calculating Fq, Fq, and Fm though, we will spend some 
time discussing results for static properties of the deuteron. 



8.3 Static Moments 

The values of the form factors Fm and Fq at — are related to the magnetic 
and quadrupole moments of the deuteron via the definitions in Eq. (|5^). If 
no meson-exchange currents (i.e. no local four-nucleon-one-photon operators) 
enter the ed scattering kernel, and there are no relativistic corrections to the 
current, then fid and Qd are given by simple integrals of the deuteron radial 
wave functions u{r) (S-wave) and w{r) (D-wave), 

Aid = Hp + Kn - -{Kp + Kn - -) I dr {r); (157) 







-j^ /"OO Y /"OO 

fjLQ — / dr r"^ u{r)w{r) — — / drr'^w'^{r). (158) 

V50 Jo 20 Jo 



Another interesting static property of the deuteron is its matter radius, 
rd- This is dominated by the long-distance part of the deuteron wave function, 
since it occurs in an integral weighted by r^: 

drr^[u^{r) +w^{r)]. (159) 

F,m ploying^ the wave functions derived above, along with three leading po- 
tential models auu leads to the results shown in Table |3[ The experimental 
data is included for comparison only, since we have used a very simple electro- 
magnetic current to generate these results. 
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rd (fm) 




MS (fm^) 


Nijm93 


1.968 


0.847 


0.271 


AV18 


1.967 


0.847 


0.270 


CD-Bonn 


1.966 


0.852 


0.270 


XPT NLO 


1.975 


0.859 


0.266 


XPT N^LO 


1.967 


0.845 


0.262 


Experiment 


1.971(6) 


0.857406(1) 


0.2859(3) 



Table 3: Static moments of the deutoron for the Nijm93, CD-Bonn, andpWlS potentials, as 
well as for the NLO and N^LO xPT potentials of Ref.ESl. 



The table shows that the deuteron matter radius is the same to three 
significant figures in all six calculations. Moreover, the results are all within 
the error bar pLihe, value extracted from hydrogen-deuterium isostope-shift 
For the magnetic moment the situation is not quite so 



16^ . 



163,169 



measurements 

good, but it is clear that Hd is predominantly determined by Kp and The 
size of the spread in tha.±iieoretical models — and of their discrepancy with 
the experimental value e3cj — is consistent with the magnitude expected for 
higher-order contributions to the deutero n cu rrent. 

Park et al.tM and Phillips and Cohent£ZI have also performed calculations 
of the static properties of the deuteron in EFT-motivated approaches. Both 
of these papers employed one-pion exchange as the long distance part of the 
NN potential — rather than chiral two-pion exchange — and various different 
short-distance regulators. The results are in accord with the above findings. 
Firstly, is extremely well-predicted in the effective theory, once As, the 
asymptotic S-state normalization of the deuteron wave function, and B, the 
deuteron binding energy, are reproduced. This is in line with the findings of 
EFT(/r) and ERT. The deuteron radius is largely determined by the tail of 
the deuteron wave function. Similarly, fid is largely insensitive to the short- 
distance physics, since it is dominated by Kp and k„. 

The large variation in Qd bctweeii the_calculations presented in Table || 
was also seen in the studies of Ref. E£j andE£3. There is considerable sensitiv- 
ity to the short-distance physics in this observable. Presumably, this is why 
the five theoretical values all differ markedly from the experimental number. 
Apparently there are significant contributions to this quantity from physics be- 
yond those of single nucleon operators. One such piece o£x)hysics is a two-body 
current resulting in the emission of a quadrupole photonLJ. As argued above, a 
natural size for the contribution of this operator to fiQ, will give approximately 
as large an effect as terms one order in the expansion beyond what we have 
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considered here. 



8.4 Electromagnetic Deuteron Form Factors 



Fc calculated with the NLO wave function of xPT is shown in Fig. ^ It is 
compared with results from the Nijm93 wave function. The agreement is good 
for photon momenta of up to ab out 500 MeV for a regulator mass A = 600 MeV 
in the exponential of Eq. (147). 



10" 



Nijm93 
NLO xPT 
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Figure 38: The charge form factor of the deuteron at NLO. The sohd line is the result 
for the Nijm93 potential, while the dot-dashed line is a calculation with the NLO chiral 

potential. 

Results are presented for the deuteron magnetic form factor Fm in Fig. 
These plots lead to conclusions similar to those inferred from the plot of Fc- 
This is a consequence of the dominance of the S-wave pieces of the wave func- 
tion in determining the magnetic structure of the deuteron at low q. 

The situation is somewhat different for the quadrupole form factor Fq. 
Since the quadrupole moments obtained in the Nijm93 and NLO xPT calcu- 
lations do not quite agree there is a discrepancy of a few percent between the 
two form factors at g = 0. On the other hand, the shape of the two curves 
is essentially the same out to g « 500 MeV. In fact, the quadrupole moments 
in both calculations differ appreciably from the experimental result. Simple 
power-counting arguments suggest that a two-body counterterm might be im- 
portant in Fq. To further examinethis problem, we have calculated the ed 
tensor-polarization observable T20 in four different ways. We have performed 
two different types of calculation for each of the two wave functions: Nijm93 
and NLO xPT: 

• A calculation with single nucleon operators; 



10 



10"' 



82 



10° 



10" 
10"' 



Nijm93 
NLO xPT 



Nijm93 
NLO xPT 



0.0 200.0 400.0 600.0 800.0 

q (MeV) 



10 



0.0 200.0 400.0 600.0 800.0 

q (MeV) 



Figure 39: The magnetic and quadrupole form factors of the deu tern n for the Nijm93 and 
NLO xPT wave functions. Legend as in Fig. tSq. 



• A calculation where we add — in an admittedly non-systematic way — a 
two-body counterterm which shifts the value of the quadrupole moment 
in both calculations to the experimental value. 

The results of this investigation are presented in Fig. We see that 

when the same currents are used the results for T20 are the same with both 
wave functions out to g « 500 MeV. This is just a restatement of the results 
presented for electromagnetic form factors above. However, we also see that 
shifting the value of Qd from the impulse approximation result Qd ~ 0.27 fm^ 
to the experimental value Qd = 0.29 fm^, has a marked effect on T20 at momen- 
tum transfers of order a few hundred MeV. Moreover, it would seem that once 
the deuteron quadrupole moment is rpprndnrpd hnth thp pntpntial mnHp1 and 



the EFT calculation describe the data of Refs.liZJ'EJiliiilJ^ESLia fairly well 



out to momenta q w 600 MeV. This suggests that the physics which determines 
T20 at these momentum transfers is the physics of one-pion-exchange — which 
the two potentials here have in common — along with a deuteron which has the 
right binding energy and the correct quadrupole moment. 



8.5 7d->7r"(i 

Pion photoproduction on the nucleon near thresholiLhas been studied up to 
0{P^) in with the delta integrated out in Ref. The differential cross- 
section for a photon of momentum k to produce a pion of momentum q at 
threshold is pure S wave and given by 

'Thanks to Matt Dorsten for his help in generating this figure. 
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Figure 40: The ed scattering tensor polarization observable T2o- The dashed (dotted) 
lines correspond to insertions of a single nucleon operator only, using the Nijm93 (NLO 
xPT) wave function. The solid (dot-dashed) line is the Nijm93 (NLO xPT) result when a 
two-body counterterm is added so th at r(|pro^ue^d. [Experimental data is taken from 

Refs.l 
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where i?o+ is the electric dipole amphtude. At LO, 0{P), there are contribu- 
tions only to charged pion channels arising from minimal substitution on the 
pion-nucleon vertex ( "KroU-Ruderman term") and from the pion pole. There- 
fore the charged channels are expected to be large compared to the neutral 
channels. At 0{P^), minimal substitution on the relativistic correction to the 
pion-nuclcon vertex contributes to photoproduction on the proton. At 0{P^) 
there are tree-level contributions involving the magnetic moment of the nu- 
cleon. Together with the previous orders, this reproduces an old "low-energy 
theorem". However, at 0{P^) there are also finite loop corrections where the 
photon interacts with a virtual pion which then rescatters on the nucleon. This 
large quantum effect was missed in the old "low-energy theorem" and in most 
models. At 0{P^) there are loops with relativistic corrections at the vertices 
together with new counterterms. ResultslHS for Eo+ up to 0{P^) from a fit 
to data above threshold constrained by resonance saturation of parameters are 
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shown in Table l'^ 



Table 4: Values for -Bo+ in units of 10 ^/m^+ at various orders in and in experiments, 

for the different channels. 





-(io--V"^-+) 


P 


pi 


pi 


pi 


Experiment 


IP - 




34.0 


26.4 


28.9 


28.2 


27.9 ± 0.5 
28.8 ± 0.7 
27.6 ± 0.3 


jn 


— > ■K~p 


-34.0 


-31.5 


-32.9 


-32.7 


-31.4 ± 1.3 
-32.2 ± 1.2 
-31.5 ± 0.8 


IP- 







-3.58 


0.96 


-1.16 


-1.31 ± 0.08 
-1.32 ± 0.08 


7n 










3.7 


2.13 





Charged-pion channels are indeed larger than neutral pion channels and 
show good convergence and agreement with experimental data. For neutral- 
pion channels convergence is less apparent, but the absolute values are much 
smaller. The O(P^) result for 7p — > n^p is in relatively good agreement with 
the recent results from Mainz and Saskatoon. To the same order there is 
a prediction for the near-threshold behavior of the 771 ir'^n reaction; the 
cross-section is considerably larger than that obtained in models that omit the 
important one- loop 0{P^) diagrams. In fact, the neutron cross-section is a 
factor of four larger than the proton cross-section, in violation with classical 
intuition. Neutral pion photoproduction on the deuteron is the best candidate 
for a test of this peculiar prediction for the neutron. 



We now consider pion photoproduction on the deuteron Consider a 
process with two nucleons, a single photon and a single pion in the initial and 
final state, respectively. Following the counting rules outlined above, the con- 
tributions are ordered as follows: 0{P): At LO, the matrix element is given by 
tree graphs with one spectator nucleon, with LO vertices, 0(P). This contribu- 
tion is nonvanishing only for charged pion photoproduction. At this order one 
finds the KroU-Ruderman term (Fig. |4l|a). 0{P^): The first corrections to the 
leading terms are the same as above but with 0{P^) vertices. For example, 
there is a 1/M correction at this order (Fig. |4l|b). 0{P^): There are three 



'"In this calculation, the dominant isospin breaking effect, namely the charged to neutral 
pion mass difference has been taken into account. A fully consistent calculation including all 
effects from virtual photons and the quark mass differences is not yet available. The result 
for Eq^^" is based on the assumption that the counterterms entering at 0{P'^) are the same 
for the proton and the neutron apart from trivial isospin factors. 
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classes of corrections at this order: (i) One-loop graphs with one spectator nu- 
cleon (Fig. |4l|c). (ii) Tree graphs with one spectator nucleon and with either 
one 0{P^) vertex (Fig. |4l|d) or two 0{P^) vertices. In the case of neutral 
pion photoproduction, the counterterm graphs with one O(P^) vertex vanish, 
while graphs with two 0{P^) vertices give rise to corrections propor- 

tional to the nucleon magnetic moments, (iii) Tree graphs with an interaction 
between nucleons with only 0{P) vertices. These graphs are of two types: two- 
nucleon graphs like the Feynman graph of Fig. |4l| e and timc-ordcrcd graphs 
such as those in Fig. p] f. 0{P^): There are many new contributions at this 
order, which can be constructed following the same procedure. Most of these 
are calculable corrections to the lower orders, but among the single-scattering 
diagrams for neutral pion photoproduction an undetermined counterterm ap- 
pears. 




Figure 41: Characteristic graphs which contribute at each of the first three orders in pion 

photoproduction. 

In the case of neutral pion photoproduction at threshold a number of the 
graphs in Fig. ^ will not contribute: those where the photon line is attached 
to a pion, and those that go like 5* • where q is the outgoing pion momentum 
and S is the covariant nucleon spin. In particular, all the leading-order graphs 
(Fig. ^a) vanish, which immediately suggests that the cross-section will be 
smaller than that for charged pion production, and more sensitive to two- 
nucleon contributions. Moreover, the time-ordered graphs (Fig. also 
vanish. It should be noted that contributions from the delta are minimized at 
threshold, and so integrating out the delta should not damage convergence of 
the theory too much. As noted above, to O(P^) there is one undetermined 
parameter, which is fitted to 7p n'^p above threshold. 
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We limit ourselves here to threshold neutral-pion photoproduction on the 
deuteronE£3. As explained in the introduction, this is of special interest because 
it can provide information on the elementary neutron amplitude. 

Consider the reaction '-f{k) + d{pi) tt'^ (q) + d{p2) in the threshold region, 
(f ~ 0, where the pion is in an S wave with respect to the center-of-mass (cm) 
frame. For real photons {e-k — 0, with e the photon polarization vector), the 
differential cross-section at threshold takes the form 



|fc| da_ 
\q\ dil 



= 7:El (161) 

l9l=0 



We now present the results of the chiral expansion of the dipole amplitude Ed 
to 0{P'^). The single-scattering contribution is given by all diagrams where the 
photon and the pion are absorbed and emitted, respectively, from one nuclcon 
with the second nucleon acting as a spectator (the impulse approximation). 
Some of these diagrams are depicted in Fig. To 0{P^) one finds 

l + m^/md [2 ^ ^ J 

^ k - I d^pp]-{pf^ + pf^^)<t>*{p)t- Uip- k/2)\, (162) 



M 7 2 
evaluated at the threshold value 

\k\ = fcthr = m^o (l - = 130.1 MeV. (163) 

V 2md / 

Here is the deuteron mass, J = (cti + (T2)/2, and is the deuteron wave- 
function. 

Two remarks concerning Eq. ( |162| ) are in order: 

(1) The elementary S'-wave pion production amplitudes to 0(P^)llZ3 are 
given in Table IJ. The proton amplitude that results from a one-parameter 
fit above threshold compares well with experiment. The neutron amplitude is 
predicted to be considerably larger than the one obtained from the old "low- 
energy theorem" that omits the important one-loop O(P^) single-scattering 
diagrams (Fig. ^l|c). 

(2) It is important to differentiate between the ir^d and the n^N {N = p, n) 
cm systems. At threshold in the former, the pioa is-produced at rest; it has, 
however, a small three-momentum in the latterll23'tZ3. Consequently, a single- 
nucleon P-wave contribution appears which is proportional to the elernentary 
amplitudes P^ ^ and P^°" that obey the P-wave low-energy theorems 
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^1 



0.480 \q\ GeV" 



0.344|g|GeV- 



(164) 



with 



M ) M 



2M 



4M 



(165) 



where p is the nuclcon momentum in the ir'^d cm system. The P-wave am- 
phtude for neutral pion production off protons holds to within 3%. Here we 
neglect the energy dependence of the elementary and S- and P-wave 
amplitudes since the pion energy changes only by 0.4% for typical average 
nucleon momenta in the deuteron. 

The two-nucleon contributions at 0{P^) that survive at threshold in the 
Coulomb gauge are shown in Fig. All vertices come from the lowest order 
Lagrangian. 




Figure 42: Two-nucleon graphs which contribute to neutral pion photoproduction at thresh- 
old to 0{P^) (in the Coulomb gauge). All vertices come from 



At 0(P*), we have to consider tree graphs with exactly one O(P^) ver- 
tex. The relevant terms in this Lagrangian fall into two classes, the first with 
fixed coefficients due to Lorentz invariance, the second with four low-energy 
constants ci , . . . , C4 that are related to pion- nucleon scattering and two other 
constants that are related to the isoscalar and isovector nucleon anomalous 

magnetic moments. The pertinent two-nucleon diagrams are shown in Fig. O, 

(2) . . . 

the blob characterizing an insertion from C\}^. In addition, there are relativis- 

tic corrections to the graphs in Fig. |4|. 

One can show that all terms proportional to the low-energy constants Ci 

vanish at threshold leaving insertions of terms ~ 1/2 A/, ~ gA/'^M and ~ ko,i. 

Counting powers of momenta, one sees that the contributions from some graphs 

at this order are divergent. This is due to the fact that in the chiral expansion 

one has truncated the pion-nucleon form factor and thus does not suppress the 
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Figure 43: Characteristic two-nucleon graphs contributing at 0{P'^) to neutral pion pho- 

(2) 

toproduction. The hatched circles denote an insertion from C^ff- 



high- momentum components.. This 
of the NN potential in xPT 
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hqnomenon also occurs in the calculation 
As was done there, we introduce an 
additional Gaussian cut-ofF factor of the form F{q'^) = exp{— g^/A^}, where 
the cut-ofF Ajr varies between the rho mass and A-^. Notice that the diver- 
gences appearing in other graphs all cancel and thus no further regularization 
is needed. 

Note that to 0(P'^) there are no four-nuclcon operators contributing to the 
deuteron electric dipole amplitude. The only possible graphs of this type that 
contribute at threshold to this order come from an NNNNtt vertex, to which 
a photon is attached cither at the outgoing pion leg or at the vertex itself. 
Both vanish for neutral pions as they result from minimal substitution in a 
pion derivative. We therefore conclude that for threshold neutral pion photo- 
production to 0{P^) there is no contribution from any four-nucleon operator 
and thus no new, undetermiTied parametei^appear. 
-Using the Argonne VlSlj, the ReidEIa, 
ll3 potentials, we find 
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the Nijmegenll£3 and the Paris 



E'/ = (0.36 ± 0.05) X 10"7m^+ 



(166) 



The P-wave contribution amounts to a 3% correction to the S wave. The 
sensitivity of the single-scattering contribution E^'^ to the elementary neutron 
amplitude is given by 



0.36 - 0.38 • (2.13 - E^^") x 10"7m^+ 



(167) 



Consequently, for — 0, we have E^'^ = —0.45 which is of opposite sign 

to the value in Eq. ( p^ based on the xPT prediction for E^^" (Table |). If 
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one were to use the empirical value for the proton amplitude (Tabic the 
single-scattering contribution would be somewhat reduced. 
We find for the two-nuclcon contribution at 0{P^) 

Ed'^ = (-1.90 , -1.88 , -1.85 , -1.88) x 10"^™^+ , (168) 



for the Argonne V18, Reid, super soft coreE£ll and Bonnlila potentials, in order. 
Further tiie result changes by Only a few percent if the chiral potential of Ref. 
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is used 
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From the graphs of 0{P ) we find 

Ed'^ = (-0.25 , -0.23 , -0.27) x 10-V"^7r+ , (169) 



for the Argonne V18, Reid and super soft core potentials, respectively, setting 
A^ = 1 GeV. These amount to corrections of the order of 15% to the 0{P^) two- 
nucleon terms. If one varies the cut-off A^r from 0.65 to 1.5 GcV, the 0{P^) two- 
nucleon contribution using the Argonne V18 potential varies between —0.24 
and —0.29 (in canonical units), which is a modest cut-off dependence. 

In summary, taking the Argonne V18 potential for definiteness, the chiral 
expansion of the electric dipole amplitude is shown in Table |^. Allowing 
E^j^ to vary between —1.0 and —1.5 and for a fixed value of SJ^" = 2.13, 
E^f varies between 0.2 and 0.4 (all numbers in canonical units). The results 
for the two-nucleon contributions are stable at 0{P^), whereas we assign a 
conservative uncertainty of ±0.1 to the corrections at next, order due to the 
cut-off dependence. We thus obtain the 0{P^) predictiovX^ 

Ed = {-l.8±Q.2)xlQ-^/m^+. (170) 



Table 5: Values for in units of 10 '^/m^+ from one-nucleon contributions (lA'^) up to 
0(P''), two-nucleon kernel (27V) at O(P^) and at 0{P^) , and their sum (lAf + 2N). 



IN 


2N 




lA^ + 2iV 


P + P^ +P^ + P^ 


pi 


pi p^ 


-P^ + P^ + P^ 


0.36 


-1.90 - 


-0.25 


-1.79 



To see the sensitivity to the elementary neutron amplitude, we set the 
latter to zero and find Ed = —2.6 x 10~^/to^+ (for the Argonne V18 poten- 
tial) which results in a cross-section twice as large as the one from the 
prediction (170). This is consistent with previous model resultsllI3ill3. 

This is a unique situation: the EFT makes a prediction that differs sig- 
nificantly from conventional models. Such a test was recently carried out at 
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Saskatoon E£3. The experimental results for the pion photoproduction cross- 
section near threshold are shown in Fig. together with the x^^^rW^dic- 
tion at threshold. Inelastic contributions have been estimated in Ref.E£3 and 
are smaller than 10% throughout the range of energies shown. At threshold, 



Ref.liiJ extracts 

Ed = (-1.45 ± 0.09) X 10" V™7r+ ■ (171) 

While agreement with xPT to order 0{P^) is not better than a reasonable 
estimate of higher-order terms, it is clearly superior to models. This is com- 
pelling evidence of the importance of chiral loops and of the consistency of 
nuclear EFT. A further test will come from an electroproduction experiment 
currently under analysis in Mainz. 
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Figure 44: Reduced cross-section an = {k/q)a in for neutral pion photoproduction as 
function of the photon-Miergy in MeV. Threshold is marked by a dptted line. Squares are 
data points from Ref. and the star is the xPT prediction of Ref. F^H Figure courtesy of 

Ulf Meifiner. 



8.6 -f*d^Tr°d 

In the case of an off-shell photon the threshold cross-section takes the form 

a{k^)^\Ed\^ + eL\LdW (172) 

with El = — (fc^/fcQ)e the longitudinal polarization for a photon with polariza- 
tion e. Now in addition to the transverse multipole E^, there is a longitudinal 
multipole Ld as well. Both of these multipoles are functions of the photon 
virtuality, k^. 
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These multipoles have been computed to 0{P^) in the chnal expansion. To 
this order there are graphs with one spectator nucleon and a contact operator, 
together with the two graphs of Fig. ^ with an ofF-shell photon. The proton 
electric dipole amphtude at the photon point is not well described to this 
order. However, the fc^-dependence of the elenientapz_amplitude for the proton 



is dominated by the third order contribution 
0{P^) results for Ed and Ld are shown in Fig. 




for a discussion). The 
In both cases, the two-body 
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Figure 45: S-wave multipoles of the deuteron. In the upper and lower panel, the electric 
dipole and the longitudinal multipoles are shown, respectively, by the solid lines. The dot- 
dashed line is the single scattering contribution. The dotted lines are obtained by changing 
the elementary neutron amplitude by ±1 X 10~^/M^+. 



contribution is large and dominates the amplitudes. Note that Ed varies more 
significantly with increasing |fc^| than Ld- As illustrated by the dotted lines in 
Fig. which have been obtained by a constant shift of the ir^n amplitudes 
by ±1 X 10~^/M^+, there is sensitivity to the elementary scattering off the 
neutron. To get a more realistic estimate for the S~wave cross-section, one 
can adjust the values of Xq_^_^ and Xq_^" {X = {E, L}) at the photon point to 
the values obtained from the O(P^ ) ca lculation. More precisely, the electric 



di pole amplitudes are taken from renlZj, the Lq^ from the best fit obtained in 



refliSa and the ij." using the resonance saturation estimate with = —125.6 
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and X' = —0.23, as detailed inE3. The so calculated S-wave cross-section is 
collected in Table |^ for a photon polarization e = 0.67. A range is given, which 
is obtained by adding a constant shift of ±1 x 10~'^/M^+ to the elementary 
n'^n longitudinal multipole (but keeping the electric dipole amplitude at the 
same value as before). This serves to illustrate the sensitivity of the S-wave 

comparison, note 



cross-section to the so far unmeasured Lq^ 



multipol^ 

that the S-wave cross-section measured on the proton 
-0.04 GeV^ lies between 0.15 and 0.45 ^b. For P a 



p For. 

184[l85 



for -0.10 <k^ < 



-0.1 GeV^, a ~ ao/2, 
with ao the S-wave cross-section for neutral pion production off the proton. 
However, the curvature of cr(fc^) is rather different from the proton case. This 
is partly due to the interference of the proton and neutron amplitudes and 
partly a kinematical effect, since for a given polarization e and virtuality fc^, 
cl is larger for the proton than for the deuteron. An O(P^) computation of 
this process is currently underway. 



-fc2 [GeV2] 


0.02 


0.04 


0.06 


0.08 


0.10 




0.11 


0.15 


0.17 


0.19 


0.20 


a range [^b] 


0.09 - 0.13 


0.12 - 0.18 


0.14 - 0.21 


0.16 - 0.23 


0.17 - 0.24 



Table 6: S— wave cross-section for the scaled single scattering amplitudes as explained in the 
text. The range is obtained by adding a constant shift of ±1 X 10~^/M^+ to the elementary 

TT^n longitudinal multipole. 



8.7 nd^ Tid 

To 0{P^) in xPT the vr-deuteron scattering length can be written asi 



(1 + m) 



(173) 



(I+m/2) 

where /i = m-^jM is the ratio of the pion and the nuclcon mass and the tt-N 
scattering lengths have the decomposition + Oirp — 2a+ — oS^°''> -I- a*^^^-* 
where is the isoscalar S-wave scattering length and the superscripts refer 
to the figures. The various diagrammatic contributions to aT^d are illustrated 
in Fig. 46 and Fig. |7[ The leading contribution. Fig. ^(a), has a vertex from 
and is therefore 0(P^). Fig. p§(b) has two nucleon propagators 



/•(2) /p2 

'-'ttN 

(P^^), one pion propagator (P^-^), three vertices from c'^^j^ (P^) and a two- 
body factor (P^) and so gives an 0{P^) contribution. One can further verify 
that the graphs of Fig. ^ are O(P^). Together, Figs. |6| and ^ are aU that 
contribute at 0(P^) at threshold. 
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Figure 46: Contributions to the -kcL scattering length at 0{P'^) (a) and 0{P^){b). The dots 



are vertices from C^^l, and the sliced blob is from C'" .,. 

71 N TTJV 



,(2) 



The contributions to a^rd from the graphs of Fig. 

, 1 > 



are: 



(2a) ^ 



>wf 



f2fc,2c) _ 



q-<Tiq- CT2 



Iwf 



(174) 



(175) 



where indicates that is sandwiched between deuteron wavefunctions. 

These matrix elements have been evaluated using a cornucopia of wavefunc- 
tions; results are displayed in Table 0. Clearly a*^^"^ dominates. This is the 
result of the shorter range nature of a^^^''^''\ It is important to stress that 
the dominant contribution from these graphs is quite independent of the wave- 
function used. This implies that xPT, which relies on the dominance of the 
pion-exchange, is useful in this context. 



-CO- 



Figure 47: Contributions to the 7r-deuteron scattering length at 0{P^). 
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for various deuteron wavefunctions in units of ^ . We 
:eV, gA = 1-32 and = 139.6 MeV. 



wf 1 








Bonn 

ANL-Vl 

ReidpSe 


17f 




-0.02021 
-0.01960 
-0.01941 
-0.01920 


-0.0005754 
-0.0007919 
-0.0008499 
-0.0006987 



Table 7: Contributions of Fig. ^ 
use U = 92.4 ^ 



To 0{P'^) in 



47r(l 



J IT 



2 

-4ci + 2c2 - — + 2c3 
4M ^ 



64^/4 ' 



(176) 



where the Ci are low-energy constants from The sole undetermined pa- 

rameter entering the 0{P^) computation of a-„d is therefore a combination of 
ci, C2 and C3: A = — 4ci -I- 2(c2 + C3). 



There is i;eceiit 



tering lengths 
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a 



gives 



aerimental information about both the tt N and 7rd scat- 
i23. Since involves constants that are not fixed by chiral 
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symmetry we can use experimental information about ttc? scattering to prec 
; the recent Neuchatel-PSTETHZ (NPE) pionic deuterium measurement 



-0.0259 ± 0.0011m: 



(177) 



For the contributions of Fig. ^ we take the average of the a*^^"-' and o^^"^) 
values in Table |: a'^'^''+^b+2c) ^ _o.0203m-i. We then find from Eq. ( ^ : 



-(2.6 ±0.5) • 10 



(178) 



Note that although a'^^'^'^-' is small, there is a strong cancellation between a^^"^ 
and flTrd which leads to a sensitivity f±a,a(^'''^'^'. Our value of a+ is not consistent 



with the Karlsruhe-Helsinki value 



19C 



= -(8.3 ± 3. 



the 



new NPE value deduced from the strong intoxariiinn shifts in pionic hydrogen 



and deuterium, which is small and positive 
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1+ = (0...5) • 10 



-3, 



The result Eq. (178)_agrees, however, with the value obtained in the SM95 



-3.0 • 10- 



partial-wave analysis 

Given the ambiguous experimental situation regarding a^, it seems most 
profitable to turn our formula around and use the ttc? scattering data to con- 
strain A. We can write 



(l + /i/2){a,d-(a(2'^)-|-a(2Mc))| 



^(1 
4M^ 



3Mm^s 
167r/2 ' 



(179) 
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Using Eqs. (|T7§) and (|T7^) we find A = -(0.08 ± 0.02) GeV"\ where 

we have taken into account the error in the determination of a^^- 

In Table || we give values of the relevant Ci obtained from a realistic_fi± 
to low-energy pion-nucleon scattering data and subthreshold parameters 
Central values lead to cr(0) = 47.6 MeV and a+ = -4.7-10" 
of the Ci give the conservative determination: A 
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These values 
-(0.18±0.75)GeV~^ Also 



Table 8: Values of the LECs Ci in GeV ^ for i = 1, . . . , 3. Also given are the central values 
(cv) and the ranges for the Ci from resonance exchange. The * denotes an input quantity. 



i 






^Kcs 


ranges 


1 


-0.93 ±0.10 


-0.9* 






2 


3.34 ± 0.20 


3.9 




2. ..4 


3 


-5.29 ±0.25 


-5.3 


-4.5. 


. - 5.3 


A 


-0.18 ±0.75 


0.8 


-3.0. 


. + 2.6 



shown in Table ^ are values of the deduced from resonance saturation. It is 
worth mentioning that an independent fit to pion-nucleon scattering including 
also low-energy constants related to dimension three operators finds results 
consistent with the fit values of Table 
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Nucleon Compton scattering has been studied in xPT in Ref. 
the following results for the polarizabilities were obtained at LO 



194,19. 



where 



384^2 j2^^ 

6 9 A 



12.2 X 10"'' fm^ 



^^=^"- 768^ - 1.2xl0-fm3. (180) 



Recent experimental values for the proton polarizabilities are ^ 



ap + Pp^ 13.23 ± 0.86l[! jJJ x 10"* fm^ 

ap-Pp^ 10.11 ± 1.74ljj2 X 10-4 fm^ (181) 

where the first error is a combined statistical and systematic error, and the 
second set of errors comes from the theoretical model employed. These values 
are in good agreement with the xPT predictions. 

"These are the result of a model-dependent fit to data from Compton scattering on the 
proton at several angles and at energies ranging from 33 to 309 MeV. 
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On the other hand, the neutron polarizabihties are difficuh to obtain ex- 
perimentaUy due to the absence of suitable neutron targets and so the cor- 
responding prediction is not well tested. One way to extract neutron 
polarizabihties is to consider Compton scattering on nuclear targets. Consider 
coherent photon scattering on the deuteron. The cross-section in the forward 
direction naively goes as: 



da 



ifTh-iap + an)L^'r. (182) 



1=0 



The sum ap + an may then be accessible via its interference with the domi- 
nant Thomson term for the proton, fxh- This means that with experimental 
knowledge of the proton polarizabihties it may be possible to extract those 
for the neutron. Coherent Compton scattering on a deuteron target has been 
measured at = 49 and 69 MeV by the Illinois group eZI and, usipg^ tagged 
photons, in the energy range Ej = 84.2 — 104.5 MeV at Saskatoonll£!l. Data 
for Ej of about 60 MeV is currently being analyzed at Lund. 

Clearly the amplitude for Compton scattering on the deuteron involves 
mechanisms other than Compton scattering on the individual constituent nu- 
cleons. Hence, extraction of nucleon polarizabihties requires a theoretical cal- 
culation of Compton scattering on the deuteron that is under control in the 
sense that it accounts for all mechanisms to a given order in xPT. There 
ns-pf this reaction in the framework of conventional po- 



exist a few calculatit 

tential models |123't2j'E^l!l. These calculations yield similar results if similar 
input is supplied, but typically mechanisms for nucleon polarizabihties and 
two-nucleon contributions are not treated consistently. We will see that 
provides a framework where this drawback is eliminated. 

The Compton amplitude we wish to evaluate is (in the jd center-of-mass 
frame) : 



'^'^ V^M' m T^^^N,,, ik',k) (P) (183) 



(27r)6 

where M {M') is the initial (final) deuteron spin state, and A (A') is the initial 
(final) photon polarization state, and k (fc') the initial (final) photon three- 
momentum, which are constrained to |fc| — \k'\ — oj. The amplitude Tj'* 



represents the graphs of Fig. 48 and Fig. |49|b where the photon interacts with 
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Figure 48: Graphs which contribute to Compton scattering on the deuteron at 0{P^) (a) 
and O(P^) (b-d). The sliced and diced blobs are from C^^^ (c) and (d). Crossed graphs 

are not shown. 



only one nucleon. Of course this amplitude must be evaluated in the jd center- 



of-mass frame. The amplitude T^^j^ represents the graphs of Fig. [49| b where 
there is an exchanged pion between the two nucleons. 

The LO contribution to Compton scattering on the deuteron is shown in 
Fig. |48|(a). This graph involves a vertex from and so is 0{P^). This contri- 
bution is simply the Thomson term for scattering on the proton. At next order, 
0{P^), there are several more graphs with a spectator nucleon Fig. ^(b,c,d) 
and as well as graphs involving an exchanged pion with leading order vertices 
(Fig. |49^) and one loop graphs, with a spectator nucleon (Fig. ^9|b). The full 
amplitudes are given in Ref.ll^. 




Figure 49: Graphs which contribute to Compton scattering on the deuteron at 0{P'^). 

Crossed graphs are not shown. 

For the wave function ip we use the energy-indepeadent Bonn OBEPQ 
wave function parameterization which is :^aufid in Ref. 
wavefunction computed to N^LO in Ref. 



138 



113 as well as the chiral 



The photon-deuteron T-matrix 
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in Eq. ( |183| ) is then calculated and the laboratory differential cross-section 
evaluated directly from it: 



Tli'VMxl'^ (184) 



M'X'MX 



where is the initial photon energy in the laboratory frame and E'^ is the 
final photon energy in the laboratory frame. Fig. ^ shows the results at 
Ej — 49, 69, and 95 MeV. For comparison we have included the calculation 
at 0{P^), where the second contribution in Eq. ( |183[ ) is zero, and the 77V 
T-matrix in the single-scattering contribution is given by the Thomson term 
on a single nucleon. It is remarkable that to 0{P^) no unknown counterterms 
appear. All contributions to the kernel are fixed in terms of known pion and 
nucleon parameters such as m^, qa, M, and f-^. Thus, to this order xPT 
makes predictions for Compton scattering. 

The curves show that the correction from the 0{P^) terms gets larger as ui 
is increased, as was to be expected. Indeed, while at lower energies corrections 
are relatively small, in the 95 MeV results the correction to the differential 
cross-section from the 0{P^) terms is of order 50%, although the contribution 
of these terms to the amplitude is of roughly the size one would expect from the 
power-counting: about 25%. Nevertheless, it is clear, even from these results, 
that this calculation must be performed to O(P^) before conclusions can be 
drawn about polarizabilities from data at photon energies of order mjr- This 
is in accord with similar convergence properties for the analogous calculation 



for threshold pion photoproduction on the deuteronE££ 




0.0 45.0 90.0 13S.0 laO.O ■ 0.0 4 5.0 90.0 135.0 130.0 



Figure 50: Results of the O(P^) (dotted line) and O(P^) (solid line) calculations for E-./ = 
49 MeV, 69 MeV and 95 MeV respectively from left to right. 

We have also shown the sb^llinois data points at 49 and 69 MeV0 and 



the Saskatoon data at 95 MeVE22l. Statistical and systematic errors have been 
added in quadrature. It is quite remarkable how well the 0{P^) calculation 
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reproduces the 49 MeV data. However, it is clear that the agreement at for- 
ward angles is somewhat fortuitous, as there are significant 0{P^) corrections. 
Meanwhile, the agreement of the 0{P^) calculation with the 69 MeV data is 
very good, although only limited conclusions can be drawn, given that there 
are only two data points, each with large error bars. 

These results are qualitatively not very different from other existing cal- 
culations. At 49 and 69 MeV the 0{P^) results are very close to those in 
■Rpf Fp and a few nb/sr higher, especially at back angles, than_thpse of Refs. 

At 95 MeV the 0{P^) result is_dose to that of Ref. H, higher by 



several nb/sr at back angles than Ref. E2ll^nd several nb/sr lower than the 



calculation with no polarizabilities of Ref. However the EFT calculation 

is the only one to incorporate the full single-nucleon amplitude instead of its 
polarizability approximation. The tendency to higher relative cross-sections in 
the backward directions is at least in part due to this feature. 

Although nominally the domain of validity of the Weinberg formulation 
extends well beyond the threshold for pion production, the power-counting 
fails at low energies well before the Thompson limit is reached. Consider the 
0{P'^) contribution shown in Fig. We can use this graph to illustrate the 
transition to the very- low energy regime ~ rn^/M. It is easy to see that this 
graph becomes comparable to the order P^ graph of Fig. |49| b when u>M. 
Here p is a typical nucleon momentum inside the deuteron and w is the photon 
energy. Since our power-counting is predicated on the assumption that all 
momenta are of order m^, we find that power-counting is valid in the region 
rn^/M ^ P <C . Therefore, in the region lo ^ B the Weinberg power- 
counting is not valid, since the external probe momentum flowing through the 
nucleon lines is of order P^/M, rather than order P. It is in this region that 
the Compton low-energy theorems are derived. Therefore our power-counting 
will not recover those low-energy theorems. Of course the upper bound on the 
validity of the effective theory should-jncrease if the A-resonance is included 



as a fundamental degree of freedom E2H. 

There are 0{P'^) counterterms which shift the polarizabilities. We use 
this freedom to fit the SAL data using the incomplete 0{P'^) calculation. A 
reasonable fit at backward angles can be achieved with a„ = 4.4 and /3„ = 10, 
which are in startling disagreement with the O(P^) X-P^^ expectations. That 
a fit can be achieved with these values of the polarizabilities is not surprising; 
continually increasing f3p + /3„ at approximately constant ap + a„ decreases 
the cross-section at forward angles and increases it at back angles. In fact, it 
seems that if (3p + /3n is sufficiently large then the character of the cross-section 



°At this energy Ref. E23 

only presents results with Oip + a-n = j3p + /3n = 0, wlpeh : 
are considerably less forward peaked than the corresponding calculation of Ref. t£j. 
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Figure 51: Interaction which contributes to Compton scattering on the deuteron at 0(P'*). 
The shced blob represents a 1/M correction vertex from C^^j- 

at 95 MeV can change completely from forward peaked to backward peaked. 
We have also plotted the cross-section at 69 MeV with a„ = 4.4 and /3„ — 10. 
This curve is then excessively forward peaked and misses the Illinois data. 
This situation poses an interesting theoretical puzzle. Of course, a full 0{P'^) 
calculation in xPT is necessary if the neutron polarizability is to be extracted 
from the Saskatoon data in a systematic way. 
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8.9 np d'y 

The chiral filter hypothesis is an idea based on current algebra 
hypothesis is that when a long-wavelength photon or axial probe interacts with 
the nucleus the dominant meson-exchange current contribution comes from the 
interaction of the probe with the soft pions. The jatguments given above for 
ed scattering suggest that, as first argued by RhoE23, this hypothesis can be 



organized into a systematic statement of the relative importance of various 
contributions to the NN system current, using Weinberg power counting. 

This can be clearly seen in the calculation of threshold np — > dj. In 
Ref.Sa Park et al. computed the Ml isovector matrix element for this reac- 
tion up to 0{P^) in Weinberg counting. At 0{P^) the LO contribution to the 
Ml transition enters. This is just the isovector counterpart of the one-body 
operator given above in Eq. (154): 



MlB 



— [(1 + Kp)ap + K„cr„] . 



(185) 



At O(P^) the two two-body currents depicted in Fig. ^ enter (although, of 
course, now a continuum np state is the final-state of the system). We will 
denote their contribution to the operator in question by /i2trec- Park et al. 
derive an analytic expression for this operator in co-ordinate space: 



M2tr 



9Aml ^ 

e—TTTT- (n X T2) < {ai X 0-2) 



-(1 + rn^r) - 1 
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f f • (cti X cTj) - ^ (o-i X (^2) 



[1 + m^r] 



e 



(186) 
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This exchange current was derived many years ago by Chemtob and Rho 
but it now has a definite place in the expansion of the two-nucleon current in 
powers of the small parameter P. 

The correction at next order, O(P^), vanishes, and so the next term in the 
series occurs at O(P^). At this order three types of term enter. 

1. Loop corrections to the ttNN and jttNN vertices. Thus, the only cor- 
rections comes from the counterterms which renormalize the divergences 
that appear in the loop graphs contributing to the jttNN vertex. Park 
et al. chose to fix these counterterms using resonance saturation. These 
low-energy constants could also be obtained by fitting to data in the 
single-nucleon sector. The topological structure of these graphs is the 
same as the left-hand graph in Fig. |5^,and these contributions are de- 
noted by /Z2vertex- 

2. There are also a number of one-loop, two-body exchange currents (not 
shown) denoted by /i2ioop- 

3. There is the right-hand graph in Fig. |5^. This is a contact interaction 
arising from the Lagrangian ( |7^ ) . Park et al. argue that this^ntribution 
is suppressed in their "hard-core cutoff scheme" (HCCS)E3a. However, 
if this process was reconsidered in the light of their more recent work, 
where a "modified hard-core cutoff scheme" (MHCCS) was employed, 
then the effect of this graph would be non-zero. 

Without this last contribution the magnetic moment operator is just the 
sum of all of these contributions, /i = /iib + /i2trGo + Ai2vortox + Ai2ioop each 
of which is calculable in terms of gA, /jr, ^jv, or low-energy constants from 
the single-nucleon sector. This operator must now be sandwiched between 
deuteron and np scattering state wave functions. For this purpose Park et al. 
employed the AV18 wave functions. The contribution from the one-body oper- 
ator gave a cross-section for thxeshold neutron capture, ctib = 305.6 mb, 
very close to the result from ERTl3. This simply means that the threshold 
capture cross-section is dominated by the physics of the NN wave function 
tail, which is the same in ERT and the AVIS potential. 

The evaluation of the two-body magnetic-moment operators is a little more 
subtle. We write the contribution of any given operator 11 to the Ml matrix 
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element as: 



(>X - j-T. \ r- (l»7j 

{d\niB\np) 



Park et al. cut-off the contribution from the operators fix at r = re ~ 1 /A, 
since the physics of this short-distance region wiU not be well-represented in 
xPT. Multiplying the operators written above by 9{r — Vc) and employing 
the AV18 wave functions they can then evaluate the ratios i5 as functions of 
Tc. The dependence on Tc is rather weak for a range of cutoffs from Vc = 
to Tc = 0.7 fm. Summarizing the results at Tc = 0.6 fmEJ J2trcc = 2.592%, 
<52vertex — 1-477%, and (52ioop = 0.362%. This leads to a theoretical prediction: 
a = 334 ± 3 mb, in excellent agreement with the experimental cross section 



Figure 52: Two two-body currents which contribute to E2s and Mlg at O(P^). 

As in the case of EFT(7f) this approach can also be used to evaluate the 
suppressed amplitudes M I5 and E2s^. The results obtained agree with those 
discussed previously, but differ somewhat in the error estimates. 



8.10 pp de^Ve 

Similar computations were performed in Ref.ll^j for the solar proton burning 
process pp — > de^Ve- The dominant piece comes from the one-body axial cur- 
rent. Park et al. show that this contribution can be evaluated with low-energy 
experimental data, p^, and r|sd^ as long as one is careful to include Coulomb cor- 
rections in the pp wave function. This gives a result: M°^g — 4.858 ±0.006 fm, 
where we have included the errors from the experimental data, added linearly. 
In contrast direct evaluation of the matrix element with the AV18 wave func- 
tion yields: M^^^^ = 4.859 fm. Neither of these numbers include a "theory- 
dependent" part of the matrix-element which comes from the interference be- 
tween the Coulomb and strong nuclear pp interactions. Using AV18 wave 
functions this additional piece of Mib is computed to be: M'^^^ = 0.011 fm. 
The contribution of vacuum polarization in the pp final state is also small, of 
order 0.5% of A^ib- Effects of two-photon exchange are smaller still. In fact. 
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all of these numbers are robust under the interchange of the AV18 wave func- 
tion with an EFT-motivated wave function, provided that in the EFT wave 
function the short-distance interaction is tuned to reproduce the scattering 
length and effective range of the pp system. 

The strength of these short-distance pieces is unconstr^ned by chiral sym- 
metry, although we expect it to be "natural" . As in Ref. tfl, in Ref. t££l it was 
argued that this counterterm contribution to pp de^v^ will be suppressed 
by the presence of the hard core in the NN interaction. It would be inter- 
estirig_,to see how the use of the MHCCS would modify the conclusions of 
Ref. 
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However, as in EFT(7f), this would require some way to determine the 
counterterm contribution. This would need to be done either by measuring 
a low-energy neutrino cross-section in the NN system, or by calculating the 
rate for the /3-decay of tritium, including the counterterm, and then fixing the 
counterterm to reproduce the extremely well-known experimental number, as 
discussed prejijiously 



In Ref. 



1^ only the axial currents of "pion range" were evaluated. These 
gave a contribution of between 2 and 5% to the overall matrix element, de- 
pending on what cutoff radius was employed. It would seem that this strong 
cutoff dependence might be due to the omission of the corresponding coun- 
terterm from the calculation. Regardless, the number resulting from this xPT 
calculation is still accurate at the lepiel of a few per cent. Park et al. quote 
results for the astrophysical S-factorL2l: 

S„{E) = a{E)Ee^'''\ (188) 

where a is the cross-section for pp — > de^Ve, and rj — Ma/{2p). The value 



S'pp(O) = (4.05 X 10"^^)(1 ± 0.0015 ± 0.0048)^ MeV barn, (189) 

where the first error is due to the uncertainty in the one-body matrix element 
and the second is Park et a/.'s estimate of the error on the two-body matrix 
element. 



In contrast the canonical value of Kamionkowski and BahcallE2a, calculated 
using a variety of NN potential models, is: 

S'pp(O) = (3.89 X ± 0.011) MeV barn. (190) 

The more sophisticated potential-model calculation of Ref. gives a much 
smaller error bar 

S'pp(O) = (3.88 X 10~^5)(1 ± 0.003) MeV barn, (191) 
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where the error comes both from a combination of smaU effects neglected in 
the calculation of the NN wave functions and from the observed sensitivity to 
the choice of modern NN potential. 



.11 Pion Production 



The pion production reactions: p + p p + p + n'^ , N + N — > d + tt, and 
N + N^p + n + TT would seem to be a natural candidate for calculation in 
xPT. However, the situation here is complicated by the appearance of a new 
scale. Since the energy of the incoming nucleons in any of theaa processes is of 
order to^ the typical nucleon momentum in the initial state isEj y^Mrn^. This 
renders the xPT expansion for these processes poorly convergent, since the 
expansion is now one in ^Jm^j^jM ^ 0.4. This makes it difficult to reproduce 
the expcrimcntaLdata, which is in fact-of very high quality, in both the case 
of neutral pion Ha and charged pion e2II production. There have been many 

to high-enough order that 
but, at this stage it seems 
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attempts to treat the S-wave p: 
some level of convergence is seen 
that the series is just not well-behaved. 

One solution to this is to resum certain mechanisms-Ln the xPT expansion 
in some non-systematic way, as discussed in Rcf. tIil'Ei3. Alternatively, one 
can turn to the production of P-wave pions. There one juight hope that the 
series is better behaved, since the LO contribution is largeEij. The results from 
early studies of the production of P-wave pions in xPT are encouraging, but the 
expansion parameter remains large. The advantage of employing xPT to study 
the production of P-wave pions is that, if the series is well-behaved enough, 
the experimental data might be able to be employed to extract information 
on contact operators involving four nucleons and a single pion. It is exactly 
this tyQ£_pf operator which enters in the novel three-nucleon forces discussed 
in Ref. EiH. The study of P-wave pion production therefore provides us with 
an opportunity to use xPT to systematically connect the three-nucleon system 
with pion-production in the two-nucleon system. 



8.12 Three-Body Forces 

In the Weinberg program the order in which the three-body forces start con- 



tributing depends on whether or not one is including_t 
us first discuss the case where the A's are present c^* 



le isobars explicitly. Let 
The three-body forces 



then first appear at order 0{P^) and are given by tree diagrams involving A's 
as intermediate states and zero, one and two pion exchanges. At order 0{P^) 
three classes of contributions appear. The first involves two-pion exchanges and 
a one-nucleon-two-pion vertex. The three independent one-nucleon-two-pion 
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vertices can be fixed by the study of pion-nucleon scattering. Tlie second class 
involves one pion exchange and a two-nucleon-one-pion vertex. The two inde- 
pendent coefficients of the two-nucleon-one-pion vertices can be determined by 
analyzing pion-deuteron scattering and/or pion production in nucleon-nucleon 
scattering. One particular linear combination of these two coefficients wasJ. 
termined by considering p— wave pion production iniV-t-A'^^iV-l-A'^-l-Tr 
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It is also known that if those coefficients have natural size they can appreciably 
affect the vector- analyzing power Ay in proton-deuteron scatt erin g that is not 
very well described by phenomenological two-body forces onlycH^. Finally the 
last category includes contact three-body pieces that can only be determined 
by fitting to three-body observables. In the version of the theory where the 
isobars are integrated out the three-body forces have no contribution at order 

and start at order P^. The contributions have the same form as the case 
with A's, but with different coefficients. 

Only recently the first calculations of three-body observables in the Wein- 
berg scheme were performedp^ calculation at the order P^ without explicit 
A's shows two main features c£2l. The binding energy of the triton is found to 
be close to the experimental value but shows a cutoff dependence that seems 
to be too large for a NLO calculation. The values obtained were — 8.284 
for A = 540 MeV and B3 = 7.546 for A = 600, but with a 10% variation in 
the binding energy for a 10% variation in the cutoff. It is presently unclear 
whether cutoff dependence is somewhat accidental or indicates that three-body 
forces that could cancel this dependence appear at lower order than predicted 
by Weinberg's power-counting, as shown to be the case in the pionless theory. 
The second feature is that spin-observables are well described by the NLO 
calculation, including the vector-analyzing power Ay . 

8.13 Formal Problems with Weinberg Power Counting 

As we have already discussed, at LO in Weinbergs power-counting there are 
(2) 

contributions to Vq from both the local four-nucleon operators, Cs,t and from 
the exchange of a single potential pion, giving a momentum space potential of 

\^fV (q2+m2) 

where C denotes the combination of Cs.t appropriate for a given spin-isospin 
channel. The LO amplitude results from solving the Schrodinger equation with 
this potential. This can be justified in the EFT only if M ^ P^^, otherwise di- 
agrams at different orders in Weinbergs expansion are numerically of the same 
size. At two-loops in the ladder sum there is a logarithmic divergence in the 
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(a) (b) 

Figure 53; Graphs with logarithmic divergences at LO in Weinberg's power counting. The 
solid lines are nucleons and the dashed lines are pions. 



graph shown in Fig. pi 3^ that must be regulated. In dimensional regularization 
the divergent part of this graph is 

which requires a counterterm with a single insertion of the light quark mass 
matrix, for instance 

_i4'")Tr[m,(E + Et)] [n^ -]^&^hr[m,{^ + {N^aNf + ... . 

(194) 

However, the coefficients of these operators must scale like M^, and since 



m 



M"^ ^ P'^ these formally higher order operators in Weinberg's power- 
counting are required at LO to absorb divergences in the time-ordered products 
of the LO potential, Vq' . Ignoring the multi-pion vertices arising from these 
operators, they can be re- absorbed into the LO operators with coefficients 
Cs,T- We are then in the situation where there is no chiral expansion, multiple 
insertions of the light quark mass matrix are not suppressed compared to LO 
interactions. 

The Schrodinger equation sums diagrams that are of different order in P. 
Paradoxically, it is consistent with Weinberg's power counting not to use the 
Schrodinger equation. In that case the near threshold bound states and the 
associated behavior of the scattering amplitudes are not recovered. It was 
shown in Ref.tS that taking M ~ P~^ corrects this problem. However, a more 
intricate scaling might in fact be in operation. 

Resolving these issues is currently the most important line of investigation 
in multi-nucleon EFT, since only then will there be an acceptable formal basis 
for the EFT at scales |k|> m^. 
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9 Gains over the Traditional Approach 



Traditionally, low-energy phenomena in light nuclei have been described using 
the Schrodinger equation with a model for the NN potential. In practice, all 
sensible potentials include the one-pion exchange potential (OPEP) as their 
long-range tail. The inclusion of OPEP is required on theoretical grounds, as 
it uniquely describes the longest-range part of the strong interaction. This 



s-essential in the description of 
lO. On the other hand, apart 



can also be seen experimentally, since OPE 
the higher partial waves in NN scattering 
from their common one-pion exchange tails, different NN potential models 
look rather dissimilar. Some include a set of single (or, in some cases, double) 
boson-exchangfiS-to generate the short-distance pip::^pQ£,the interaction, as in 
the NijmegenOilSy or the original Bonn potentials 
phenomenologir a 1 ..short-distance structure, e.g 



146,222 



Others ha 



Paris potential or the Argonne potentials&c22lQ, 



ilie,Reid potential 



ial^lil^. 



more 



the 

In all cases the short- 
distance potential must be rich enough to describe the existing data. Once 
a particular form is chosen the Schrodinger equation is solved and the free 
parameters are fit to the NN scattering data and certain deuteron bound- 
state properties. This is now an extremely sophisticated enterprise. Potentials 
exist which fit the np and pp scattering data up ta.laboratory energies of 350 
MeV with a per degree of freedom close to onefflQ. 

The treatment of short-distance physics employed in many of these poten- 
tial models can be motivated from an EFT point of view — Potentials ranging 
from the modern AVISH to the venerable Reid potential tI3 are based on the 
idea that the details of the short-distance physics in the NN problem are not 
important. It is important to introduce a structure for the short-distance op- 
erators that parameterizes our ignorance about that part of the interaction, 
but the detailed form of that parameterization is not crucial. This is per- 
haps not surprising: physically it is just not possible to resolve the details of 
short-distance physics using low-momentum observables. Therefore attempts 
to get the "right" short-distance NN potential using low-momentum data are 
futile. Potential models have predictive power precisely because the details of 
the short-distance physics do not matter greatly. 

It should therefore be clear that the goal of a xPT NN potential is not to 
learn about the short-distance structure of the NN interaction. The chiral La- 
grangian quite deliberately makes no dynamical assumptions about dynamics 



PThe CD-Bonn potential nominally involves single-meson exchanges at short dis- 

tances, but the mass of the scalar-isoscalar meson varies from partial wave to partial wave. 
Therefore in classifying OBE potentials as we have done here it is not entirely clear which 
class to place it in. 
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at scales of order A^, and so, by construction, the xPT NN potential contains 
only gross information about the physics which takes place at this scale. There 
is, however, the chance to learn something about the pion physics which is the 
provenance of xPT. 

Of course, "traditional" potential models are not just employed to fit NN 
scattering data. They also do a very good job in describing a wide array of 
low-energy data obtained using electroweak probes. To calculate the response 
of the nucleus to these probes, single-nucleon currents and potential-model 
wave functions are often employed. In some cases such a construction does 
not lead to a conserved electromagnetic current, and so additional "model- 
independentlLtwo-body currents are added to restore current conservation (see, 
for instancetSj). Conversely, transverse currents are not constrained by current 
conservation and so are always "model-dependent". The usual practice is to 
model these currents via something akin to the resonance saturation hypoth- 
esis used to estimate counterterms in the single-nucleon sector of xPT. The 
results for these transverse currents reflect dynamical assumptions about the 
important physical mechanisms, e.g. that pionic, p-meson, pn^f, luttj, Attj, 
and Ap7 meson-exchange currents dominate the transverse currents. 

Employing xPT to calculate these current operators leads to two gains 
over these more standard approaches. Firstly, as in the case of the NN poten- 
tial, it simplifies the picture of the dynamics down to its essential features: the 
"long-distance" part of the current operator comes from the physics of pions 
and nucleons interacting with the probe, while the short-distance part must be 
parameterized in some way. Secondly, and perhaps of more practical impor- 
tance, the calculation of current operators in xPT systematizes the treatment 
of the current operators in the NN system. In particular, xPT constrains the 
number and form of current operators allowable at a given order and mandates 
the type of pion-loop and pion-exchange mechanisms which should be included 
in the calculation. This means that, as in all EFT calculations, errors can be 
estimated in a reliable way. Such error estimation is possible in potential-model 
calculations (see, for instancelld), but it is significantly simpler in xPT. 

10 The Nuclear Many-Body Problem as an EFT 

In order for the EFTs that we have discussed in the previous pages of this 
chapter to have any impact whatsoever on nuclear structure, a way of solving 
the many-body problem must be developed, that is consistent with the power 
counting. As we do not know exactly what the final EFT will look like in 
the momentum region relevant to nuclei, we now turn our attention to the 
problem of solving the nuclear many-body problem for the hamiltonian given 
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in Eq. (|l|). The hope is that the EFT program wiU aUow enable us to match 
onto Eq. (|l|) and provide the tools to systematically calculate corrections to 
it. Equally important is that the many-body problem has to be solved anyway 
at some point! Although this problem has consumed nuclear physicists for 
several decades, the tools are still not available to make reliable calculations in 
bound systems with more than ~ 10 nucleons. Furthermore the techniques em- 
ployed in systems of 6-10 nucleons - Green's function Monte Carlo (GFMC) 
and variational methods - are sufficiently technical and challenging numeri- 
cally that very few theorists have access to them. In still heavier systems the 
techniques of choice, such as the shell model, are phenomenological, provid- 
ing no basis for systematic corrections or evaluating errors. For these reasons 
there is some interest in the possibility that the ideas of EFTs may point the 
way to new many-body techniques that are both simpler and more system- 
atic and controlled. Although the first such steps in this direction have been 
taken only recently-Jiere we summarized some recent work (that of Haxton 
and collabor|atop£23EH3) (for others investigations of multi-nucleon dynamics 
see e.g. Ref. cHZl'cHa) which suggests that conventional nuclear models do have 
controlled effective theory counterparts in which modern techniques, such as 
the renormalization group, can be employed to simplify the classical nuclear 
physics problem of point nucleons interacting through a static potential. A 
related and very interesting development that we will not discuss imzDhies sim- 



ilar progress in systematically treating low-density nuclear system£23'E23. The 
use of dimensional regularization in this scenario leads to great simplifications 
in various calculations. 

The canonical microscopic model for nuclear structure - the shell model 
(SM) - is motivated in many textbooks by Brueckner's treatment of nuclear 
matter. The language is at least superficially reminiscent of an EFT. The 
starting point is a collection of nonrelativistic nucleons interacting through a 
potential, e.g., of the Nijmegen or Argonne/Urbana type within an infinite 
Hilbert space 

1 ^ 

ij = l 

where Tij is the relative nonrelativistic kinetic energy operator and Vij the 
nucleon-nucleon potential. The related SM Hamiltonian acts in a restricted 
space and employs a softer "effective" potential. 



110 



Motivating Hs m is the notion that the determination of V^^^ might be simpler 
than solving the original A-body problem: the foundation of Brueckner theory 
is that high-momentum contributions to the wave function might be integrated 
out in a rapidly converging series in Pnuciear or, equivalently, in the number 
of nucleons in high-momentum states interacting at one time outside the SM 
space. This assumption - sometimes termed the hole-line expansion - clearly 
has some connections with our discussion of the power-counting for three-body 
forces, and thus may find some justification through EFTs. 

As will become apparent later, a typical SM space will contain explicitly 
~ 60% of the wave function that resides at long- wavelengths, thus guarantee- 
ing at least a qualitative description of A-body correlations important to soft 
collective modes. Implicitly the high-momentum components are swept into 
a rather poorly defined "effective interaction," often determined empirically. 
The strength of the SM resides in the first of these two aspects: the technology 
developed for direct diagonalizations in large SM sjpaces is quite remarkable, 
including recent progress inJjanczos-based methods E^, in treatmimts of light 
nuclei involving many shellsE^, and in Monte Carlo samplingE3E2j. Its weak- 
ness is the numerous uncontrolled approximations that become apparent when 
one tries to view the shell model as a faithful effective theory (ET). How- 
ever, we show below that the same numerical strides that have advanced shell 
model diagonalizations now allow us to remove these uncontrolled approxima- 
tions. The resulting ET has many differences with the shell model and many 
similarities to the EFTs discussed previously. 

Among the SM uncontrolled approximations are the following: 
1) Even in LO, where only the pairwise interaction of high- momentum nucleons 
is included in W^^^ , the functional form of the resulting effective interaction is 
not as simple as assumed in the SM, 

{\H^ff\)sM = {am^^^m (197) 

where the Greek symbols label single-particle shell-model states. As nucleons 
in high-momentum states can scatter pairwise repeatedly, this expression can 
at best be valid in some LO sense: beyond LO three-, four-, and higher-body 
body operators must be successively induced. Furthermore the energy denom 



inators describing this scattering depend on the total many-body energ]c2j, 
not simply the single-particle states participating in the scattering (apart from 
special cases, like a LO calculation restricted to a single shell). Thus the SM 
simplification of H'^^f is ad hoc. 

2) SM wave functions are orthogonal and normed to unity. In ET the effec- 
tive wave functions are naturally defined as the restrictions of the true wave 
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functions to the model space 



I*,) (5A/|*,)|5A/) = (198) 
ET 

Thus the norms are less than unity and orthogonality, which holds for the full 
wave functions, is lost for their restricted counterparts. 

3) Shell model interactions frequently depend on fictitious parameters such as 
"starting energies," introduced to adjust the energy denominators in the high- 
momentum ladder diagrams. Answers depend on these parameters, as well as 
on others (e.g., the oscillator strength if the SM is defined in terms of harmonic 
oscillator Slater determinants, and the number of included oscillator quanta) 
that define the division between the SM and excluded high-momentum spaces. 
Clearly the predictions of a correct ET must be independent of all such phys- 
ically irrelevant parameters. 

4) Perhaps most serious, the important issue of effective operators is almost 
never addressed in a meaningful way. In many cases practitioners adopted a 
phenomenological H"^^^ which, while successful in producing spectra, provides 
no diagrammatic basis for calculating effective operators or wave function nor- 
malizations. Even in cases where H"^^^ is derived from some underlying NN 
interaction, the practice is generally to then employ bare operators. In some 
well-studied cases, such as allowed /? decay in the Ip and 2s\d shells, it is 
then recognized that a phenomenological renormalization (e.g., qa 1) of 
operators greatly improves agreement with experiment. But the origin of this 
renormalization and its evolution with momentum transfer q are left unclear. 
The situation is very unsatisfactory and undercuts the shell model as a pre- 
dictive tool. 



10.1 Self- consistent Block- Horowitz Solutions 

We consider the canonical nuclear structure problem of nonrelativistic point 

3gen 
po- 



nucleons interacting through a realistic NN interaction, such as the Niun 
potential discussed previously, or the similar Argonne V18E| or ReidOS''^^' 



tentials. The question is whether the uncontrolled approximations in the shell 
model can be removed, leaving a more complicated but still tractable effective 
theory. The approach involves three major steps: 



• Formulating a treatment of effective interactions and operators that exploits 
the basic assumption in Brueckner theory — that interactions at high momenta 
can be integrated out in a cluster expansion (essentially an expansion in pa'^, 
where p is the nuclear density and a an interaction range) — but is otherwise 
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exact. The convergence of the expansion could then be tested numcricaUy and 
should depend on the operator under study and the momentum transfer. Al- 
ternatively (and preferably) one might be able to make connections to EFTs 
that would justify the expansion. 

• To find numerical tricks for implementing this formulation, demonstrating 
their validity in cases (e.g, A=2,3,4) where the expansion can be carried to all 
orders, so that the answers should then agree with Faddeev and other exact 
methods. 



• To imbed the formulation in a heavier nucleus, where the cluster expansion 
can be carried out only partially. 

There is some reason for optimism that if the first two goals can be 
achieved, the third might yield very accurate results: the Argonne group clus- 
ter variational Monte Carlo effort on ^^O appeared to yield nearly exact results 
when clusters up to A = 5 were included. Below we describe progress toward 
the first two goals, contrasting an exact FT of the deuteron and '^He with the 
shell model to illustrate the shortcomings of the later. 

The approach is sketched in Fig. The Hilbert space is divided into a 
long- wavelength "SM" space, defined by some energy scale Asm, and a high- 
momentum space that extends to some scale Aqo . Potential models effectively 
contain ultraviolet cutoffs embodied in their hard but finite short-range cores. 
Thus if Aoo is set to ~ 3 GeV/hui ^ 140 for the potentials mentioned above, 
ground state energies will be accurate to better than 1 keV. That is, for suffi- 
ciently large Aoo 

Hll/iAsM,Aoo) H'^iAsM). (199) 

Aoc> large 

All correlations within the "SM" space are included, but the high-momentum 
correlations in the excluded space are limited to n-body, where n is the cluster 
size. Thus the lowest order effective interaction is 

<4) - HliU.,) (200) 

It corresponds to embedding the A-body ladder diagram of Fig. |5^(b) between 
SM states: A-2 of the nucleons are spectators, with the remaining pair scat- 
tering via a two-body ladder. The notation (n=2,0) states that the two-body 
cluster has no explicit dependence on the nuclear density, varying as The 
n=3 A-body ladder of Fig. ^(c) is similarly 

Hff, = Hff , , + Hff (201) 

(n=3) {n=2,l) {n=3,0) ' 
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where H^^f{n = 2, 1) is the two-body part of the three-body ladder 

(«/3|^^(i£2,i)l"'/3'> = E ("/37|if(i£3)l"'/^'7 - «'7/?' + •••), (202) 

where kp denotes the Fermi level. This decomposition which is done only 
to emphasize the content of the cluster expansion illustrates that H^J^^^^ 

contains o) ^^^^ ^ correction to the two-body interaction that de- 

pends linearly on the density and is obtained by identifying one ingoing SM 
single-particle leg of the three-body ladder with one outgoing leg, summed over 
all occupied states. It also contains a true three-body piece H'^^L^ q) , where 
three SM ingoing single-particle states connect to three distinct outgoing states 
after undergoing a series of scatterings outside the SM space. The point is sim- 
ple pedagogy: treatments of successively larger clusters in the high-momentum 
space correct the lowest-order two-body H^^^2 o) adding terms proportional 
to the p, p^, etc., in the spirit of Brueckner theory. It also adds true three- 
body terms, true four-body terms, etc. Thus an expansion through four-body 
clusters yields H'^nL2 2)' two-body interaction corrected through order p^, 
-^(^=3 1) ' three-body interaction through order p, and H^^^^ , a density- 
independent true four-body interaction. We emphasize that the convergence 
of such an expansion is still a matter of conjecture, despite its physical plau- 
sibility and some justification from EFT power counting. The procedure has 
not yet been given a chance to succeed because the individual steps - certainly 
those beyond the two-body level - have not been carried out. We now discuss 
procedures that might make such an expansion feasible. 

The calculation begins with a definition of the "SM" space. The goals of 
handling bound states and of generating an effective interaction that is trans- 
lationally invariant leaves one sensible choice, many-body states constructed 
from harmonic oscillator Slater determinants. To exploit the relative/ccnter-of- 
mass separability of harmonic oscillator Slater determinants, one must separate 
the SM and high-momentum spaces so that all configurations satisfying 

E < AsMfico (203) 

are retained in the former. For example, a SM calculation of with Asm ~ 
4-1- Ag, where Aq is the number of quanta in the ^^O closed shell, would include 
all 4?iti> configurations, e.g., OpO/i, 2p2h. and ApAh excitations of nucleons from 
the Ip shell into the 2sld shell, Iplh excitations of a Is shell nucleon into 
the is2dlg shell, etc. One can define the projection operator onto the high- 
momentum space by 

QsM = Q{AsM,h)- (204) 
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Figure 54: Cluster expansion of the effective interaction. 



where b is the oscillator parameter. Thus the included or "SM" space is defined 
by two parameters, Asm and b. The preservation of translational invariance is 
also important numerically, as it reduces the two-body ladder to an effective 
one-body problem, etc. 

is then 



The resulting Bloch-Horowitz (BH) equation 

1 



236 



H-ff 



H + H- 



tQsmH 



E - QsmH " 
= S|*sm) |*sa/) = (1 - QSA/)!*) 



(205) 



where |^) is the exact wave function and H\'i>) = The difficulty posed 

by this equation is the appearance of the unknown energy eigenvalue in the 
equation for H'^'^f . Thus this system must be solved self-consistently. Note 
that there is no explicit reference to the harmonic oscillator in this equation: 
it enters only implicitly through Qsm in distinguishing the long- wavelength 
"SM" space from the remainder of the Hilbert space. We emphasize that a 
proper solution of the BH equation must yield results (energies and operator 
matrix elements) that are independent of Qsm- 

There is an extensive literature on this and similar equations, often involv- 



115 



iog-cLxlivision of H into an unperturbed i/o and a perturbation Hi = H — Hq 
|237|j238| xhere are well-known pathologies with this division ipznlying the ef- 
fects of near-by intruder states on the perturbation expansionc^jElJ. Here we 
explore another approach that is nonperturbative and involves, in effect, a com- 
puter summation of diagrams. The method is based on the Lanczos algorithm 
and offers a remarkably simple solution to the issue of self-consistency. 

In the Lanczos algorithm a basis for representing a Hamiltonian is formed 
recursively in such a way that the resulting Hamiltonian is tridiagonal. Given a 
Hermitian operator H and and an initial normalized vector \vi), the successive 
steps are 



H\vi) = ai\vi) + (3i\v2) 

H\V2) = f3l\vi) + a2\v2) + I32\V3) 

H\v3) = P2\v2) + aalva) + I33\v4) etc. 



(206) 



so that the H takes the form 



H 



( ai f3i 

Pi 




V 







^2 f32 
P2 "3 













1^3) 


/ 





(207) 



The remarkable property of this algorithm has to do with truncating the pro- 
cess in Eq. (206) after n steps, where n can be much smaller than the dimension 
of the Hilbert space. The resulting truncated matrix in Eq. (207) then con- 
tains the information needed to reconstruct the exact 2n-l lowest moments 
of H over the eigenspectrum. As extremum eigenvalues are crucial to higher 
moments, one common application of the Lanczos algorithm is in determining 
such eigenvalues and their associated eigenfunctions. Another is to begin with 
the vector \vi) — 0\g.s.) and then use the algorithm to calculate the moments 
of the response of the ground state \g.s.) to the operator O. A small number of 
moments, e.g., ~ 100, often is sufficient to construct a response function with 
a numerical resolution com parable to that achieved experimentally. 

A third application 
One finds 

1 



is in constructing fully interacting Green's functions. 



E-H 



\vi) = gi{E)\vi) + g2{E)\v2) + 



(208) 



where the gi{E) are continued fractions that depend on ai,f}i and where E 
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appears only as a parameter. For example, 



9i{E) = 



1 



(209) 



E — a\ 



E-a2- 




It follows that the BH equation can be solved self-consistently with only a 
single solution of the effective interactions problem, even in cases where mul- 
tiple bound states are needed. The procedure is: 

• For each relative-coordinate vector in the SM space I7), form the excluded- 
space vector \vi) = QsmH\')) and the corresponding Lanczos matrix for the 

operator QsmH. Retaining the resulting coefficients on, Pi for later use, con- 
struct the Green's function for some initial guess for E and then the dot product 
with {i\H to find {i\H<'ff {E)\-f). 

• Perform the "SM" calculation to find the desired eigenvalue E' which, in 
general, will be different from the guess E. Using the stored ai,pi, recalculate 
the Green's function for E' and H'^f ^ {E') then redo the "SM" calculation. Re- 
peat until convergence, i.e., until the input E' in the Green's function equals 
the desired output "SM" eigenvalue. 

• Then proceed to the next desired bound state, e.g., the first excited state, 
and repeat the above step. Note that it is not necessary to repeat the H^^^ 
calculation. The eigenvalue taken from the "SM" calculation is, of course, that 
of the first excited state. The procedure then generates distinct H*^^^ {E')s for 
each desired state. 

The attractiveness of this approach is that the effective interactions part of 
the procedure, which is relatively time consuming as it requires one to perform 
a large-basis Lanczos calculation for each relative-coordinate starting vector 
in the "SM" space, is performed only once. The diagonalization in the model 
space is generally much faster: modern workstations can handle even large- 
dimension shell model calculations (sparse matrices of ^ 10^) quickly (^ 30 
minutes). In practice we found that self-consistency is achieved easily: six to 
eight cycles is typical. (More cycles are required for states with small binding 
energies.) Thus it is quite practical to derive the exact H'^^^{E)s for a series 
of bound states. 

In our view the Lanczos approach to the effective interactions problem 
appears to be remarkably simpler than the standard procedures of the field, 
particular in view of the need for extensions to multi-nucleon ladders. The 
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traditional approach divides the effective interactions problem into an energy- 
dependent piece (often called the Q-box) and an energy-dependent one (rep- 
resented by folded diagrams). Frequently the energy dependence is removed 
by making a unitary transformation to a nonHermitian effective Hamiltonian. 
The Lanczos procedure is also much more in the spirit of EFT: the full Hamil- 
tonian in the high-momentum space is never constructed. Rather, it is replaced 
by a much smaller Lanczos matrix which contains exactly the most relevant 
long-wavelength information of the full matrix, the 2n-l lowest moments. Thus 
the procedure can be viewed as a numerical ET in which this information is 
recursively extracted. 

Now we discuss the results of applying this procedure to the simplest nuclei, 
d and ^He, carrying the above process to completion (two-body and three-body 
ladders, respectively). The motivation is two-fold: demonstrate the numerical 
procedures we described above, and provide for the first time exact effective 
theory results that can be compared to those of the shell model. 

The harmonic oscillator mode expansion must be sufficient to represent 
both the long-distance tails of bound states and the short-distance "hard core" 
scattering predicted by realistic NN potentials. (The Argonne V18 po- 
tentials are shown in Fig. |5^.) Inclusion of high- momentum states through 
Aoo ~ 50 yields a deuteron binding energy accurate to ~ 60keV; extending 
this to ^ 140 produces a result accurate to one keV. Fig. |5^ shows the rate 
of convergence as a function of b and Aoo. The convergence for the binding of 
■^He (not shown) is slightly more rapid, a feature related to the deeper binding, 
as will become clear later. 




Figure 55; The Argonne V18 ^Sq potentials. 
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The binding energies and operator matrix elements for simple systems like 
■^He can, of course, be calculated by other methods, e.g., Faddeev techniques 
or GFMC. We thus want to stress that the point of the following discussion is 
to do analogous calculations in the context of an effective theory, so that we 
begin to see the shortcomings of conventional techniques like the shell model 
as well as possibilities for overcoming those shortcomings. A first test of the 
techniques outlined above is to solve the BH equation for some SM-like space 
to then see if the resulting self-consistent energy is, indeed, the correct value. 
For model spaces of 2, 4, 6, and 8/iw in the case of the deuteron we obtained a 
binding energy of -2.224 MeV (using y^fe = 1.6/ and Aoo = 140). The exact 
result is -2.2246 MeV. Our ^He result for same SM spaces and Aqo — 60 is 
-6.87 MeV, in agreement with the corresponding GFMC result of -6.87 ± 0.03 
MeV. 




Figure 56: Deuteron ground state energy convergence as a function of Aoo for several choices 
of the oscillator parameter b. This b is defined as in the independent-particle SM: the 
corresponding b for relative motion is \/2 times the values shown. It is clear that if an 
inappropriate size is chosen for the basis states (e.g., b=1.9), the rate of convergence can be 

greatly slowed. 

More interesting is the evolution of the wave functions, shown in Tables || 
and |l^. The tables were generated as above: for some large value of A^a 
SM wave functions were generated for various choices of Asm by interating 
the SM calculation until the self-consistent energy is fully converged. The 
deuteron and "^He calculations involve two- and three-body ladder sums in the 
excluded space, yielding sets of two- and three-body "SM" matrix elements 
of H'^^f for the model spaces. The deuteron calculation is rather trivial; for 
Aoo ~ 50 the ■^He BH calculation involves a dense matrix of dimension ~ 10'', 
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Table 9: ET ra oulito f pr the deuteron ground state wave function calculated with the Argonne 
V18 potential r^^rr^^ . The columns on the right correspond to different choices of the ET 
model space, the analog of a SM space. The rows correspond to the resulting amplitudes for 
the designated, selected configurations | n, ^ Lj). The quantities within the parentheses are 
the square of the norm of the effective wave function, e.g., the probability that the deuteron 
resides in the corresponding "SM" space. 



basis state 


amplitude 




2nuj 


Ahuj 




8 hu 


exact 


(65.9%) 


(79.5%) 


(86.1%) 


(91.3%) 


(93.0%) 


(100%) 


1 1, ^^i> 


0.81155 


0.81154 


0.81155 


0.81155 


0.81152 


0.81155 


1 2, ■''S^) 


0.00000 


-0.31483 


-0.31483 


-0.31483 


-0.31482 


-0.31483 


1 1, 'D,) 


0.00000 


0.19524 


0.19524 


0.19524 


0.19523 


0.19524 


1 3, 


0.00000 


0.00000 


0.24945 


0.24945 


0.24944 


0.24945 


1 4, ■''S^) 


0.00000 


0.00000 


0.00000 


-0.20851 


-0.20850 


-0.20851 


1 5, 


0.00000 


0.00000 


0.00000 


0.00000 


0.12596 


0.12596 



still rather modest by current SM standards. The matrix is dense because 
we work in relative Jacobi coordinates, r athe r than anjibscheme, utilizing 



standard Talmi-Brody-Moshinshy methods E13. (See Ref. El^ for details.) The 
results in Table ^ were obtained with approximately 100 Lanczos iterations: 
it is apparent that the convergence is then quite good. The wave functions 
must be normalized according to Eq. (205): this involves calculating unity as 
an effective operator. We will return to this point below. 

The tables show the lovely evolution of the wave function in ET, an evolu- 
tion quite unlike that of typical shell model calculations. The wave functions 
obtained in different model spaces agree over overlapping parts of their Hilbert 
spaces. Thus as one proceeds through 2huj, Ahuj, 6huj, ... calculations, the ET 
wave function evolves only by adding new components in the expanded space. 
The normalization of the wave function grows accordingly. Thus, for '^He, the 
Ohuj ET calculation contains 0.311 of the full wave function in the effective 
space; the 0+2+4?kj result is 0.700. 

This evolution will not arise in the standard SM because the wave function 
normalization is set to unity regardless of the model space. It will also not arise 
for a second reason, illustrated in Table The matrix elements of H'^^^ are 
crucially dependent on the model space: the listed results for •^He show that 
a typical matrix {a\H'^-^^l3) changes very rapidly under modest expansions 
of the model space, e.g., from 2?ill} to Ahuj. Yet it is common practice in 
the shell model to expand calculations by simply adding to an existing SM 
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Hamiltonian new interactions that will mix in additional shells. We suspect 
the behavior found for "^He is generic in ET calculations: it arises because a 
substantial fraction of the wave function lies near but outside the model space 
(e.g., see Table ^0|). An expansion of the model space changes the energy 
denominators for coupling to some of these configurations, and moves other 
nearby configurations from the excluded space to the model space. Naively, 
relative changes in effective interaction matrix elements of unity are expected. 



Now we turn to the question of operators. The standard procedure in the 
SM is to calculate nuclear form factors with bare operators, or perhaps with 
bare operators renormalized according to effective charges determined phe- 
nomenologically at = 0, using SM wave functions normed to 1. As we now 
have a series of exact effective interactions corresponding to different model 
spaces, we can test the validity of this approach. The results for the elastic 
magnetic form factors for the deutcron and "^He are shown in Fig. Is^. One sees 




3 456 01234 

q(f') q(f') 



Figure 57: The magnetic elastic form factor for the deuteron (left panel) and ^He (right 
panel) calculated with the exact H"-^-^ , SM wave fun qliu i i i ih i i ormalized to unity, and a bare 
operator are compared to the exact result (solid line) P^5H226| effective operators and 



the proper wave function normalizations are used, all results become identical to the solid 

line. 



in each case that by the time one reaches a momentum transfer q ^ 2.5/ f, 
random numbers are being generated: bare operators used in conjunction with 
exact effective wave functions generate results that differ by an order of mag- 
nitude, depending on the choice of the model space. This is not surprising, of 
course. If one considers the operation 

e^«^-VT3|5.s.) (210) 

at momentum transfers ^ ikp, where kp is the Fermi momentum, most of the 



121 



Table 10: As in Table |9| only for ^He. The basis states are now designated somewhat 
schematically as | N,a), where N is the total number of oscillator quanta and a is an index 
representing all other quantum numbers. 



state 


amplitude 


OHlj 


2hLj 




6huj 


8 huj 


exact 


(31.1%) 


(57.4%) 


(70.0%) 


(79.8%) 


(85.5%) 


(100%) 




0,1) 


0.55791 


0.55791 


0.55791 


0.55795 


0.55791 


0.55793 




2,1) 


0.00000 


0.04631 


0.04613 


0.04618 


0.04622 


0.04631 




2,2) 


0.00000 


-0.48255 


-0.48237 


-0.48243 


-0.48243 


-0.48257 




2,3) 


n nnnnn 

u.uuuuu 


n nn79Q 

u.uu t 


u.uu t Ol 


u.uu ( ou 


n nn79Q 

u.uu I 


u.uu / 




2,4) 


0.00000 


0.16707 


0.16698 


0.16706 


0.16706 


0.16708 




2,5) 


0.00000 


0.00566 


0.00564 


0.00565 


0.00565 


0.00566 




2,6) 


0.00000 


-0.00017 


-0.00017 


-0.00017 


-0.00017 


-0.00017 




4,1) 


0.00000 


0.00000 


-0.02040 


-0.02042 


-0.02043 


-0.02047 




4,2) 


0.00000 


0.00000 


0.11267 


0.11274 


0.11275 


0.11289 




4,3) 


0.00000 


0.00000 


-0.04191 


-0.04199 


-0.04208 


-0.04228 




4,4) 


0.00000 


0.00000 


0.28967 


0.28978 


0.28978 


0.29001 




4,5) 


0.00000 


0.00000 


0.01059 


0.01059 


0.01059 


0.01059 




4,6) 


0.00000 


0.00000 


-0.00213 


-0.00212 


-0.00211 


-0.00210 




4,7) 


0.00000 


0.00000 


0.00998 


0.01000 


0.01000 


0.01000 




4,8) 


0.00000 


0.00000 


-0.11319 


-0.11327 


-0.11330 


-0.11335 




4,9) 


0.00000 


0.00000 


0.08446 


0.08447 


0.08446 


0.08448 


1 4, 10) 


0.00000 


0.00000 


-0.08613 


-0.08626 


-0.08632 


-0.08638 


|4,11) 


0.00000 


0.00000 


-0.00210 


-0.00211 


-0.00211 


-0.00211 


1 4,12) 


0.00000 


0.00000 


-0.00252 


-0.00254 


-0.00256 


-0.00257 


1 4,13) 


0.00000 


0.00000 


-0.00020 


-0.00020 


-0.00020 


-0.00020 


1 4, 14) 


0.00000 


0.00000 


-0.00010 


-0.00010 


-0.00010 


-0.00010 


1 4,15) 


0.00000 


0.00000 


-0.00012 


-0.00013 


-0.00012 


-0.00012 



Table 11: Selected BH 3-body effective interaction matrix elementsc23EHj for ^He, in MeV, 
illustrating the strong dependence on the "SM" space. The Argonne V18 potential was used. 





2huj 




6hLU 


8 huj 


(0,1 1 H-ff 1 2,1) 


-4.874 


-3.165 


-0.449 


1.279 


(0,1 1 H-f^ 1 2,5) 


-0.897 


-1.590 


-1.893 


-2.208 


(2,1 1 H-f^ 1 2,2) 


6.548 


-2.534 


-4.144 


-5.060 
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resulting amplitude should reside outside the long-wavelength model space, in 
any simple view of the nucleus. That is, the strength resides entirely in the 
effective contributions to the operator. If these components are ignored, the 
results have to be in error. 

Clearly the effective interaction and effective operator have to be treated 
consistently and on the same footing. If O is the bare operator, one finds 

(*/|6|*,) = (*}^-^|0^^^|*f^) (211) 

where 

O^ff = (1 + HQsMj^ ^77^)^(1 + IF-^ ttQsmH) (212) 

and where the effective wave function normalization of I^E'^''^''^) and \^'^^^^), 
mentioned earlier, must be determined using the effective operator 1, e.g., 

1 = = {^ff\{l + HQsM- ^7^)(1+ ^ \ „ QsMHmff) 

(213) 

These expressions can be evaluated with the Lanczos Green's function methods 
described earlier. When this is done, all of the effective calculations, regardless 
of the choice of the model space, yield the same result, given by the solid lines 
in Fig. 1^. 

We would argue, based on this example, that many persistent problems in 
nuclear physics — ranging from the renormalization of in (3 decay to the 
systematic differences between measured and calculated Ml electromagnetic 
form factors — very likely are due to naive treatments of operators, treat- 
ments that fail to satisfy the basic rules of ETs. It should be apparent from 
the above example that no amount of work on H'^^f will help with this prob- 
lem. What is necessary is a diagrammatic basis for generating H'^^f that can 
be applied in exactly the same way to evolving O'^^^. From this perspective, 
phenomenological derivations of H^^^ by fitting binding energies and other 
static properties of nuclei are not terribly helpful, unless one intends to simul- 
taneously find phenomenological renormalizations for each desired operator in 
each range of interest. 

The point of this discussion has been largely pedagogical. To our knowl- 
edge the examples given above are the only ones in the classical nuclear physics 
literature in which a "model space" calculation has been formulated in a way 
that satisfies the basic rules of ET. The results — particularly the invariance 
of energies and operator matrix elements under changes in how we define Q 
— are obvious from the perspective of ET. The fact that standard techniques 
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like the SM violate so many of the rules of ETs is cause for optimism: much 
can be done to improve the rigor of such nuclear physics tools. As we have 
demonstrated here, new numerical methods can be developed to handle the 
full ET problem at the cost of only modest additional effort. 



10.2 Crossing the Many-Body Border 

The effort to derive SM effective interactions directly from realistic NN interac- 
tions was largely abandoned in the early 1970's when it became apparent that 
the high-momentum scattering problem was nonperturbive. A series of efforts 
involving perturbative expansions in either the potential V or the correspond- 
ing G-matrix (the sum of the two-body ladder diagrams) yielded discouiaging 
results, with higher-order terms often dominating over lower-order one£lfl. In 
particular Shucan and WeidenmuUer showed that any overlap in the spectra of 
indn ded ari d excluded, high-momentum states would lead to a nonperturbative 
serie£2SElj. Thus, the uncontrolled approximations in the SM approach grew 
out of necessity: it became a largely phenomenological tool in part because 
early efforts to provide a more rigorous basis for the model ended in failure. 

Given our two SM-like ET calculations, the deuteron and ^He, it is in- 
teresting to return to this old problem of the nonperturbative behavior of 
the nonrelativistic many-body problem, using these examples as laboratories. 
Today's challenge, in fact, is more modest than that of the 1970's: modern 
computers coupled with new algorithms, like the Lanczos solution of the BH 
equation introduced above, allow us to solve rather complex nonperturbative 
problems. To extend the calculations above to ^He and heavier nuclei one has 
to achieve only modest increases in efficiency, such as a substantial loweidng of 
Aoo- Below we summarize some very recent work by Haxton and LuiElj that 
suggests that the goals of the 1970's — a perturbative treatment of the nonrel- 
ativistic nuclear physics problem — may not be out of the question. Through 
renormalization group and other techniques, the most severely nonperturba- 
tive contributions to the potential can be systematically extracted and treated. 
The approach is quite similar to that of EFTs. We conclude our discussion of 
the many-body problem by sketching this effort to cross the border between 
EFTs and conventional nuclear models. 

The deuteron calculations were performed for Aqo ~ 140, a value where the 
hard core of the AV18 potential is resolved sufficiently well that the resulting 
binding energy is good to 1 keV. Results are variational in this parameter: as 
Aoo is increased, the correct binding energy is approached from above. The 
evolution was explored in Ref.l22a and found to vary as 

Egs 2.175MeV -I- 3.178MeVe~^ °^^(^~/^°)' (214) 
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to very good accuracy, a functional dependence that already suggests some 
nonperturbative behavior. 

The BH calculation determines the full set of matrix elements {a\H'^^f\l3) 
for the SM space. These individual matrix elements carry information about 
the convergence not easily seen in Egs . Their behavior in perturbation theory 
can be examined by expanding the BH propagator 

1 1 1 . X 1 . X 

-Q{V-V^)— — + (215) 



E — C^E[ E — Hq E — E — 

evaluating the series term by term. Here and Vo are the harmonic oscillator 
Hamiltonian and potential. We thus lower Aqo and test whether states above 
Aoo can be treated perturbatively. 

Once a scale of K^o ^ 70 is reached, striking differences appear in the 
rate of convergence of different matrix elements. All are greatly improved in 
first and second order perturbation theory. Those matrix elements well away 
from the "boundary" at Aoo converge quickly to their correct values in this 
way. However, those matrix elements {a\H'^^f\l3) where either \a) or \(3) re- 
sides in the last shell (with TV = Aoo) converge only very slowly to the correct 
values, after the initial improvement in low-order perturbation theory. Typ- 
ical ~ 10'^ orders of perturbation theory are required to produce the correct 
value. This clearly suggests that slow convergence is associated with the rela- 
tive kinetic energy operator QT contribution to QH ^ as the only transitions to 
states outside the Hilbert space generated by this operator have Art = 1. At 
large r the strength of this transition becomes quite large, —QVo{r)\a) , reflect- 
ing the unphysical asymptotic behavior of harmonic oscillator wave functions. 
This amplitude propagates nonperturbatively as {V — Vo)/Hq ^ 1. Eventu- 
ally enough high-momentum harmonic oscillator wave functions are coupled 
together to produce the softer asymptotic fall-off characteristic of the correct 
bound-state wave function- 



As discussed in Ref. the required resummation is guided by the ob- 
servation that the true potential, V{r), falls off properly at large r. Thus a 
reorganization of the BH equation in which the propagator is always sand- 
wiched between V{r) should remove the unwanted propagation. This leads to 
the following recasting of the BH equation: 

{a\H^ff\(3) - {a\m + i{a\~{a\)E~Tm^\f3)) 

+ i&\V + V ^^^^ QV0) (216) 

where 

|a> = ^\a) (217) 
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If a and (3 are not in the last shell, |d) — \a) and \(3) — so that the 
above rewriting of the BH equation reduce s to the original form. Otherwise, a 
modified wave function generated by Eq. ( |217] ) must be used. 

This discussion is a necessary diversion from the principal topic of this 
section: this source of nonperturbative behavior dominates one that is of more 
interest (and more difficult to remove), hard core scattering. It should be 



clear the the lowering of the scale Aqo at the cost of an evaluation of Eq. (217) 
represents a tremendous savings: in place of a dense matrix QH whose elements 
have to evaluated numerically, we have a sparse matrix QT whose matrix 
elements are known analytically. The kinetic energy Green's function is thus 



trivial to evaluate via Eq. (208). For example, in the case of the deuteron, QT 



is tridiagonal in the harmonic oscillator basis. Thus the Green's function can 



be expressed analytically in the form of Eq. (208). 

Once the replacement \a) — > \a) is made, can be lowered to ~ 40 while 
maintaining 1 keV accuracy, if one works to third order in perturbation theory. 
But with further lowering, errors arise in the perturbative expansion that per- 
sists even if calculations are carried to very high order. At Aqo — 30, third order 
perturbation theory reduces 10% errors in (n' = 11' = Q\H'^^^n = 1/ = 0) 
to ~ 0.2%. But an error in excess of 0.1% — corresponding to 50 keV — 
persists after 10 additional orders of perturbation. Numerically one can verify 
that the nonperturbative tail is generated by the scattering at very small r. 
As expected, the nonperturbative contributions to s — d matrix elements are 
much smaller than those for (n' = W = Q\H'^^f\n = 1/ = 0) and other s — s 
matrix elements. 



We sketch an elegant solution to this problem described in Ref.Elj. The 
most singular short-ranged contributions are removed from the potential ana- 
lytically through a renormalization group (RG) procedure. The method starts 
with the well-known Talmi integral expansion of harmonic oscillator radial ma- 
trix elements and involves a systematic removal of the most singular of these 
by introducing local, separable operators familiar from EFT treatments. The 
coefficients of these operators are then run by the RG equations. The result 
establishes a very close connection between potential descriptions and EFT 
treatments, and also allows a clean separation of long-distance contributions 
to the potential that can be treated perturbatively. 

As our goal here is to illustrate the procedure, we work to LO. (The NLO 
and N^LO results are given in Ref. clj.) One begins with the Talmi integral 
expansion for radial harmonic oscillator matrix elements 

/■oo 

{n'l'\Vsr{r)\nl) = / R^,v{r)Vsr{r)Rni{r)r^dr 
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= B{n', I' -,71, 1; p)Ip{b) (218) 
p 

where 

The B{n' ,l';n,l;p) are known analyticalljEij. Here ysr(?') denotes the short- 
range contributions to the potential (we return to this point below). The 
integrals Ip can be viewed as a systematic expansion of the nonperturbative 
hard core in terms of the parameter (rc/6)^, where Tc is a distance associated 
with the size of the hard core that remains unresolved at scale Aqo • Thus the 
LO term in the expansion is 

Q -1 /"OO 

4°io = 7.7^7^/ V^rirydr (220) 



r(3/2) 63 

while in NLO one includes 

r(i) _ ^ W T/ ^.^.2 



2 1 r°°.. _ 



2 1 



In N^LO one obtains 

4% - ^^^l VArydr. (222) 

and so on. 

Now the radial component of the s-wave potential VoS{r) is 

Vs'\r)^Vo^Sir). (223) 



Substituting this into Eq. (218) yields 

{n'l' = 0\v}°\r)\nl ^0) = B(n', 0; n, 0; 0)4(6) (224) 
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where 

It foUows that the contact i ntera ction wiU produce the exact lowest-order Talmi 
integral contribution (Eq. (|220|)) to Vgr provided 



Vb = 47r / Vsrirydr. (226) 
Jo 

Using a normalization that will be convenient later, we can then rewrite 
Vsrir) as 

Vsrir) = Vi'\v) + fe3^;i^<(Aoo)<5(r) (227) 

where vjP is a new potential whose leading-order behavior is of order {rdh)'^ 
relative to Vsr (determined by subtracting from the original Talmi integral 
expression for Vsr the contribution from Eq. ( |220| )). The coefficient ag''(Aoo) 
is a dimensionless coupling (Note that throughout this discussion, we have 
suppressed the spin and isospin of Vsri^)- Thus Vsrir) represents the radial 
function obtained after spin and isospin matrix elements have been taken. For 
example, the components of the AV18 potential contributing to s — s transitions 
is 

Vsr (r) = Vi (r) + V2 (r)fi • fa + X^s {r)Si • ^2 + 14 {r)Si -a^Ti- fa (228) 
so that 

Vsr (r) = Vx (r) - W2 (r ) + ^3 (r) - (r) (229) 

for the deuteron (Z = 0, s = 1, i = 0).) 

The above procedure is, in fact, general. The NLO contribution (Eq. ( p2lD ) 
can be removed from s — s wave deuteron matrix elements of Vsr (r) by intro- 
ducing a second, higher-order contact operator 

«r(Aoo)^(^'<5(r)+^(r)^') (230) 

where 



while the N^LO contributions can be removed by 

af (Aoo)(^'<5(r)^' + ^(^'<5(r) + <5(r)^')) (232) 
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where 



ar(Aoo) - TT-TfT- / V,r{rydr. (233) 
'dnb'Tiuj Jq 

Likewise there ax^coeficients of additional contact operators (denoted a|'*, al'^, 



and af in Ref.lflj) that remove the Talmi integral contributions to s — rf and 



d matrix elements to N^LO. 
Up to this point we have simply rewritten the "bare" potential — the 
potential that acts in the space defined by Aqo — in an entirely equivalent 
form, exploiting the Talmi integral expansion and the fact that the unique 
translationally-invariant two-body contact operators generate those Talmi in- 
tegrals. Now we consider, in LO, the effects on aff of integrating out the single 
shell at Aoo, thereby mapping the original problem into an effective one with 
Aqo — > Aoo — 2. We require that some low-energy matrix element 

{n'0\Vsr{r)\nO) (234) 

remain invariant, so that ag*(A) must run. The relevant equation for ag** then 
is the single-shell limit of the BH equation, evaluated in LO (which eliminates 
sd-BJid dd contributions to transitions and the Green's function). The result 



ar(A - 2) = ag«(A) + -L^ (235) 

(2)! Eo-{A+l)-a^o'{A)^^ 

where Eq = E/hu; is a dimensionless energy. (E is the BH energy being 
determined selfconsistently.) This leading-order RG equation is independent of 
the matrix element chosen to evaluate the running: it is an operator equation, 
independent of the running scheme. Thus, to this order, the most singular 
Talmi integral can be subtracted exactly from all matrix elements in the low- 
momentum space. This difference equation thus tells us how to run ao^(Aoo) 
to any desired new scale A. 

While the RG equation is a difference equation, the connections with EFT 
become even clearer by introducing the natural definition of the derivative 

rf<(A) _ ar(A)-ar(A^2) 

dA - 2 ■ ^^^^^ 

As will be apparent from results given below, the terms in the denominator of 



Eq. (235) that compete with A are one power lower in A. Thus for large A the 



RG equation can be written in a more conventional form 

dag^(A) 1 /A ^,,,,,2 



A-^-2V-<(^) (237) 
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which has the solution 



'(A)~-j|. (238) 



As A is a dimensionless nonrelativistic energy, this is equivalent to running as 
1/p, where p is the momentum, a result familiar from EFT. Similarly, al'* runs 
as A~^/2 and af'' as A~^/^. 



This procedure is followed through order (r/b)'^ in Ref. cl3, which involves 



a series of coupled RG equations for ag'*, al*", af, , al , and a'l . Beycpc 
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lowest order the running is no longer scheme independent; the fit in Ref. 
is to the longest wavelength properties of the low-energy system, the set of 
matrix elements coupling the n = 1 and n = 2 oscillator shells. Thus the 
analogy with standard EFT approaches is very close. In effect, the primary 
difference is that in EFTs the strengths of the local operators are determined 
by directly fitting data. Here a phenomenological potential — itself fit to data 
— determines the initial value of the local operators that are then run by 
the RG equation. It is also apparent that the detailed short-range form of 
the potential is irrelevant to shell-modcl-inspired effective theories: the needed 
quantities are Talmi integrals over the short-range part of the potential. Any 
change in the detailed radial behavior of the potential that preserves the value 
of the Talmi integrals will also le ave l^he low-energy theory unchanged. 



The numerical results in Ref. E13 are quite promising: in N^LO, the hard 
core (below 0.7f) nonperturbative contribution to the scattering can be stepped 
down to A ^ 10 — that is, to nearly a typical SM scale — with a loss of 
accuracy of ~ 7 keV. The portion of the core one can treat accurately increases 
as the order of the RG calculation is increased. As the potential is known, one 
can — if one is willing to invest sufficient effort — carry this procedure out to 
an arbitrary number of orders. The numerical results are given in Fig. One 
sees that the accuracy with which one can run to small A continues to improve 
according to the work invested (LO, NLO, and N^LO) in the RG equations. 

The detailed use of this Talmi integral subtraction and related RG proce- 
dures will not be discussed here, though it is treated in Ref.Ll3. The essential 
idea is to exploit the RG equations to treat the nonperturbative scattering due 
to the hard core, then use perturbation theory to handle the remaining softer 
parts of the potential. Thus Vgr (r) is defined in terms of the Talmi integrals 
for V{r)9{ri. — r), where is the radius of the unresolved hard core. The mar- 
riage of perturbation theory and RG techniques can be done in a number of 
ways: many aspects of this question remain to be explored. However it is clear 
these two tools are complementary. If one starts at very high Aoo, the part of 
the hard core that cannot be treated by perturbation theory is limit to a very 
small region near r ~ 0. However to minimize the effects of omitted higher 
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Figure 58: The results of integrating out the hard-core contribution (vc = 0.7f) to the AV18 
potential by the RG method, in LO, NLO, and N^LO. These orders correspond to the 
removal of successively less less singular contributions of the leading-order Talmi integrals, 
and thus can be viewed as an expansion in (rc/6)^. One can lower the scale of the low-energy 
theory from Aoc = 50 to 10 analytically with the loss of ~ 7 keV in the accuracy of the 
deuteron ground-state energy in N'^LO. 



order terms (i.e., Talmi integrals) in the RG equation for matrix elements with 
large values of n, Tc must be made quite small. As Aqo is lowered, Tc can be 
extended without loss of accuracy, but perturbation theory — an expansion in 
1/A — becomes less effective. 

The extensions to three- and four-body clusters have not yet been carried 
out, though nothing in the procedures is specific to the deuteron. The RG 
equations will, of course, involve the full set of general three- and four-body 
translationally invariant contact operators that can be constructed, just as 
EFT treatments of the three-body problem were required to introduce three- 
body contact terms. Thus the RG equations become more complicated, but 
not conceptually more difficult. 



11 Conclusion 



In this chapter we have given a somewhat detailed overview of the progress 
that has been made in applying the ideas of EFT to nuclear physics. This 
will ultimately provide the bridge between QCD and the rigorous study of 
multi-nucleon processes. The field is still somewhat embryonic, but as we have 
shown, impressive (and rigorous) results have already been obtained. There is 
much more to come! 
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